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SUMMARY 


Methods have been developed for calculating the incompres 
sible flow past airfoil sections, taking full account of the boundary 
laver and wake. It may be hoped that systematic use of these 
will enable quantitative predictions to be made of the influence 
of transition position, Reynolds Number, thickness ratio and 
trailing edge angle on the reduction of coefficients of lift, ete., 
below their values for ideal fluid 

The circulation and the surface pressure distribution are found 
by potential theory. The boundary layer and wake are repre- 
sented in the calculation by a source distribution, whose effect is 
found by transforming the flow to the plane of the complex po- 
tential of ideal flow past the airfoil. An equation expressing the 
constancy of total pressure outside the boundary layer and wake 
then gives the circulation 

With measured boundary layer data for two airfoils of differ- 
ent shapes, the method predicts lift coefficients agreeing closely 
with those observed. A priori predictions of boundary-layer 
growth and thence of lift coefficients, pitching moments, and 
hinge moments on plain flaps of several sizes, have been made for 
a Joukowski airfoil at three Reynolds Numbers. Direct experi 
mental verification of the results has not been possible, but they 


accord well with expectation 


(1) INTRODUCTION 


i he FIND THE WAY in which real fluid flows past an 
airfoil one must proceed by iteration between two 
types of calculation: (I) for boundary layer and wake 
growth with a known external pressure distribution, 
and (II) for the inviscid flow outside the boundary 
layer and wake of an airfoil on which the displacements 
due to these regions are known. This paper formulates 
the latter type of problem for two-dimensional incom- 


pressible flow, and describes a fairly simple analytical 
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method of solution: The subsequent iterative pro 


cedure is also outlined. There are several methods for 
the first type of calculation (including those by the 
present author for turbulent boundary layers! and for 
turbulent wakes”). These leave certain gaps, notably 
about the location of the transition points, but within 
their limitations the full flow pattern may be obtained. 
In particular, sectional characteristics such as the lift 
curve may be predicted without prior knowledge of the 
boundary layer. 

The analysis mainly concerns conditions close to 
the airfoil trailing edge, as is to be expected when one 
considers the relation of flow with small finite viscosity 
to that of inviscid fluid. Neglect of viscosity removes 
from the equations of motion the terms of highest order 
nevertheless 
become predominant Without 
these terms a unique solution does not exist, but when 


which, however small their coefficients, 


near a_ singularity. 
their coefficients are very small their influence on the 
flow may be neglected away from the singularity, pro 
vided conditions in the immediate neighborhood of the 
singularity have been correctly fixed. There is thus an 
arbitrary parameter, most conveniently represented in 
the subsonic case by the circulation, in irrotational in- 
viscid flow The Kutta-Joukowski 
hypothesis that the trailing edge velocity is finite, which 
fixes the circulation uniquely, is thus a suggestion as to 
the influence on the flow pattern of the highest order 


past an airfoil. 


terms when their coefficient tends to zero. 


1.1 Previous Work 

The existence of inviscid circulatory flow with stream- 
lines slightly displaced from those of ideal flow outside 
a finite boundary layer and wake was demonstrated by 
Bryant and Williams‘ in 1926. In the same paper G. I. 
Taylor showed that in a real fluid lift is given by 
the Blasius relation L = pl I of potential theory, 
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provided I’ is the circulation round a large circuit en- 


closing the airfoil and cutting streamlines in the wake 


He suggested the formula, which may 
would be ap- 


at right angles. 
be called the “‘lift-circulation theorem,’ 


proximately true when such a circuit was shrunk to the 


immediate neighborhood of the airfoil. Its exact truth 


under these circumstances has subsequently been 
proved by Preston’ and extended by Temple® to com- 
pressible flows in which disturbances are of order (dis- 
tance from airfoil) ~' at most. 

It follows that equal amounts of vorticity, of oppo- 
site signs, must be shed into the wake from upper and 
lower boundary layers at their separation points. 
This was used by Howarth’ to calculate the lift coeffi- 
cient of a thin elliptic cylinder. The rate of vorticity 
shedding at a separation point S where the free stream 


velocity is q is 


"Ou 1. 
u dy = -q° 
J0 Oy 2 


Howarth therefore fixed the circulation by applying 
the condition g," = qs" on the velocities at the separation 
points 5S), S», these points being obtained by boundary 
layer calculations using the potential flow pressure 
distribution. The final result was thus reached by 
iteration. 

Piercy, Preston, and Whitehead* made an empirical 
allowance by means of sources downstream for a wake 
in Howarth’s method, and Preston’: subsequently 
considered the case of an airfoil whose boundary layers 
are known. He modified the airfoil shape by the addi- 
tion of the displacement thickness 6* on the surface and 
along a line extending to infinity downstream, and at- 
tempted to calculate the potential flow about the new 
body by breaking 6* into symmetric and antisymmetric 


l l 
parts 6,* = > (6.* + 6:*), 6-* = 5 (6.* — 6,*) (suffices u, / 


refer to upper and lower surfaces, respectively). Con- 
tributions from these parts were dealt with separately. 
6.* represents a cambered displacement of the airfoil 
center line, equivalent to a reduction of incidence, and 
6,* can be represented by a symmetrical source dis- 
tribution. Attractive as this approach is, it leads to a 
computation of great complexity, and renders unavoid- 
able some empiricism in determining the final circula- 
tion so as to satisfy the vorticity condition. 

In the present method the boundary layer displace- 
ment effect is represented at the outset by an unsym- 
metrical source distribution, and the circulation which 
satisfies the necessary condition is found in a single step. 
This enables one to demonstrate analytically the mecha- 
nism whereby boundary layers determine the external 
flow, and leads to fairly simple numerical calculations. 
A “change of shape’ method possesses no obvious ad- 
vantage over the immediate introduction of sources in 
incompressible flow, although necessary for the corre- 
sponding supersonic problem, in which the flow is deter- 
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mined by effective shape, and circulation has less rele 
vance. 

The present calculations refer to flows to which the 
ideal fluid theory provides some approximation. Phys 
ically this is the case when the boundary layers do not 
separate far ahead of the trailing edge. The analysis 
is then conveniently carried out in terms of the moder- 
ately small correcting term y = 1 — I°/I'y (where I is 
the actual circulation and I’) that of ideal fluid theory 
The question of flows which involve large separated re- 
gions (e.g., near the stall) is briefly discussed in section 
S. Existing boundary layer theories are inadequate to 
deal with the dead water regions, but no new difficul- 
ties of principle arise in calculating the circulation. 
The forward influence of separated regions on the pres- 
sure distribution is less easily found. The method of 
free streamlines bounding regions of constant pressure 
(Kirchoff flow) seemed in the past to provide the ob- 
vious approach to this type of problem, and a com- 
bination of this with the present type of iterative pro- 
cedure should prove successful, although it presents for- 
midable difficulties. 


(2) List oF SYMBOLS 


0 = suffix referring to ideal flow (Kutta 
Joukowski circulation ) 

1.2 = suffixes referring to points ?;, P» at edges 
of upper and lower boundary lavers at 


trailing edge 7 


t = coordinate in plane of airfoil 
ie = ¢t — b = coordinate referred to airfoil 
trailing edge t = b 
n = from 7 distance measured normal to ideal 
flow streamlines 
§ = lengthwise coordinate along airfoil sur- 
face and wake 
= ae’ + = coordinate in plane of circle 
= 5 — a = coordinate referred to rear stag- 


nation point on circle 
Uy, U = main stream velocity at infinity and at 
edge of boundary layer (not necessarily 
at P,), respectively 


= incidence 


a 

I = circulation round airfoil 

¥ = circulation defect = 1 — ['/T 

Ww complex potential of flow with circula 
tion T 

UW = @ + t¥ = complex potential of ideal 
flow (circulation I’, ) 

q = resultant velocity of potential flow 

Aq = Au — iAv = increment in complex veloc 
ity due to source distribution 

Ww = complex potential of flow due to circula- 
tion + sources 

y* = displacement flux of boundary layer and 
wake = U’6* 

@ = pressure rise coefficient (py — p)/(1/2)pq° 

bi; Go, G—4,.« = coefficients of transformation from / to ¢ 
planes 

Oi, C46 Br = coefficients in intermediate transforma- 
tions via uw and ¢ planes 

Ay Se 5 Mg 3 = constants of airfoil shape in expansion of 
t’ in terms of s’ 

CC De, . = constants derived from A, A’, A” 
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Xx = tan (@ — )r/W1 
C. Cw, Ce = lift coefficient, pitching moment, hinge 


moment 


E = (flap/chord) ratio 

R = Reynolds Number based on chord and 
velocity at infinity = Uyc/v 

T = trailing-edge angle of airfoil 

w = 2-—1r/"r 

K = curvature of ideal flow streamlines 

6, 6*, J = boundary layer thickness, displacement 
thickness, momentum thickness 

H 6*/3 = boundary layer form parameter 


(3) GENERAL THEORY 


3.1 Displacement of Streamlines by Boundary Layer and 

Wake 

Potential theory gives the flow which actually occurs 
outside the boundary layer and wake, provided that at 
the edges of these regions the pressures and streamline 
directions are correctly chosen. The first of these re- 
quirements may ultimately be satisfied by iteration, 
calculating boundary layer growth from the external 
pressure. The second is here met by treating the flow 
as potential everywhere. 
placement outside the region previously occupied by 
boundary layer and wake may then be brought about 


The correct streamline dis- 


by means of a distribution of sources: 
Suppose A is a point of the airfoil surface, B the point 
at the edge of the boundary layer on the normal to the 
surface at A, lL’, u the velocities at height y from the 
surface of potential and boundary-layer flows respec- 
Then the flux across AB in a potential flow is 

*B 
u dy. 


tively. 


2B 
l’ dy, whereas the actual flux is } Thus 


| JA 


e 


for the streamline through B to maintain its position 
when the region previously viscous now forms part 
of the potential flow, an additional displacement flux 


*B 
iin } (U 
1 


at B) must cross AB. 
ducing a source distribution of linear strength dy*/ds 
on the surface of the airfoil, since the nondimensional dis- 
placement flux y*/U pc is small. Similarly in the wake 
dy*/ds correspondingly defined must be introduced 
along the streamline of the real flow which originates 
In this position it experiences 


— u) dy U\6* (where LU’; is the velocity 


This may be secured by intro- 


from the trailing edge. 
no reaction. 


3.2 Criterion for Circulation 


The source distribution having been introduced, the 
potential flow is fully specified by its circulation. This 
must be such as to make pressures at the edge of the 
boundary layer agree with those which actually occur, 
and may be determined from the consideration that 
total pressure p + (1/2)pqg? is constant everywhere 
outside the boundary layer and wake in a steady motion. 
In particular, if suffixes 1, 2 refer to points P;, P, at the 
edges of the upper and lower boundary layers on nor- 
mals from the trailing edge 7, 
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Pi + (1/2)pqi> = pe + (1/2) pg (1) 
On the conventional assumption of boundary-layer 
theory that transverse pressure variations are negligible 


P, = Py = P, 2 


and 


which provides an alternative interpretation of How 
arth’s vorticity condition. In section 4 this is shown to 
lead to a unique circulation when the boundary layer 
thicknesses 7‘P; = m, TP2 = no Clearly, 
however, the criterion throws too much weight on the 


are known. 


positions of P,, P2, which could be varied quite suffi 


ciently to aitect the calculated circulation without 
greatly affecting the neglected static pressure rise or 
the loss in total head if the points were just inside the 
boundary layer. This difficulty is automatically over 
come by calculating coefficients of pressure rise through 


the boundary layer, which may be defined as 


l | 
@ = (P, — P,) . py, We (P, Py) | Pa pq." 
/ =< fz 
(4) 
In terms of these 
gi“/q2” = (1 — @)/(1 — @) = 1+ A@ D) 
where A® = @ — @. A@is small, of the order of 0.01, 


and adjusts the circulation criterion to correspond to 
the positions of P;, P: actually used. 

For convenience these positions will be taken as the 
points at the edges of the boundary layers along the 
normal through the trailing edge to the streamlines of 
ideal flow (i.e., of flow with Kutta-Joukowski circula- 
tion). The pressure rise along this normal may be 
calculated by a method given in section 5. m, m2 are 
used to denote distances measured along the normal. 


3.3 Position of the Wake 


The wake may be thought of as an infinitely long 
flexible tail attached to the airfoil, and the source dis- 
tribution representing it must be placed on either side 
of the streamline which leaves the trailing edge. The 
position of this streamline presents a difficulty in prin 
ciple, since there is no possibility of computing it once 
the full flow pattern has been simplified to a potential 
flow. But it will be seen in section 4.5 that the ve 
locity increments at P,, P, due to sources in the wake 
receive negligible contributions from points more than 
a few boundary-layer thicknesses up- or downstream of 
T. In this distance it is thought justifiable to identify 
the streamline from the trailing edge with the bisector 
of the trailing-edge angle, which is the streamline in 


ideal flow. 
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(4) SOLUTION AS A POTENTIAL FLOW 


In this section conformal transformation is used to 
calculate the velocities at points P;, P, at the edges of 
upper and lower boundary layers at the trailing edge due 
to the source distribution which represents the bound- 
ary layer and wake. The complete potential flow pat- 
tern is then obtained by superposing a circulatory flow 
round the airfoil satisfying the condition (g;/q2)* = 
1 + Aw. The airfoil shape and boundary layer quan- 
tities y*(s), Aa, 7P; = m, and 7P, = ny are for the 
present all assumed known. 


Since the sources are distributed on the surface and 
along the rear stagnation streamline of potential flow 
with the Kutta-Joukowski circulation, a convenient 
transformation is provided by the complex velocity 
potential, w) say, of that flow. In the wy plane the 
configuration becomes a semi-infinite slit along the real 
axis, with a source distribution on either side. The 
transformation scheme is shown diagrammatically in 
Fig. 1, using / for the coordinate in the airfoil plane, and 
z= ae in the plane of the circle through which 


the transformation must proceed. 


Suppose w(/) is the velocity potential of flow round the 
airfoil with an arbitrary circulation I’, I(t) the ve- 
locity potential of this flow combined with that due to 
sources, Agi, Ags the complex velocity components in 
the w, plane due to the source distribution alone. 
(Suffixes 1, 2 are used throughout to indicate evaluation 
at P,, Ps, respectively, or, in the z and w planes, at the 


points which correspond to P\, P2:.) Then 
div dw dw 
= + Aq 
dt, dt, dt, 
dwy ( dw ) - ‘ 
= - Ae (ba) 
dt, dwy 1 ‘ 


Here the term (dw/dwy), depends linearly on I’; the re- 
maining terms are independent of I. Similarly in the 


equation 


div : dwy (=) a 
dts 7 dts Aw 2 " 


Hence I’ may be found 


(6b) 


only (dw/dw»)2 depends on I. 
from Eq. (5), rewritten as 
dW |? dW * 
(1 + Aw) (7) 
dt dts 
This equation is a quadratic in IT, with one root, 
which must be reiected, corresponding to the flow in 
which the velocity vectors at P;, P, are in approximately 
opposite directions. Discussion of the terms may be 
simplified by considering the orders of magnitude in- 
volved. If boundary-layer separation does not occur 
on the upper surface, the incidence a and the defect 
in the circulation below that of ideal fluid theory, 
These 


namely 7 | I'/I, are moderately small. 


may be treated as quantities of the first order of small- 
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nondimensional _ trailing-edge 


ness, ~0.2 say; the 


boundary-layer thicknesses 1,/c, m2/c, a typical value 


for which is 0.04, are then of the second order. 


4.1 Flow in Plane of Circle 


Calculation of flow in the ¢ plane of the airfoil must 
proceed via a circle in the z plane, to which any airfoil 
by the Theodorsen-Garrick 
The complex velocity potential for flow 


may be transformed 
method.!! 
with circulation I past a circle |z 
locity at infinity in a direction making an angle a with 


= a having unit ve 


the real axis is (see, e.g., reference 11) 


4+ (a?/z)e* — (iI'/27) log z (8 


WwW = 2 


If the point z = a, corresponding to the trailing edge of 
the airfoil, is a stagnation point, I has the Kutta- 
Joukowski value 

I) = —47a sin a (9 


and the complex potential is 


WwW = ze ‘* + (a*/s)e* + 2asinalogs (10 
Writing z = ae’ * ”, w = o + ty, Eq. (10) becomes 
o/2a = cosh dX cos (@ — a) — Osina (11 

y/2a = sinh X sin (@ — a) + A sina (12 

y = constant gives the streamlines, @¢ = constant the 


equipotentials. Atz = a,\ = 6 = Oand ¢/2a = cosa. 
Thus normals to the streamlines at the trailing edge 


are given by 
cos a = cosh \ cos (@ — a) — @sina (13) 
Expanding in powers of i, 6 


(A? — 67) cos a + (1/3) (3A? — 67)8 sin a + 
higher powers = 0 (14) 


For small \, 6, the second term is small compared with 
the first, especially since a is small also. Neglecting 


this term, the normals are given by 
@= +X (15 


in the upper and lower half planes, respectively, and do 
not depend on a. It will be seen in the next sub-section 
that the trailing-edge boundary-layer thickness m is of 
the order of \*. Thus J is of the first order of smallness, 
and the last equation is correct to the second order. 
The distances in the wy plane from the trailing edge 


along these normals are obtained from Eq. (12): 


y/2a = +r? + O(A4) (16 
The points in the circle plane corresponding to P;, P: 
» . (i 4 rv l Ar 
are of the form z = ae om, 2 = ae v2 respec 


tively. Postponing to section 4.2 the question of find- 
ing Ay, A» from m, M2, we now write down the velocity 
l 


at the point z = ae ‘ when the circulation is 
(1 — y)I'y. The complex potential is then 
w= ze '* + (a*/z)e* + sina (1 — y)alogz (17 
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and the complex velocity is 
du a Z a — 
= e _ e + (1l— y2smea (18) 
dz z La Z 
lo get some insight into the process whereby the 
boundary-layer controls the external flow, the power 


series expansion 


s/a = 1+ (1+ id + id? + O(A) (19) 


may be substituted in the expression for velocity, and 
the first two terms retained. Coefficients of higher 
powers of \ in the resultant expression rapidly become 
complicated, and to obtain greater accuracy in a nu- 
merical calculation than is given by the formulae be- 
Results 


obtained in this way will be as correct as those from 


low, it is best to use the exact value of z/a. 


Eq. (15), 1.e., to the second order of small quantities. 


Substitution of Eq. (19) into Eq. (18) gives 


u a 
- x 
(1+ 2)2A\ cos a + 2A* sin a — y:2isin a] 
oe i a . 
= (1 + 2) 2A cos 
a—~ x (20) 


In particular, putting y = 0 
= (1+ 7) 2d cosa x 


l 
l + d tan a (1 — 1) | 
a 


1+ 72)2\ cosa X (21 


l 
l1—A (1 ates tan a) — 
‘ l 
wi 1+ tan a) 
Z 


du , 
: = SA- cos" a {1 — X(2 — tan a) + O(A*) | 22) 
az 
and 
dw = 5 ?; 
=]|]- tan a (1 + 2) + ()(A\-) 93) 
dw é / 


4.2 Transformation of Airfoil to Circle 


Every airfoil in the ¢ plane is connected uniquely 
with the circle | = ain the s plane by a transformation 
of the form 


= F(z) = aya +a + a_)/2 + a_2/s? +... (24) 


where F(z) is analytic on and outside the circle 2) = 
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a, except possibly at a, which corresponds to the trail- 
ing edge ¢ = F(a) = bsay. There F’(a) = 0. 

When the transformation is known, the values of \ 

+ . ” 1 d . . 
at the points ae which correspond in the z plane 
to P;, P2, may be found from 

m,, mn. = |F(ae' ***) — b 25) 

The modulus of transformation 


dt 
dz 


= V2 ae |F'(ae" 26) 
may then be obtained directly. 

In the general case of an airfoil of arbitrary shape the 
transformation (24) is not known, although the corre- 
spondence between points of the circle and points of the 
airfoil may be found with any desired accuracy by the 
Theodorsen-Garrick!'! method. Sufficient information 
can however be obtained by considering the form of the 
transformation in the immediate neighborhood of the 
trailing edge. 

Let us write 2 — a = 2',t — 


, ? 


b = t’, and denote the 
trailing-edge angle of the airfoil by r. The transforma- 
tion which takes the trailing-edge region of the airfoil 
surface into the tangents at the trailing edge is of the 
form 


u(t’) = {+ ht’? ae ee Zé 


say. Similarly the transformation which takes the 
circle in the neighborhood of s’ = O into its tangent at 


that point is of the form 


Pie’) = ¢g ~ 15-> IS) 


The relation between the tangent pair in the uw plane and 
the single tangent in the ¢ plane is of the Schwarz- 
Christoffel type 


-- r 20) 


that is, 


(1+ dt’ +... z’ “= (1+ d,s’ + 
30) 
where the series in parentheses are analytic near the 
origin. Inverting the series on the left-hand side of 
Eq. (30), we may write 


where w = 2 — r/r(>1). TheconstantsA,A’,A”,... 
may readily be found from the known correspondence 
between points of the circle and points of the airfoil. 
(If the real axis in the ¢ plane bisects the trailing-edge 
angle, A is purely real. For a symmetrical airfoil, 
Ade ad 

In the following, the symbols B, B’, C, C’, D, D’, 
D",... are used to denote quantities which can be ob- 
tained from the A’s by simple algebra. 
are not carried out here, since only the forms of the re- 


are also real.) 


The derivations 


sults are of interest. 
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If z’ lies on a normal to the streamlines at the trailing 
edge, then 
tA = 


l + 
s'| = ale —-lj= 


l l 
V2an(1+ 54+ a+...) (32) 


Moreover, from Eq. (31) it follows that 


Als’\* [1 + Bis’ 


| + (Q(/z’|°) 


Therefore 
(3.3) 


n = A(V2adx)° [1 + BA + O(A*)] 


7 is usually small, so that ” is approximately of the order 
of d°. 
To obtain the velocity of ideal flow in the airfoil plane 


must be evaluated: 


dt dt’ ‘iy ia I rt tw 
= ~ = Aws”™ lL +-{1+ A's’ + 0(z") 
dz dz w 


Therefore 
dt 


dz 


= Cr°~ '[1 + Cr + O(A’)] (34) 


+ 1)Xd 


> > << ’” I 
From Eq. (22) we have, for s = ae 


Iwo l 
: =2NV 2 cos a} -_ A(1 — > tan «) +004 | 
dz 2 


Therefore at the corresponding point in the airfoil plane 


the velocity is 


dwy 2V2 cosa FX 
dt| C 


l 
[ < (C” + tan a)A + 00%) (35) 


Expressed in terms of distance m in the airfoil plane, 


this is 
dw a/eum , 9 Ses 
i le Dn” *"[{1 + (D’ + D” tan a)n’ * + O(n)] 
cf 
(36) 
Il — £)A 


The expression for velocity at a point z/a = e 
on the lower normal is obtained by changing the sign 
of tan a in this expression, D’, D” being unaltered pro- 
vided the airfoil is symmetrical. 


4.3 Velocity in w, Plane Due to Source Distribution 

The correspondence between the surface of the air- 
foil and the real axis in the wp plane is obtained via the 
z plane, in which the stagnation streamline behind the 
trailing edge is given by 


6 = a — sin~'(X sin a/sinh \) (37) 


from Eq. (12), since Y = 0. From this relation the co- 
ordinates of the stagnation streamline in the airfoil plane 
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may be found using Eq. (31), and thus the distance s 
along the streamline in terms of A. Substitution of 6 
from Eq. (37) into Eq. (11) gives ¢ in terms of A, and 
hence of s. The airfoil surface corresponds to \ = 0, on 
which the relation is found more simply. 

P, corresponds to the point ¢ = 2a cos a = @r (say), 
y = 2a\," cos a = yY. The complex velocity at this 
point due to a source of strength m at a point (@¢, 0) of 
the real axis (which must be regarded as a solid bound- 


ary) is 
— iy 


wT (dr — o)* + he 


m (dr — @ 


Writing 
dy* ds dy* 


m = = 
ds do do 


and integrating along the real axis, the velocity at P, is 
Aq = An, = 1Av, 
i dy* or—o— if 
d do (¢r — o)? + Ww 


(38 
us 


A more useful form is obtained by substituting @ = 
gr + y tan x: 
After a little reduction 


1 pre dy* 
An, — iAr, = f a 
TH J (2/2) Ix 


I : 
( = sin 2x + 1 costx dx 


Integrating by parts 


l 
x 
wy 


Bos 
lv = sin 2x + 7 cos? x)| sis 


. fe as 
| y* (cos 2x — isin 2x)dx (39 
Ty, ) 


The first term vanishes at both limits, and we are left 


Au, — iAvy, = — 


with 
1 r/2 
Am, = y* cos 2x dx 
TY —(r/2 
(40) 
l w/e ; 
Av, = y* sin 2x dx 
mY —(r/2 
Thus differentiation of y* may be avoided. Am, 


Av, are clearly determined by the values of y* in the 
immediate neighborhood of the trailing edge, for 90 
per cent of the integral is contributed by the range 
o — dr < dW, ie., in a chordwise distance of the 
order of 5 boundary-layer thicknesses. 

From the form of Eq. (40) it may be concluded that 
Aq depends on the shape rather than on the size of the 
distribution of y*. y*/y, is a measure of the ratio 
5*/n, of displacement thickness at a point of the air- 
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foil surface to total boundary-layer thickness at the 
trailing edge; thickening or thinning of the whole 
boundary-layer does not affect the order of this ratio, 
but if accompanied by change of displacement thick- 
ness rate of growth and of velocity profile shape at the 
trailing edge, as would normally happen, it changes the 
value. A typical distribution of y* is shown in Fig. 1. 


4.4 Equation for Circulation Defect 


Eq. (7) may now be used to calculate the circulation 
defect y using the power series expansions derived in 
sections 4.1 and 4.2 for the various terms. Only quan- 
tities of the first order of smallness will be retained. 
To this order the expressions for quantities at P; be- 
come, from Eqs. (23), (35) 


dw * S cos* a 
dt 7 C? 
\-"’" [1 + (C” + tan a)A] 
dw 7 ‘ 
+ Ag = 1 — tan a(l +12) + 
dw ~A\ 
Au; — 1Av, 
dw . 7 
F + Aq} = 1+ 2Au — tana (41) 
adWo “1 
Therefore 
dW ? S cos” a —* _ 
= = AL 1+ (C° + tan a)d;] x 
dty ( i 


[ + 2Au, — ~ tan a (42) 
Ay 


The corresponding expression for P: is obtained by writ- 
ing A. for Ay and —a@for a. (For an asymmetrical air- 
foil C” must also be altered, but for simplicity this is not 


done here): 


dW ? S cos” a i a 

- = —— \g””’* [1 + (C” — tan a)A2] X 

dts Se 

[ + 2Au, + — tan a (4:3) 
Ay 
; ; dW |? |dW?. ; 7 " 
Equating ; in accordance with Eq. (7), 
~ | dt dts 


(“)" 
1— Aa) 
No 


[1 + Ch —_ A2) + tana (Ay + da) | = 


l l 
1 — 2(Au, — Auz) + ¥ tan «( T \ ) (44) 
Al A2 


The equation may be written 


(“)’ 7 
(1 + e) — 
rs 


1 + 2(Au,; — Aus) 


an equation for y. 


] l 
ytan a ( T ) 
Ay A2 


(45) 


where 
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(and clearly, « < 1). 

This equation provides a theoretical explanation of 
certain qualitative boundary-layer effects on lift which 
are well recognized experimentally. For an airfoil of 
finite trailing edge angle 7, the effect of thickening the 
upper boundary layer, keeping the lower constant, 
The 
to decrease the 
circulation and the lift, as expected. The 
effect should occur if \; were decreased (e.g., 


would be to increase \; and probably also Au. 
effect would thus be to increase jy, 1.e., 
converse 
by in- 
creasing the amount of laminar flow on the upper 
surface). The effects of 
A», Which represents the lower boundary-layer thick- 


increasing and decreasing 
ness, are not clear from the equation, and would appear 
to depend on actual rather than on relative thicknesses. 
It might be expected that the larger the trailing-edge 
angle, the more marked will be such changes in lift ob- 
tainable via trailing-edge thickness (e.g., from transi- 
tion movements or changes of Reynolds Number). 

Joukowski's hypothesis, in the form y —~ 0, may be 
deduced from Eq. (45) if one assumes that Ao, Am, 
Au. and the ratio \;/A2 remain finite as the viscosity 
and therefore boundary-layer thickness tend to zero. 
For then, since A;, A» > 0 as vy > 0, the coefficient of y on 
the left-hand side becomes infinite, while the right-hand 
side remains finite. 
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In streamline coordinates (s, 7) the relation 


1 Op 


Kwu= —- 
p On 


(where A, is the curvature of a streamline) expresses 
the balance between pressure gradient and centrifugal 
force in a curved boundary layer. When 7 is measured 
along the trailing-edge normals as defined in section 
9 


) 
eds 


1 Og 
q on 


i. = 


where g is given by Eq. (36). Thus 


T l Tv , ” 1 
iS — =)" tana) 


OT n ¥ 


te] 


(46) 


upper and lower signs corresponding to upper and 
lower normals, respectively. 

To estimate the pressure rise through the upper 
layer, the power law profile (see, e.g., reference |) 

‘ 4 \ (1/2 = 1 

U; ( ny ) 
(where m; = boundary-layer thickness) may be used. 
(The line along which m is measured, despite its in- 
finite curvature at the surface, is sufficiently close to a 
normal to the mean flow in the boundary layer to justify 


the use of this interpolation.) We have then 


oe n H l 
Pa A: t= 2f x,(“) dn 
l 0 ny 


g Ur 
= = UT, — 1)! + 
WT 
» l l 
n° (D’ + D" tan a) (11 ~ - .) (47) 
Ww Ww 


Similarly for the lower layer 


9 


ji< — Gh = 2)" + 
WT 
1 1 
ns" * (D’ — D" tan a) (1. + = 1) (48) 
w w 


From the first term of each expansion it would ap- 
pear that pressure variation is greater when a boundary 
layer is remote from separation at the trailing edge 
(i.e., when /7 there is low) than when separation is im- 
minent. This prediction, at first surprising, expresses 
the fact that in the former case the centrifugal force on 
the boundary layer is greater because of the higher 
velocities close to the wall where curvature is greatest. 
The prediction is modified somewhat by the second 


term, which is larger for a thicker layer. However, it 
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may be concluded that in general pressure rise is greater 
through the lower than through the upper layer, and 
A® = @ — @ 1s negative, so that (g,/q2)* is slightly less 
than unity (from Eq. (5)). The relative importance 
of the first term decreases with trailing-edge angle -. 
and vanishes for a cusped airfoil (r = 0, w = 2). In 


that case 
o= 2 


:(D' + D"” tan a) [H — (1/2)]-" (49 


and the pressure rise through the upper boundary 
layer might be expected to exceed that through the 


lower. (q:/g2)” would then be slightly greater than 
unity. 
Static pressure in a curved stream is difficult to 


measure, and the predictions of this calculation cannot 
be tested conclusively. However, the experimental 
data discussed in the next section suggests that the cal- 
culation is satisfactory. Qualitatively, the relative 
pressure rises are as predicted for the cusped and finite 
angle trailing edges of the airfoils referred to; and it 
may well be that calculation is more accurate than 
measurement since it does not disturb the flow. 


(6) CALCULATIONS USING MEASURED BOUNDARY 
LAYERS 


At the outset it was desirable to test the accuracy of 
the method by computing the lift for an airfoil on which 
the boundary layers had been measured. The data 
obtained in Preston’s experimental investigation! 
was suitable tor this purpose, and convenient in pro- 
viding a comparison between the present method and 
his theoretical work. Boundary-layer growth and lift 
coefficient were measured on a symmetrical Joukowski 
airfoil 11.8 per cent thick and on a Piercy 12-40 airfoil, 
both at 6 deg. incidence at a Reynolds Number of ap- 
proximately 0.5 X 10°. These airfoils were chosen as 
having approximately the same thickness ratio’ but 
widely different trailing-edge angles, 0 deg. and 22.10 
deg., respectively. 

For the Joukowski airfoil, the measured lift coefh- 
cient was 0.88 C;,. The measured (q:/q2)* was 1.025. 
The present method, using this ratio, predicts C, = 


0.91 C,,. Preston deduced the same value. However 
from the boundary-layer thicknesses the present method 
predicts (g;/g2)? = 1.031, on which basis C;, = 0.92 C7z,. 


The difference between experimental and measured 
C,’s may be due in part to the difficulty of constructing a 
cusped trailing edge. An imperfect cusp would account 
for the observed loss of lift due to the boundary layer 
being greater than that predicted. Similarly, the 
measured A® = 0.025 compared with the predicted 
value of 0.031 would be explained by finite trailing-edge 
thickness, but it is in any case doubtful whether better 
agreement could be expected from measurements of 
Ao. 

For the Piercy airfoil the present method predicts 
30.6 per cent loss of lift, which is in exact agreement 





wit 
cul 
wa 
be 

me 
circ 
—() 
whi 
P,, 

ure 
per 
valt 


iS 4 
kno 
lift 


“I 
— 


T 
give 
tion 
(i 
the 
Ui; 4 
cula 
on t! 
used 
velo 
faces 
0.9L 
(ii 
wake 
The 
to c 
the k 
valu 
the ¢ 
edge 
powe 
y* ii 
whic! 
valid 
point 
neces 
(iti 
Stagn 
the r 
the ré 
Au, . 
on th 
(iv 
edge 
(v) 
tainec 


2reater 
Tr, and 
tly less 
rtance 


igle A 


indary 
rh the 
than 


ult to 
‘annot 
nental 
1e cal- 
‘lative 
finite 
ind it 
than 


RY 


icy of 
which 
data 
) i? 1 

pro- 
1 and 
d lift 
ywski 
irfoil, 
f ap- 
en as 
but 
22.10 


oeffi- 
025. 


ever 
thod 
4 o = 
ured 
ing a 
ount 
ayer 

the 
icted 
edge 
-tter 


s ol 


licts 
nent 





CHARACTERISTICS OF 


with observation. This prediction is based on a cal- 
—0.040. The observed Aw, however, 
It is uncertain how much weight should 


culated Ao = 
was —0.005. 
be attached to this value; it is not consistent with the 
measured boundary-layer thicknesses and measured 
circulation, since these must inevitably predict Ao = 
—().040. Moreover it is unlikely that the points at 
which the velocities were measured were exactly the 
P,, P» to which the calculations refer. Using the meas- 
ured ratio (g,/q2)*" = 0.995 leads to a prediction of 27 
per cent loss in lift. With the same data Preston's 
value was 27.5 per cent. 

From this agreement it is concluded that the method 
js accurate when the boundary-layer thicknesses are 
known, and that the accuracy of a priori predictions of 
lift depends only on that of boundary-layer growth. 


(7) A Priort CALCULATION OF FLOW PATTERN 
7.1 Iterative Procedure 

The iterative procedure for calculating the flow at 
given incidence a, Reynolds Number U c/v and transi- 
tion positions, is as follows: 

(i) With an assumed circulation, say (1 — yo) times 
the Kutta-Joukowski value, the potential flow velocity 
U, at the edge of the boundary laver and wake is cal- 
culated. Except close to the trailing edge, the velocity 
on the surface and on the stagnation streamline may be 
used; at the trailing edge there is a singularity, and the 
velocity at the edge of the boundary layer for both sur- 
faces must be faired to a finite value of the order of 
0.9L’. 

(ii) The growths of 3, 77 in the boundary layer and 
wake are calculated from this velocity distribution. 
The thickness at the trailing edge may be used at once 
to calculate the potential flow velocity at the edge of 
the boundary laver and if necessary to correct the faired 
value. The boundary-layer quantities required are 
the distribution of y* = U\6* = U,8H and the trailing- 
edge thicknesses 1), 22 which may be obtained from the 
power law relation 6/3 = J7(/7 + 1)/(/7— 1). Strictly, 


v* is given not by U,é6* but by | (UW) — u)dy 


0 

which is slightly different at the trailing edge; but the 
validity of the form parameter equation for // at such a 
point is itself uncertain, and a correction would not 
necessarily improve the estimation of y*. 

(iii) The trailing-edge region of the airfoil and the 
stagnation streamline of ideal flow is transformed to 
the real axis in the plane of w) = @ + ty by means of 
the relations of section 4.3, and the velocity components 
Au;, Aus in the @ direction due to a source distribution 
on the axis calculated by Eq. (40). 

(iv) Coefficients of pressure rise through the trailing- 
edge boundary layer, and hence the difference Ao = 
® — @ are calculated by the method of section 5. 

(v) Using the values of Ax, Ae, Am, Aus, A® now ob- 


tained, Eq. (7) is solved for y. Calling the value so ob- 
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tained 7, the calculation must then be repeated until 
convergence is achieved, i.e., until y, = Yn, Say. 
The following qualitative argument based on experience 
in using the method shows that if yo is an overestimate, 
y: will be an underestimate of the true y, and vice 
versa: The value calculated from Eq. (7) depends 
mainly on the trailing-edge boundary-layer thicknesses 
represented by i, Ae; variation of circulation affects 
the pressure distribution and hence boundary-layer 
growth on the upper surface much more than on the 
lower surface, and thus broadly speaking affects \, but 
not Av. 
on the upper surface will be less severe, and therefore 


If yo is an overestimate, the pressure gradient 


the boundary-layer thickness and \,; smaller, than is 
actually the case; if A, is smaller than its true value, 
Eq. (45) leads to underestimation of y. 

This suggests an interpolation which shortens the 
iteration: if starting from yo leads to y;, and starting 
from y; leads to yo, then starting from Ay) + (1 — 
A)y: might be expected to lead to Ay; + (1 — A)y; 


choosing A to make these equal, ie, A = (ys — y;) 
(yo — 271 + Y2), gives the value y = (yoyo — yi? 
(yo — 2y¥1 + y2) as an estimate of the true y. 


7.2 Lift Coefficient and Sectional Characteristics ° 


(i) The lift coefficient is given at once by the lift- 
circulation theorem ZL = pliI’. Thus C, = (1 — 
y)C,, where C,, is the lift coefficient in ideal fluid. 

(ii) The pressure distribution on the airfoil surface 
can be found by a fairing method similar to that used 
for the velocity at the edge of the boundary layer: 
having found y we know the velocity at P; exactly, and 
therefore the static pressure p,; hence, from the pressure 
rise coefficient w,, the pressure at the trailing edge. 
The theoretical pressure distribution on the airfoil sur- 
face for potential flow with circulation (1 — y)I'y may 
thus be faired to the correct trailing-edge value. It 
may be noted here that Pinkerton'* has found experi- 
mentally that the surface pressure on an airfoil agrees 
with that of potential flow with the circulation inferred 
from the observed lift coefficient, except near the trail- 
ing edge. 

Strictly, the fairing process could be replaced by 
computation, as all the quantities involved are now 
known, but increased accuracy would be unlikely. 

From the pressure distribution, pitching and hinge 
moments may be calculated by quadratures. 


7.3 Predicted Characteristics for a Joukowski Airfoil 


Calculations have been carried out for the lift co- 
efficient and pitching moment, and for the hinge mo- 
ments on plain flaps, for a simple Joukowski airfoil 11.8 
per cent thick at 6° incidence at three Reynolds Num- 
bers, 10°, 4 X 108 and 4 X 10%. This airfoil was 
chosen because of the simplicity of the transformation 
which relates it to the circle. 





586 JOURNAL OF THE AERONAUTICAL SCIENCES—SEPTEMBER, 1954 


= I 



















forw 
forn 
by t 
4 me to te 
ind i 
| Tl 
2 f a first 
3 # the I 
itera 
pr lat 
2 = ~ favol 
ok TI 
three 
a al Jouk 
wer SRS ome 
\ < f com 
cient 
Or / # : “ and | 
/ R=4x1O per ¢€ 
/ R=4x10° coeffi 
“a a R=0o 4 = culat: 
(JOUKOWSKI os 
CIRCULATION ) \. that \ 
a 2 p _ 0, - 
’2 Us 
a g . J , 8.1 Be 
Svs 
be a) 


4b (Y , istics 

x specif 
S VARIATION OF SURFACE PRESSURE. | J sist 
i & WITH REYNOLDS NUMBER ON 7 |“, 
S and st 


ms A _JOUKOWSKI _AEROFOIL ] | 


be de 











transi 


as wel 


-7e —— SURFACE PRESSURE IN IDEAL FLOW 1 iio. 


WITH ACTUAL CIRCULATION $2 Se 
~§r Cak 











separa 

———— SURFACE PRESSURE ALLOWING FOR the tra 

<< e on able o: 
9 x BOUNDARY LAYER. investi 
Wife displac 

1.0 l | l l 1 | | a nag 
2 3 5 6 7 8 9 O _ 

Fic. 2. region. 


ent me 
the nei 

















CHARACTERISTICS OF 

Laminar boundary layers were calculated from the 
forward stagnation point by means of the Thwaites™ 
formula 3? = 0.45v U-6 , i L* ds, and turbulent layers 
by the method of reference 1. Transition was assumed 
to take place at the suction peak on the upper surface, 
and at 0.4c¢ on the lower. 

The circulation derived from the boundary layer 
first calculated was used to provide a second estimate of 
the boundary layer and thus of the circulation; further 
iteration was then avoided by using the linear inter- 
polation outlined in section 7.1, a trial in the most un- 
favorable case having shown this to be very accurate. 

The calculated surface pressure distributions for the 
three Reynolds Numbers and for ideal fluid (Kutta- 
Joukowski circulation), which corresponds to R = &, 
are shown in Fig. 2. Fig. 3 shows the ratio between the 
computed and ideal fluid theory values of the coeffi- 
cients of lift, pitching moment about the leading edge, 
and hinge moments on plain flaps of 20 per cent, 30 
per cent, and 40 per cent flap,/chord ratio. The lift 
coeflicient is obtained directly from the calculated cir- 
culation, the remaining coefficients by integration of the 
computed chordwise loading. As far as can be judged, 
the results are correct qualitatively, and of the order 
that would be expected. 


(8) Discussion 

8.1 Boundary-Layer Transition 

Systematic use of this type of calculation appears to 
be a perfectly feasible means of finding the character- 
istics of an unstalled airfoil section with transition at 
specified points (either fixed, or related to the pressure 
distribution in a known way, e.g., at the suction peak 
or at the laminar separation point). Inability to pre- 
dict transition conditions accurately is perhaps the 
most serious remaining gap in boundary layer theory, 
and such phenomena as laminar separation followed by 
turbulent reattachment—the so-called ‘‘bubble’’ transi- 
tion found on thin airfoils approaching the stall—cannot 
be dealt with at present. This particular type of 
transition has a local effect on the pressure distribution 
as well as influencing it via the trailing-edge boundary 
layer. 
8.2 Separated Flow 


Calculations for incidences high enough to produce 
separation of the upper boundary layer well ahead of 
the trailing edge would require knowledge not yet avail- 
able of the flow in the separated region. The present 
investigation has shown, however, that the streamline 
displacement due to the boundary layer is not the 
most important factor in producing loss of lift, and mod- 
erate accuracy could be expected from a calculation 
which ignored the displacement due to the separated 
region. Thus a combination of Howarth’s and the pres- 
ent method could be used to obtain lift coefficients in 
the neighborhood of C_,,,,., circulation being determined 
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from a condition g,” = g.*, where q, is now the potential 
flow velocity at the separation point, and g, that at the 
edge of the lower boundary layer at the trailing edge. 
Such a calculation would be carried out directly in 
terms of I, without any reference to ideal flow condi- 
tions. In calculating the boundary layer however, 
one must have some idea of the influence of the sepa- 
rated region on the upstream pressure distribution, and 
this could only be got from the type of mixed calcula- 
tion, making some use of free streamline methods, 


which was suggested in the ‘“‘Introduction.”’ 


8.3 Uniqueness of Solution 


The argument of section 7.1 may be used to infer 
that even when separation occurs ahead of the trailing 
edge only one C, is possible for a given type of transi- 
tion, since the iteration process must converge. For 
if Cy, is overestimated, the adverse gradient leading to 
separation will be too severe, and separation will occur 
earlier with consequent reduction in lift. However at 
a given (high) incidence, two C,’s may be possible 
because of different types of transition before and after 


a wing has been stalled. 


8.4 Compressibility Effects 

For thin airfoils at small incidences at Mach Num- 
bers less than about 0.7 the use of the Prandtl-Glauert 
Law should be justified, characteristics computed by the 
present method being multiplied by (1 — 47*%)~' *. 
Data regarding compressible turbulent boundary layers 
are rather limited, but it seems likely that at these 
Mach Numbers, which are well below those at which 
shock waves appear on a thin airfoil, the effects of com- 
pressibility in the boundary layer are slight. 

At transonic and supersonic speeds, however, the 
problem becomes one of change of effective airfoil 
shape complicated by the presence of shock waves, 


(Continued on page 620) 








The Secondary Flow About Struts and 
Airfoils’ 


W. R. HAWTHORNE? 


Cambridge University 


SUMMARY 


The flow about a strut when the approaching velocity varies 
in the spanwise direction has been examined theoretically and 
experimentally. The theory supposes that the fluid is inviscid 
and incompressible and that the flow is a small disturbance of 
the two-dimensional potential pattern appropriate to the profile. 
It is found that the shape of the profile at the nose influences 
the size of A rounded 
causes such large disturbances that the original assumptions of 


the disturbances appreciably nose 


The disturbances are much less with 
A method 


the theory are invalidated. 
a wedge-shaped nose and least with a cusped nose. 
of designing profiles to minimize the disturbances is described 
and an example The flow about 
the bicusped strut is examined in detail and the energy in the 
secondary flow normal to the direction of the approaching flow 
is determined for one type of initial velocity distribution. This 
energy is found to vary as the fourth power of the thickness of 
the strut and to be a maximum at a certain value of the displace- 
ment thickness of the original nonuniform stream. The second- 
ary flows produced by a biconvex profile are shown to be small. 
Calculations for a thin circular are airfoil are presented to show the 
effect of lift and the existence of vortex sheet trailing off such an 


a bicusped strut—is given 


airfoil is indicated. Ling’s experiments with a bicusped and con- 
ventional profile placed in a nonuniform stream are summarized 
to show their confirmation of the major effects predicted by the 
theory. The drag due to secondary effects is appreciably less 
with the bicusped strut and so is the general disturbance of the 
flow as shown by wake traverses. The energy in the secondary 
flow is of the order of that predicted by the theory (for the bi- 
cusped strut) but the drag is much greater than would be de- 
duced from measured or calculated energy. 

Some experiments with struts placed in an hydraulic channel 
with a sandy bottom show that the scouring of the beds of rivers 
may be attributed to secondary flows. The use of a bicusped 
profile reduces the scouring considerably compared to an elliptic 


profile of the same thickness and chord. 


INTRODUCTION 


y I (ig FLOW AROUND airfoils and struts placed in an 
approaching stream with a nonuniform velocity is 
an instance of three-dimensional rotational flow which 
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occurs frequently. An example is the flow at the 
junction of wings and tail planes with an aircraft fuse- 
lage, where the wing is swept by the nonuniform flow 
in the boundary layer on the fuselage. Another example 
is in axial compressors and turbines, where the boundary 
layers on the inner and outer walls of the annulus pass 
through rows of stationary and rotating blades. A 
bridge pier in a river in which the velocity varies from 
top to bottom is a further example of a strut in a non- 
uniform flow. 

The flow around a strut projecting from a wall along 
which fluid is flowing is complicated. Fig. 1 shows this 
flow diagrammatically. Far upstream the pressure is 
uniform, but the approaching velocity U varies in the 
direction of the span of the strut. The disturbance to 
the flow may be thought of as due to the variation in 
stagnation pressure at the leading edge, which causes 
a downward spanwise velocity. This spanwise flow 
converts the comparatively simple upstream flow into 
a complex three-dimensional flow, with velocities nor- 
mal to the main flow which persist downstream of the 
strut. These secondary velocities bear some resem- 
blance to the velocities induced by the trailing vorticity 
from a finite wing. 

The analysis of this type of flow has been discussed 
by von Karman and Tsien,! Squire and Winter,’ Haw- 
thorne,* and Kronauer.‘ The effect of viscosity has 
been neglected, except for its influence in establishing 
the nonuniform initial flow. von Karman and Tsien 
have discussed the behavior of a wing represented by a 
lifting line in a nonuniform flow. Squire and Winter 
analyzed nonuniform flow in a curved passage when the 
disturbance of the initial flow was small, and showed 
the important growth of the component of vorticity 
in the direction of flow. Hawthorne derived a general 
expression for this component of vorticity and showed 
that it could be used to predict approximately the char- 
acter of the flow in bends in which the original motion 
was greatly distorted. Kronauer analyzed the flow 
through cascades of airfoils and struts and discussed 
the behavior of the flow at the nose of the airfoils and 
in the cascade passages. While all the experimental 
work of these authors shows that the distortion of the 
flow tends to be considerable, some success has been 
achieved with the inviscid theory by regarding the 
three-dimensional flow as a perturbation of the initial 
quasi two-dimensional flow approaching the body or 


bend. 
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SECONDARY FLOW 


In this paper the results of this perturbation theory 
are derived and examined to determine their limitations, 
and the theory is used to evolve strut profiles which 
minimize the distortion of the initial flow. 

The paper is divided into three parts. Part I con- 
tains a derivation of the equations for the components 
of vorticity in the flow. This repeats Squire and 
Winter’s analysis in a slightly different form. 
sions for the total secondary circulation induced in the 


Expres- 


flow are then derived and used to show the influence of 
the shape of the leading edge. 

In Part II the results of applying the theory to bi- 
The 


flow about a thin circular arc airfoil is also examined to 


cusped and biconvex strut shapes are discussed. 


show the effects when lift is present. 

Part III contains the results of some experiments in 
which a bicusped strut is compared with more con- 
ventional shapes. 

Part I—GENERAL THEORY BASED ON THE SMALL 
PERTURBATION OF A TWO-DIMENSIONAL FLOW PATTERN 


(1) A Derivation of the Equations for the Components of 

Vorticity in the Flow+ 

Far upstream of the body, let the velocity be U in 
the x direction (Fig. 1) varying only in the z direction. 
The basic velocity will be that corresponding to the 
two-dimensional flow around the body in an xy plane 
If V is the vec- 
tor of this two-dimensional motion when the far up- 
stream velocity is unity, then V will be invariant with 
cand the vector UV will describe the basic quasi two- 
To this basic flow will be added a 


in whic the approach velocity is U. 


dimensional flow. 
small velocity whose vector is v so that the velocity 
anywhere in the field is UV + v. 


t In Squire and Winter’s paper the following analysis is pre- 
sented by writing down the equations of motion in curvilinear 
coordinates, two of the coordinates being the streamlines and po- 
tential lines of the flow. The longer 
method adopted here is offered as an alternative and is directed 
toward obtaining a particular form of the result which is useful 


basic two-dimensional 


| inthe subsequent analysis 





The vector 


U’ grad gq? X grad U = sU(0q°/O0n) (OL’/Os) — nU(0g*/ds) (OU/dz) 


The vector 


(@-V)UV = [&(0/ds) + n(0/On) + £(0/0z)]sUq = s[é 


ince $ and g are invariant with z and U with s and n. 
vector V along the streamline and hence 


Os /Os 


ABOU 


n(06/Os) 
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The flow is assumed to be inviscid and incompressible. 
The equations of motion will be solved for the pertur- 
bation vorticity @ = curl v by noting that in inviscid, 
steady flow the cross product of velocity and vorticity 
is equal to the gradient of a scalar (1.e., the total pres- 
sure divided by the density) and therefore the curl of 
this cross product must equate to zero. 

Since V is a vector of a two-dimensional potential 
motion, div V = 0 and curl V = 0. 

From continuity, for the three-dimensional flow, 


div (UV + v) = 0 
or LU div V+ V-grad U + divuv = 0. Then div v 
0. Nowcurl (UV + v) U curl V + grad l kK VG 
curlv. Hence 
curl (UV + v) = grad U X¥ V+ 1.1 


Now the cross product of velocity and vorticity is given 


by 


(UV + v) X (grad U X V+) = 
Uq* grad U + (v-V) grad lL’ — 
(v-grad LU’ V+UVXet+vXeo 


where g is the scalar of V. 

If v is assumed to be small compared with lg, the 
second and third terms on the right-hand side of Eq. 
(1.2) are small compared with the first, and the fifth is 
small compared to the fourth. Hence, only the first 
and fourth terms need be retained. 

Taking the curl of the remaining terms in Eq. (1.2), 
simplifying, and noting that div v and div @ are zero, 


U grad g* X grad U + (@-V)UV — 


(UV: Vio = 0 (1.3) 


Expressions for the components of w will be obtained in 
three directions, in the z direction and along and normal 
to the streamlines of the basic motion. These will be 
designated as vorticities ¢, n, and ¢ in the directions of 
the streamline, normal, and z-axis, respectively, with 
corresponding unit vectors s, n, and z. 


(1.4 


L'(Og Os) + nU(Og/On) + Fq(0U/dz)] + 


Ug[é(Os/Os) + n(Os/On)] (1.5) 


Now 0s/0s is the rate of change of the direction of the 


(1.6) 


where @ is the angle between the direction of V and the direction of the basic motion far upstream, i.e., the x-axis, 


Fig. 1. 


To obtain 0s /0n consider the continuity equation div V = 0. 


This may be written 


div sg = 0 
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or 
S:(0/Os)Sq + n-(OSg/On) = O whe 
or 
Thi: 
$-S(0qg/Os) + S-q(OS/Os) + n-s(0g/On) + n-qg(Os/On) = 0 1.7) 
By Eq. (1.6) the second term is zero. The third term is also zero, being the dot product of two mutually per- 
pendicular vectors, and hence Not 
Os/On = —(Nn/q) (Oq/Os) (1.8 
Substituting expressions (1.6) and (1.8) in (1.5) 
or 
(o-V)UV = s[EU(0g/0s) + nl (Oqg/On) + fq(OU’/dz)] + n{ L’&q(00/Os) — Un(dg/ds) | 1.9 
The vector 
: i ss This 
(UV-Vjo = Ug(0/ds) (sé + nn + 26) WI 
ice : Sat ae = 1€ 
= sl/q(O&/Os) + nl’q(On/Os) + 2Uq(0¢/ds) + Ug E(Os/Os) + Ug n(ONn/ds) 1.10 
since 02/Os = 0. 
The rate of change of the direction of the normal to the streamline along the streamline is given by sad 
on/Os = —s(00/0s) 1.1] 
Substituting expressions (1.6) and (1.11) in (1.10) 
(UV-V)@ = sUq[(OE/Os) — n(06/0s)] + n Ug[(On/ds) + £(08/Os)] + z Uq(O¢/ds) 1.12)] Hence 
The expressions (1.4), (1.9), and (1.12) may be inserted in Eq. (1.3). First for the z direction 
Uq(o¢/os) = 0 ow 
whence 
r=0 (1.13) 
and 1 
since ¢ = Ofarupstream. In the » direction 
The | 
et 
| | pressure 17 Th 
gradient a 
\}at LE | wens 
| locity 
) J » 
| total circulation pressi 
| (sr static 
a plane 
Three 
Two ¢ 
oe ™“e c. q ol j 
he to th 
“ey Sie 
due t 
pertu 
basic 
the sz 
in the 
Sstagn: 
the pe 
Thi 
: ; limite 
Fic.1. Strut in nonuniform flow: A = unperturbed streamline and B = streamline deflected by secondary effects. 








1.10 


(1.12) 
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— U(0g*?/ds) (OL’/Oz) + UlEg(08/Os) — n(Og/Os)] — Ug{(On/ds) + £(086/ds)] = 0 


whence 


(0/Os) (gn) = —(Og?/Os) (OL’/Oz) (1.14) 
This may be integrated directly along the streamline from far upstream where ¢ = 1 and n = 0 to give 
gn = (1 — gq’) (OU/Oz (1.15) 
Note that 7 will be of the same order as OL’/Oz, unless g—~ 0. In the s direction, noting that ¢ = 0 
U'(Oq?/On) (OU /dz) + EL'(Og/Os) + nU(Og/On) — Lg [0E/Os) — n (00/Os)] = 0 
or 
2q(Oq On) (OL /Oz) — g?(O/Os) (E/g) + nlOg/On + g(08/ds)] = 0 (1.16) 


This may be further simplified by noting that, since the basic two-dimensional flow is irrotational, curl V = 0. 


Whence 


s X (0/Os)sq + n X (OSq/On) 


= 0 


and omitting the terms that are cross products of vectors of similar direction 


Swe 


Ss X g(Os/Os) +nXS (1.17) 
On 
Hence from Eq. (1.6) 
q(06/Os) (Og/On (1.18 
Substituting in Eq. (1.16) for 0g ‘On from Eq. (1.18) and for n from Eq. (1.15) 
—q*(0/Os) (E/g) + 2(0U/0z) (06/O0s) = O (1.19 
and integrating along a streamline from far upstream where £ = 0, g = 1, and @ = O7 
t/q = 2(0U/0z) / (d6/q") (1.20 


The positive direction of & is in the sense of a right-handed screw advancing downstream. 


+ This result can also be obtained by starting from Squire and Winter's Eq. (13 


The analysis, therefore, describes a flow which is 
basically in sheets parallel to the xy plane. The ve- 
locity in each sheet is given by LV plus a small velocity 
v. Neglecting v in comparison with lg the total 
pressure in each sheet is p + | .pl*, where p is the 
static pressure far upstream. Hence, the sheets form 
planes or Bernoulli surfaces of constant total pressure. 
Three vorticity vectors appear, all lying in the sheet. 
Two of them are normal to the basic velocity, the first is 
g(0l’/Oz) (i.e., the grad LU’ X V term in Eq. (1.1)) due 
to the original vorticity and the second is n, Eq. (1.15), 
due to the small velocity v. The third, also due to the 
perturbation, is £, Eq. (1.20), in the direction of the 
flow. All these vorticity components are of 
the same order, provided that l’g does not tend to zero 


basic 
inthe flow. This latter condition implies an absence of 
stagnation points in the basic two-dimensional flow, if 
the perturbation theory is to remain adequate. 

This restriction makes the perturbation theory of 
limited value, since most shapes used in aerodynamics 


and hydraulics have stagnation points at their leading 
edges. The theory, however, becomes useful when 
considering means of minimizing the three-dimensional 
disturbances caused in the flow by the presence of a 
body. The results of the perturbation theory are ap- 
plied in the following section to the flow around a few 
shapes in order to calculate values of £, even though 
its use in some instances is obviously extended beyond 
the limits of the above analysis. 


2) The Secondary Circulation Far Downstream 

The values of the vorticity components and £ can 
be obtained for various shapes by using the solutions 
for two-dimensional potential flow around them. For 
instance, the values of £ obtained by integrating Eq. 
(1.20) along streamlines in the flow around a circular 
cylinder are shown in Fig. 2. At infinity downstream 
£, unlike , does not return to its original upstream 
value of zero. The effect of the disturbing body is to 
produce a vorticity in the direction of flow, which per- 
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sists downstream of the body. The velocities normal 
to the flow thereby obtained are similar to those in- 
duced by the trailing components of vorticity shed 4 


4508 





off wings of finite span. They appear, however, in the 





absence of any lift on the airfoil. 
The total circulation produced in the region between —¢ 


the y- and z-axes (Fig. 1) far downstream is a measure of 2 

















the disturbance caused by the strut. It is given by 2 
5 ; je 882 
Piotal = / J £.. dy dz (2.1) 
0 0 ‘ ~0 964 
where £., denotes the value of at x > ©. Jette 
The circulation in a segment of unit height in the z 7 
direction is given by : a 4 t z + + 
ra 
e «x a 
TT = ¢ , 29 ‘ ES tee r ra 
== f, E.. dy \<-<) Fic. 2. The variation of secondary vorticity along streamlines 
in potential flow around a cylinder. 
Now from Eqs. (1.18) and (1.20) — = component of vorticity in direction of flow 
(2 = initial vorticity, normal to flow 
, e+ co F : Y = stream function 
t.. = 2(dl dz) f (1/g*) (Og/On) Os (2.3) U = initial velocity 
; a = radius of cylinder 
, : : x = coordinate in direction of initial flow 
where the integral is taken along a streamline. Values for x = + are shown at right 
Introducing potential and stream functions for the basic two-dimensional irrotational motion 
Op = GOs Oy = gOn(= Oyvatx = +-—) 
’ a a . 7 
6r = 2(0U/0z) f f (1/qg*) (Og/Oy)O0d OW = —(OU, oz) | [0(1/g2) Y= 5 O@ 
F = 
= —@U oz) ff [1 — (1/q?)] do (2.4) 


where q is the value of g along streamline y = 0. 

The time taken for a particle to pass a distance ds along its streamline = ds/q¢ = d¢/g?. Hence the integral in 
Eq. (2.4) is the difference in time taken by two particles, one on the streamline y = ~ and the other on y = 0, to 
pass from a certain value of ¢ far upstream to another value of ¢@ far downstream. The value for 6I will be large 
if this difference is large. This may be visualized by considering a vortex filament or tube initially upstream and 
at right angles to the flow. As it travels downstream, the end near the body may get left behind and the vortex 
tube stretch until it is nearly parallel to the flow. As it stretches, its cross-sectional area decreases and by Helm- 
holtz’s theorem the vorticity increases. 

Applying Eq. (2.4) to the flow about an airfoil with lift 


wre = 22 tayguis — aren? 
a; = ‘G a) ms u) p 
a i\l/g 1/q ¥ 
= 2.5) 
- ou [* i al ' 
or, = + > [(1/qor7) — (1/q:*)] O¢ 


where subscripts u and / alone refer to flow over the upper and lower surfaces, respectively, and the subscript © de- 
notes the flow along a streamline an infinite distance from the airfoil. The net circulation 


Tet = — = if (1/Geu") — (1/guo1?) dg + J (dg/qr) — / (d¢ at (2.6) 


At a given plane normal to the flow far downstream the value of ¢ for the upper surface region is greater than that 
for the lower region by A, the circulation around the airfoil. Also, due to the circulation around the airfoil, 
Jou = 1 + e and g.; = 1 — e€, where ¢ represents the velocity due to the circulation which tends to zero at 


infinity. Hence the first integral in Eq, (2.6) is not zero. 
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On integrating between two planes normal to the flow, one an infinite distance upstream and the other an in- 


finite distance downstream 


r ° 


Gon) — (ie 


f- (d¢g qi) ~f- dg "}) = | 


The two integrals cancel out except on the surface of the 
airfoil, hence 
K} (2.8) 


Ol’ Oz } £ (ds q) 


where the integral is taken around the surface of the 


Oo. = — 


airfoil, passing from the leading edge to the trailing edge 
along the lower surface and back along the upper sur- 
jace. The positive direction of 6I is in the sense of rota- 
tion of a right-handed screw advancing in the down- 


stream direction. 


3) The Effect of Leading Edge Shape 


Since at a stagnation point Lg — O and & tends to 
infinity, only profiles with cusped leading edges may 
be analyzed by the perturbation theory. Two different 
types of profile with stagnation points at the leading 
edge may, however, be distinguished. In one the value 
of 6f calculated from Eq. (2.4) becomes infinite and in 
the other 6% remains finite. It may be expected that 
the latter profile offers less disturbance to the flow than 
For profiles with a leading edge angle less 
When the leading edge angle 


the former. 

than z, 6I’ remains finite. 

is equal to or greater than 2, 6I becomes infinite. 
This may be shown by examining the value of {“d¢ + 


w° near the leading edge. The procedure is to trans- 














Y 
y 
Yq 
Xx 
0 re 
- { ‘Weo 
3 v 
W = CONSTANT 
y 
nl 
pro 
Fic. 3. Mapping of flow around circle in inverse hodograph 
plane 
Z= X +1) 
= velocity at any point 
6 = angle between velocity vector and direction of initial 
velocity 
¥ = stream function 


Initial velocity is taken as unity. 


dg = f- [A/geu) — (1/gei)]ds-—K 


t 
a 
| 
- 
A, 
| 
>>) 
| 
| 
to 
a 
7 
~m 
to 
~1 


form the flow at the leading edge of a circle in the ¢ plane 


to the flow around a leading edge of angle \ in the z 


plane. If q’ represents the velocity on the stagnation 
streamline in the ¢ plane 
S (do go-) = -_ (dz/dt) *(d@ q 3.1 
The flow in the ¢ plane about a circle of radius / is given 
by the coniplex potential 
W = [t+ (7/t)] 3.2 


and the Karman-Trefftz transformation will be used, 


ViZ., 

(s+ anl)/(s — nl) = (t+ )/(t—D]" (3.3) 
where n = 2 — (\/m) and X tis the angle of the leading 
edge. For flow near the leading edge where ¢ = —/ 
ands = —nl 

t= —/ + re’ 
. i 34 
s= —nl + se 
where ry and s are small. Substituting in Eq. (3.3) 
gives 
= a 97)! = 
S€ = ey el NM \ Zi é oO. 
Whence 
(dz/dt)| = (ds/dr) = n*(r/21)"~° 3.0 
both when a = az, 8 = wand when a = 7/2, B = 


\/2. Near the stagnation point substitution for ¢ in 
Eq. (3.2) gives 
@ = — 2] —/(r/l1)*? cos 2a 3.4 
whence 
do@/dr = —2(r/l) when a = r 
and 
= 2(r/l) when a = 7/2 3.8 
Hence the velocity gy’ is given by 
Q’ = 2Ar/l 3.9 


Substituting the expressions (3.6), (3.8), and (3.9) in 


Eq. (3.1) and integrating from a distance ¢ ahead of 


the stagnation point at a = 7 toa distance ¢ along the 
circle from the stagnation point at a = 7/2 
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Y whence from Eq. (4.1) 
dz = Z(dW/dZ)dZ 
yee 
° = Z((1/2V R) (Z -—1)°"% - 
Yo (V R/2) (Z — 1)~*] dz 
IRCLE OF 
RaDIUS= R 
Yq . and 
C] 
a 
=. = , ~ ¢= 2 + 2)/31+/iz — 12/2 + 
oo ; ; 
(2—2Z2)VR/(Z—1) (44 
& On the circle 
\ iy ta ~- 
Z=1+ Re (4.5 
, lca whence the circle maps into a profile whose coordinates 
Fic. 4. Circle in inverse hodograph plane : ; 
in the z plane are given by 
- . Aer x = (4/3)R cos® (a/2) + 2(1 — R) cos (a/2 
"do ln4 a 3 “ale “al” - 
~e q? ™ g2n—l A ] ( 1. y = — (4/3)R sin* (a/2) (4.6 
e [y\2"-3 Ly In4 © /y\2"-3 dp The ratio of the dimensions R ‘I in the Z plane deter- 
j d ] = 520-2 ] ] (3.10) mines the thickness-chord ratio, ¢/c, in the z plane. 
0 - 0 . . - ~ 2 
By inspection of Eq. (4.6 
> inteor. comes infinite when 2n — 32 — 
and the integral becomes infinite when 2” 3. Ll. t/c = 2R/(3 — R) 


i.e., When A > 7. 


Part I]—APPLICATION TO SOME SPECIFIC PROFILES 


(4) A Bicusped Profile 


The form of Eq. (1.20) suggests the mapping of the 
basic two-dimensional motion in the inverse hodograph 


plane. The transformation is 


Z = (dz/dW) = (1/qje~” (4.1 
where II’ is the complex potential, z the physical plane, 
and Z = X + 1Y, the coordinates in the inverse hodo- 
graph plane. 

Fig. 3 shows the streamlines of the flow about the 
bottom half of a circle mapped in the Z plane. The 
streamlines cover the entire plane to the right of 
X = '/, and form closed curves beginning and ending 
at X =1, Y = 


streamline is given by 


A = (1/2) [, 


A profile for which the secondary effects are to remain 
small must have a small enclosed area when mapped in 
the Z plane, unlike the circle for which the streamline 
yY = O encloses an infinite area, Fig. 3. 

Suitable profiles are, therefore, obtained by mapping 
curves enclosing small areas in the Z plane into the z 
plane. The simplest of such profiles is obtained from a 
circle of radius R in the Z plane with its center at 
(+1, 0), Fig. 4. 

The flow, which is inside the circle, is given by 


0, and the area, A, enclosed by a 


(d6/q*) = &. /[4(dU/dz)] (4.2) 


W=V(Z—1)/R+ VR/(Z—-1) (4.3) 


whence 
R = (3t/c)/[2 + (t/c)] (4.7) 
The chord, c, in the z plane is given by 
c= 41 — R/S 
whence 


c = 4/[1 + (1/2) (¢/c)] (4.8 





0-4 FIG. 5(a) 


















FIG. 5(b) 
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Fic. 5. Bicusped profiles; (a) shapes and (b) surface pressure 
coefficients. 

t = thickness 

¢ = chord 

R = radius of circle in inverse hodograph plane 
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and the thickness x/e = (1/2) cos (a/2) [1 — (t/c) sin? (a, 2)] 
t = 8R/3 = (4t/c)/[1 + (1/2) (t/c)] 1.9 ye = —(1/2) (t/c) sin? (a/2 
The coordinates of the profile may be written as The velocity on the surface, go, is given by 
l 7 1 + (t/c) + (5/2) (to)? + 3(t/c) [1 + 2) (t/c)] cos a 
q - [1 + (1/2) (t/c)] 
Profiles and the surface pressure coefficients are shown in Fig. 5. 
5) Bicusped Profile—Downstream Vorticity 
Writing Z = 1 + re and substituting in Eq. (4.3), the stream function of the flow inside the circle is 
¥ = (Vr/R — VR/Pr) sin a«/2 
whence on a streamline 
r/R)? — 2(r/R) [1 + (1/2)? cosec? a/2] + 1 = 0 


Substituting 


and 


r/R = 1+ 2[(1 — k?)/k?] sec? d — IV 1 — k2/k?2] sec? 6V 1 — R? sin? @ 


This relation between 7 and ¢ on a streamline may be used to find the area, A, enclosed by a streamline, for 


. . 


A= oR: f (r/R)? d(a/2) = -2R: | (r/R)? de 
0 0 


The solution is obtained by substituting the elliptic functions 


























cos @ = cnu; V1 — kk? sin? d@ = dnu 
and 
V1 — &? = R’ 
——— - - 
4m l I | ae | 
7 EXACT SOLUTION 
——— EXPONENTIAL APPROXIMATION | 
\ 
4 T 0.0 
10 $+—+- | + + + ; + os +—_+-—_+—_ 
\ 
\ 
8 a ee ee | J - + + Ee 
4$e0(1+3 % LW | | | | anaes 
tz 2 aU N 
9(%) 63 ‘2 See e 2a 0.02 
\ 
ae «Be ae a f(n 
‘ (n) 
e118. —_ 
N 3 
pe + + + +—_—4 9 
i oo! 3 
——1 pt = eS 
i — 
+ - : . + « 
end 4 ee eae Coe ee | 
=> 2 a 
° 4 4 hat = 
°o ' \ 
w= 4(%)/(i+¢ %) 
. : , ° ss ° — 
Fic. 6. Bicusped profile: secondary vorticity far downstream. 0 :; 
Exponential approximation: & 
Em [1 + (1/2) (t/c)]? _ Or exp \ 10x(y/C) Fic. 7. Bicusped profile in exponential boundary layer 


(t/c)? (dU/dz) 





oe ‘ Vari- 
4 4[1 + (1/2) (t/c)]\ ation of f(m) |Eq. (6.15)] with 6*/c for t/c = 0.05 and 0.25, 





596 JOURNAL OF THE AERONAUTICAL SCIENCES SEPTEMBER, 1954 
where 
“6 do 
“= 
o V1—k* sin’ 
Then 
r/R = 1 + (2k'?/k*cn*u) — (2k'/k*) (dnu/cn*u) 5.5 
and the area becomes 
A= one | ,dnu — (4k'/k?) (dn*u/cn*u) — (Sk’ */k*) (dn?u/entu) + 
0 
(Sk’ 4/k*) (dnu/cntu) + (Sk’ 2/k®) (dnu/en2u){ du (5.6 
and on integration 
A = 2R?} (9/2) + (4/3) (R’/R*) (4 + RE — (4/3) (R'/R*) (4 — RK! 57 
where A and E are the complete elliptic integrals 
. ete do : wis : 
i = / —= , E= / V1 — k* sin? ddd 
0 V1— k* sin’ @ SV 
Hence from Eq. (4.2) 
£ = {(t/c)?/[1 + (1/2) (t/c) 2} (dU/dz) {9a — 12y V1 + (2/4 [(¥? + 3A — (¥? + 5)E]} .S 


where y = y or substituting for the chord from Eq. (4.8) 


The net circulation 


or = J &. dy = (1/4)c{l + (1/2) ¢¢ al f £. dy (5.9 


and by direct integration of Eq. (5.8) this gives 


a = c-(dU/dz)2-4} (t c)?/{1 + (1/2) (t/c)]j (5.10 


Fig. 6 shows graphically the relation for £,, derived above and an approximate exponential relation which has the 


same value of 6f and of £, aty = 0. This approximate relation is 


&. [1 + (1/2) (¢/o) 
(t/c)? (dU /dz) 


j Lor(v/c) ) 
= 9rexp )— - = (5.11 
(  4[1 + (1/2) (t/o]f 


Fig. 9 shows a curve of 61'/(c-d/dUZ)], vs. thickness-chord ratio ¢/c. 





(6) Bicusped Profile—the Energy in the Downstream 

Secondary Flow 

Far downstream of the strut the secondary vor- 
ticity, £, induces secondary velocities v and w in the y 
and z directions, respectively, in a manner similar to 
those induced by the trailing vortices shed off a wing of 
finite span. The energy in this flow is of interest since 
in wing theory it is related to the induced drag of the 
wing. 

The equation of continuity for the flow far down- 


stream is 
[o(U + u)/dx] + (Ov/Oy) + (Ow/dz) =O (6.1) 


where u is the small perturbation of the velocity U’ in the 
x direction. 

The pitot pressure far downstream is p + 
(1/2)p{(U + u)? + v? + w*} and this must equal the 


corresponding pitot pressure far upstream, viz., p + 
(1/2) pU?. 

As the pressure p is the same this relation shows 
that the velocity, u, is of lower order than v or w. 
Now 


0U/dx = —(w/U) (0U/dz) (6.2 


and when w and OL’/dz are both small, this term also 
is of lower order than Ov/Oy or Ow/dz. Hence the 
first term in Eq. (6.1) may be neglected and the equa 
tion satisfied by introducing a stream function x, so that 
w = Ox/Oy and v = —(0x/0z). The flow is then 
given by 


V*x = &, (6.3) 


Now the total energy is given by 
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U = U,[1 — exp (—2/6 


*)] 6.6 


where 6* is the displacement thickness of the boundary 
) is the 


D.* = 1/2 p / / \v° + ww dy dz 
. VU . 
1/2 | J (Dy /dy)? >y/d2)2]l dy d layer which this expression represents, and [ 
(l/2)p y)° + Z)-;aydaz 7 . ‘ : : a 
"i MI. . si value of the approaching velocity U asz—~ oo. Then 
6.4) dU/dz U,/6* exp (—2 6* (6.7 


and from Eq. (5.11 


Integrating by parts with x 0 on the boundaries and 
Onr(t/c)?U, 
“ a it \ je OS 


at infinity 


DD? == » | / x | 0*x/Oy") + (0°x/0z ] dy dz ; 
e~* "exp(— — 
-" ( 1[1 + (1/2 


y’ and 2/né* 2’, where 


tied / / E,xdydz (6.5) 
ve Writing y/né6* = - 
where 6 is the thickness in the z direction of the non- n = 4f1 + (1/2) (t/c)]/(1546*/c), 
This in- 
9r(t/c)7l 5 z 
€ é é 6.9 
1 + (1/2)(t/c) }*6* 


uniform portion of the approaching stream. 
tegral has been solved numerically for the experimental 
A particular solution of Eq. (6.3) is then 


conditions described in Section 10.3. 
To investigate the effect of 6 and the thickness-chord 
ratio an analytical solution has been obtained using 
the exponential approximation to ~, [eq. (5.11)] and sa 162° Uo I en-' e 6.10) 
n approaching velocity of the form "25m 8* (1+ n? 
To satisfy the boundary conditions x = Oat y = Oand s = O the following solution of V*x = 0 is introduced 
16¢2 Us l 2 { : 5 : , : 
X = =— . e sin As’ dX e “sin Au du + 
257 6* 1t+n? x tJo JO 
/ e~* sin Ay’ dv / e “sin Au dp (6.11 
J 0 / 0 
rhe energy in the secondary flow from Eq. (6.5) is therefore 
10 _ 9 ~ ~ 
- es 144(t/c)! 4 9" p— 202" dwt ge! 
yy =—peUge — gas ‘ ee ay @=— 
the © 25[1 + (1/2)(t/c) 21 + n? Jo So 
> : E , , ; 
J / / e~ ANY e-™ sin \s’ — - dd dy’ dz’ 
7 T Jo Jo Jo A* -P 8" 
2 F . A+ 2 4 A , a > 
| e er dy dy’ dz { (6.12) 
"JO FO fe rl 
The integral terms in the braces become 
l 2 : d” 2 ; Ye 
a : . dx — / ; - dX = 
ws 47 T Jo (A +1) (A° +n tT Jo A +n) (A? + 1)° 
ww. , 
l 21 «x {n* — 1\? l1— n’ l 
- Sarah + — log. n + (6.13) 
$72 am Lain \n- + 1 (l +n l+n 
9 
Hence the total energy may be written 
Iso 
yLy7c? 144(t/c)4f (2) 
he a 6.14 
25[1 + (1/2) (t/c)]? 


la- 
. | (6.15 


iat 
where 
en 
n° 24a /n? — 1\? 1 — n° 
f(n) = z : -+- — log. n + > a 
2) 1+ n?La l4n \n? + 1 (i = 2°)* l + n* in | 
Che value of f(m) is given in Fig. 7 as two curves plotted against (6*/c) for t/¢ = 0.05 and 0.25. 
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2.0 The potential and stream functions are given in terms of 
the coaxial coordinates by the 
ie C sinh 2V/n alc 
= 2 ; to 
iT n cosh 2Y/n — cos 2X /n . 
iro 
16H ° , 
( sin 2X /n x-d 
2 L= : : (7.3) 
n cosh 2V/n — cos 2X/n 
! 4 rt . . . 
Ihe velocity is given by 
‘ ant 
2 16 ( cosh } — cos X ) 
! } , = _ 
3 oa = (4.4 
é n*\cosh 2Y/n — cos 2X/n 
=> 
Ae) 
1.0 and the angle @ by 
(1/2)@ = tan (cot X/2-tanh Y/2) — 
0.8 tan! (cot X/n-tanh Y/n) (7.5) ar 
CYLINDER | a eee eee ie) have thee termine 
: mamee $s Some values tor &./(dU/dz) have been determined 
oe SX! BICONVEX 5.5 graphically by plotting 1/q* against @ for selected values 
4 BicusP 5.5 th ae aia oni ie en 
. BICONVEX 7.596 of the stream tunction, y, and measuring the area en- 
SICUSP = 7-596 closed. This has been done for two thickness-chord Th 
pak ratios and the results are shown in Fig. 8, for compari- 
§ son with the bicusped profile. Curves for the two pro- 
Po files only differ at values of the coordinate y less than 
the thickness of the profile (¢/c < 0.2). 
A more significant comparison may be obtained from 
"5 02 4 nY3 oY] io 7 ra the value of 6F since this gives some measure of the Sin 
Wt order of the secondary velocities and the energy of the larg 
Fic. 8. Secondary vorticity, &., for various profiles sece mdary flow. exp 
Q = initial vorticity To determine 60 analytically, Eq. (2.4) is used: pos 
y = coordinate normal to original flow z in F 
t = thickness of profile re dl t | : 
5] = <- ] — dd = I 
dz. qo’ 
T ‘ a . : : au | = 
Ihe energy thus varies as the fourth power of the =. 1——)d@ (7.6) 
thickness and, from Fig. 7, reaches a maximum when dz Jo qo" 
6*/c is approximately 0.1. This maximum is explained — gincee the profile is symmetrical about the y-axis. indi 
physically by noting that when the boundary layer is 
thin the vorticities are large but the region of flow af- V 
fected is small. When the boundary layer is thick, the 0-3 a sho’ 
region of flow affected is large but the vorticities are of 6 
small. valt 
: BICONVEX 8) 
7) The Biconvex Profile 
A profile made up of two circular ares shaped like 0-2 fF T 
a biconvex lens has the advantage that it is easier to elliy 
: ; : . ¢ 
manufacture than a bicusped profile. With a sharp ér BICUSPED So! 
leading edge it gives a finite value for 6% and hence cH P ol ¢ 
finite secondary velocities, even though £ tends to abo 
infinity on the stagnation streamline. 0"! file 
The plane, potential flow over such a section is given aie 
by Milne-Thomson® rat 
valu 
a vA ic | 
, r - To 
W = cot —, z= 5 cot 2 (7.1) “8 
n n 2 ERE a ee a adec 
0 0-1 re 0-3 ; 
where c is the chord, z = x + 7ty coordinates in the CG 9) 7 
physical plane and Z.= X + 7Y coaxial coordinates. Fic. 9. Values of 6&/c(dU/dz) for bicusped and biconvex pro- T 
rr . . ° iles ni Ss s , +4 7 
The thickness, f, is given by files. Points O show values of 
sr 8 e423 ' { ( t\2 resu 
t = c cot n(7/4) (7.2) cdU/ds) 3 ¢ ) " 3 Xe ) lar a 
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Now the flow along the stagnation streamline where On the profile Y = n(a/2) and 
the velocity is go may be divided into two parts, that 
along the profile from the mid-chord position, « = 0, @ = c/n tanh V/n 
to the leading edge, x = c/2, and that along the x-axis _ 
from the leading edge to infinity. For the part on the 
x-axis X = 0 and from Eqs. (7.3) and (7.4 aha 16 [cosh Y — cos n(x/2)] - 
Pere ey n} t cosh? V/n 
~~ Along the profile from « = 0 to c/2, Y goes from 0 to 
+o, Along the x-axis from the leading edge 
alien) 16 sinh? Y/2 (7.7) © 2 to +, VY goes from + to 0. Hence Eq. (7.6 
n*sinh* Y/n may be written, using Eqs. (7.7) and (7.5), 
duc *0 n+ sinh? \V/n 
sr —2 i-—- — -— : )a (coth Y/n) + 
dz n | ( 16 sinh* Y/2 


, n' t cosh? V/n : ae 
l— - d (tanh V/n | 1.9 
0 16 [cosh Y — cos n(2/2) |? 


The integrals are determined by contour integration in the Appendix, with the following result: 


6r "\* nT 2a 25 nT nT a (2 n 2r 1 /n\? 
: = cosec- l— cot - cot - _ cot — (.10 
c(dU /dz) 2 2 n n 2 ae 3 \n 2 n 3\2 


Since 1 < t/c < O, therefore 1 <u < 2andnis1.6or in Fig. 10. The values of £, for the flow over the upper 
larger for t/c of 0.3 and less, this expression may be surface are of opposite sign to those for the flow over 
expanded in terms of powers of [1 — (m/2)] for pur- the lower, i.e., concave, surface. The difference is due 
poses of calculation (see Appendix). Values are shown _ to the larger value of g on the upper surface compared 
in Fig. 9. with that on the lower surface. Downstream of the 
Points showing values of the formula trailing edge, spanwise velocities of the same signs but 
a differing magnitudes are induced by the secondary vor 
a Pe. (t/d211+ . (t/c)2 (7.11 ticities. Hence a vortex sheet must pass from the trail 
c(dU/dz) 3 3 ing edge downstream in a similar manner to the trailing 
— Tener oe ‘ , vorticity shed from a wing of finite span. The strength 
indicate that this is a satisfactory approximation. age ie : 
of this vortex sheet will be determined by the magni 
Values for the bicusped profile, Eq. (5.10), are also tude and distribution of the secondary vorticity & 
shown in Fig. 9. The biconvex profile has a value In the example shown in Fig. 10 the sign of the vorticity 
of 6’ somewhat but not very much greater than the in this vortex sheet is in the same sense as that shed 
value for the bicusped profile. . off a finite wing whose circulation is decreasing in 
the direction of L’ (the approach velocity) decreasing. 


8) The Elliptic Profile 


-The results of calculations of €. for a circular and an 2-05 

elliptic profiie are also shown in Fig. S. The values of CAMBER = 10% 

§./(dU/dz) were determined by graphical integration — 

of curves of 1/g? vs. 6 for various streamlines in the flow j 

about the profile. The values for £_ for the elliptic pro- hes | snide eiieuniimnesin 

file are appreciably greater than the values for bicon- Taal 4 — 

vex or bicusped profiles of the same thickness-chord f 

ratio. The rounded leading edge implies an infinite 0-4} ~ : 

value for 6I, and for this profile as for others with ae SURFACE i eee 

rounded leading edges, the perturbation theory is in- se aa 08 3 ; 

adequate. 2 cosp = 

9) Thin Circular Arc Airfoil at Zero Incidence wes. 10 — betty socte 2 ‘om ben a els 
To show the effects on a profile designed to give lift, 8 = 11.31° camber = tan 8 = 10% 

results of calculations of &,, by Wilson® for a thin circu- an — normal to flow 


lar are airfoil of 10°7 camber at zero incidence are shown Q = initial vorticit 








600 JOURNAL OF THE 


AERONAUTICAL 


SCIENCES SEPTEMBER, 1954 








Fic. 11. 
of sand in cms. 


The net secondary circulation, 1.e., the sum of the 
circulation due to the downstream secondary vor- 
ticity in the flow over both upper and lower surfaces, 
may be obtained from Eq. (2.8), i.e., 


i sew dU j ds Les K| 
dz | q f 


The Joukowski transformation may be used to evaluate 
s and q on the airfoil surface, giving the result 


(2.8) 


a dl 
Of nt = C 27 tan B[1 — 1/(cos! B) | (9.1) 
dz 
where tan § is twice the camber and cis the chord. This 


expression is negative, agreeing with inspection of Fig. 
10, from which it is obvious that the vorticity in the 
flow over the lower (concave) surface dominates that 
in the flow over the upper surface. The sign of the vor- 
ticity in the vortex sheet trailing from the airfoil is 


positive. 


For airfoils with rounded leading edges, the vorticities 
and secondary circulations in the flow above and below 
the airfoil tend to infinity and the perturbation theory 
fails. However, the will still remain 
finite and may indicate whether the secondary vor- 


value of 62 ye 
ticities (and the corresponding trailing vortex sheet) 
will induce large or small asymmetric velocities outside 
the boundary layer. When wings are immersed in a 
boundary layer, as, for instance, at the junction of wing 


Scouring contours for elliptical and bicusped profiles at zero incidence. 
g I is E 


Figures on curves give depth 


Depth shown negative, height positive, and thickness of profiles is 1 in. 


and fuselage on an aircraft, it is evident that the dis 
turbance or interference will be affected by the shapé 
of the wing profile, i.e., differing profiles designed to 
give the same lift will not necessarily produce similar 
disturbances. 


Part III—EXPERIMENTAL RESULTS 


Two experiments were designed in which the bicusped 
profile was compared with other profiles. In the first 
simple experiment the phenomenon of scouring in a 
water channel was used to indicate that the use of this 
profile lessens the disturbances caused in the flow of a 
nonuniform stream about a strut. In the second, the 
drag of a bicusped profile in a nonuniform flow was 


compared with that of a conventional strut. 


(10.1) Experiment in a Water Channel 


A strut of an elliptical profile and a strut of a bicusped 
profile of the same thickness (1 in.) and chord (5.5 in. 
were placed side by side in an open channel along 
which water was flowing at an average velocity of be- 
tween '/, and 1 ft. per sec. The channel length was 
large enough to ensure that there was a velocity gra- 
dient, the velocity being lowest at the bottom. The 
bottom of the channel was covered with sand to the 
depth of 2 in. and the depth of water above the sand 
was 3in. The water was allowed to run for a period of 
half an hour before being drained off, in order to permit 
observable changes in the level of the sand about the 


struts to develop. The results of a typical set of meas- 
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urements? of the level of the sand are shown in Fig. 11 
in the form of a contour map. The sand about the 
bicusped strut is markedly less disturbed than that 
about the elliptical strut. The experiments were re- 
peated with the struts yawed, and showed a less niarked 
but nevertheless definite reduction in the disturbances 
about the bicusped profile, Fig. 12. 

The results of this simple experiment suggest that the 
scouring in rivers and water channels is at least par 
tially a secondary flow phenomenon and that the per- 
turbation theory offers a method of reducing such 
effects. 

10.2) Measurements of Drag 

Two struts 2.4 in. thick, 12-in. chord and 42-in. span 
were used by Ling’ for experiments in the student wind 
tunnel at the Massachusetts Institute of Technology. 
The profile of one strut was based on NACA 16-006 and 
the other on Eqs. (4.7) for the bicusped profile. They 
are compared in Fig. 13(a). The drags were measured 
first with a uniform flow approaching the strut and then 
with a flow with a spanwise velocity gradient. The 
total drag on either profile was expected to be about 
1 lb., but the difference between the drag in the uniform 
flow and that in the nonuniform flow was expected 
to be of the order of | per cent of the total drag. A four- 
bar balance system was designed to keep the strut at 
zero angle of attack and to eliminate all gravity forces. 
The drag was measured by summing the readings of 
eight strain gages, four on each set of four bars. To 
eliminate effects due to the unsteadiness of the flow in 
After 

was 


the tunnel, oil-immersed dampers were used. 
careful calibration an accuracy of +0.0025 Ib. 
thought to be achievable. 

The nonuniform velocity was produced by a series of 

s-in. thick pieces of plywood spaced !/2 in. apart with 
their trailing edges 12 in. in front of the strut profiles, 
Fig. 13(b). With the struts removed from the tunnel, 
the velocity distribution in the wake produced by the 
plywood assemblies was measured by a _ pitot-static 
rake placed 12 in. downstream, i.e., at the position nor- 
mally occupied by the leading edge of the strut. A 
typical traverse in the wake of an assembly is shown in 
Fig. 13(c). Different plywood assemblies were used 
to produce wakes of different displacement thicknesses. 

To ensure that turbulent separation was achieved at 
all conditions of test, a trip wire of 0.0045 in. diameter 
was placed 8 in. from the leading edge of the bicusped 
profile and a wire of 0.007 in. diameter 4 in. from the 
leading edge of the NACA profile. With these wires 
in position there was only a very small variation in the 
normal drag coefficient from Re = 5 X 10° to Re = 9 X 
10°. 

In addition to the drag measurements, explorations 
of the wake from the strut were made with a five-hole 


tube. From these it was possible to calculate the 


+ Made by A. T. Ling, Gas Turbine Laboratory, Massachu- 


setts Institute of Technology. 
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values of the induced velocities downstream of the 
strut. 
If D, is the drag of the strut of chord c in the uniform 


stream then 


dD, = (] 2) pUo?Cps.c. 


where s is the span and L% the approach velocity. 

In the nonuniform stream there is an increase in 
drag D* due to the secondary effects and a decrease in 
drag AD, due to the decrease in the approach velocity 
in the wake of the assembly producing the shear flow 
upstream. This decrease in drag has been estimated 
from the experimental results by assuming that the two- 
dimensional drag coefficient is unaltered by the up- 
stream wake and 


(1 2)oCoe f 
a Span 


AD,/D, = f [1 — (U 


AD, = (U)? — U*) dz 


Uo)?] d(z/s) 


The results of the experiments are summarized in Table 
1. The values of 6* shown in the table are values of 
the displacement thickness of half the wake produced by 
the plywood assemblies. 

The increase of drag due to the secondary effects 
D* is very small, being only a per cent or two of the 
total drag, and frequently smaller in magnitude than 
the decrease in drag AD, due to the reduction in ve- 
locity of the main stream. The bicusped profile has, 
however, a markedly smaller secondary drag D* than 
the NACA profile. This is shown in Fig. 14, in which 
a secondary drag coefficient Cp* for half the wake, 
defined in terms of the thickness of the strut from the 
equation 

(1/2)D* = (1/2)pUo?Cp*?? 
is plotted for the two profiles. The estimated range 
of experimental error is also shown in the figure. 


10.3) Energy in the Secondary Flow 


Measurements of the velocities in the wake of the 
struts at a station 0.1 chord downstream, were made 
with a wind speed of 72.5 m.p.h. and 6*/c = 0.0496. 
From these measurements, values of the energy in the 
secondary flow normal to the main flow were obtained 
by graphical integration over the entire yz plane. 
Values of J’ J (1/2)p(v? + w*)dy dz, where 2 
the velocities in the y and z directions, respectively, 
were 0.00037 for the bicusped profile and 0.0107 for the 
NACA profile. The values are probably too small be- 
cause the coverage of the flow field was inadequate. 


and w are 


Simple airfoil theory suggests that, apart from viscous 
integral should equal the induced or 


There is, however, some discrepancy 


effects, this 
secondary drag. 
as the corresponding measured drags were 0.012 and 


0.028 for the bicusped and NACA profiles, respectively. 
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Fic. 15(b). Secondary flows 0.1 chord downstream from trailing 


edge 


of strut profiles. Bicusped profile, 6*/c = 0.05. 


The energy in the secondary flow for the bicusped 
profile was then determined theoretically from Eq, 
(6.5). 

For the above experimental conditions, x for the bi- 
cusped profile has been determinedt by a numerical 
(relaxation) method from the relation for —, given in 
Eq. (5.8) and the upstream velocity distribution used 
in the experiment. The result for D,* was 0.0008, 
which is of the right order of magnitude when com- 
pared with the value of 0.00037 obtained from the wake 
traverses, but is considerably smaller than the value of 
D* = 0.012 deduced from the drag measurements. 

This suggests that some other factor such as a dis- 
turbance of the boundary layer on the strut must be 
considered if the secondary drag is fully to be accounted 
for. Similar difficulties have arisen in cascades and 
bends in attempting to estimate the losses from the 
energy in the secondary flow alone. 

The results of measurements of the induced velocities 
at 0.1c from the trailing edge, with a wind speed of 72.5 
m.p.h. with 6*/c = 0.05, are shown in Fig. 15. The veloc- 
ities measured with the nonuniform flow were corrected 
by deducting velocities measured in uniform flow, and 
the corrected results only are shown in the figure. 
Actual measurements are shown by arrows whose length 
is proportional to, and which are oriented in the direc- 
tion of, the corrected velocity. Other lines are intro 
duced to indicate more clearly the direction of flow. 
Fig. 15(b) shows the existence of four distinct vortices. 
In Fig. 15(a) for the NACA profile there appear to be 
only two distinct vortices of much larger magnitude. 
These two vortices also appeared one-half chord down- 


stream. 


(10.4) Discussion of Experimental Results 


The growth of secondary vorticity in the flow about 
a body results in appreciable disturbances, as shown 
by the scouring of sand in the water channel and the 
flows downstream of the struts. These disturbances 
may be reduced by using shapes which reduce the 
secondary vorticities. 

There is a significant difference between the secondary 
drags of conventional and bicusped shapes. This dif- 
ference would probably be larger for a conventional 
shape with a larger radius of leading edge. The results 
show that the drag of a body immersed in a nonuniform 
stream can be reduced below its drag in a uniform 
stream of velocity equal to the maximum velocity in 
the nonuniform stream. 

Although the secondary drag is small, it is appreciably 
larger than the energy in the secondary flow. This dif- 
ference may be due to a change in the boundary layer 
on the profile, caused by the secondary flow. The 
secondary velocities tend to sweep the boundary layer 
in a spanwise direction along the strut, causing it to 

+ By W. D. Armstrong, Cambridge University Engineering 
Dept. 
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SECONDARY F 


Drag in 
Uniform 


LOW 


ABOUT 


TABLE | 


STRUTS 


AND 


Drag Measurements on Strut Profiles 


Displacement 
Thickness of 


Estimated 

Fractional 

Decrease, 
in Drag Due 


Measured 


Drag in 


Nonuniform 


AIRFOILS 


Drag Due to 


Secondary 


Effects, D* 


Normal Drag 
Coefthicient, 


Stream, half the Wake, 
D,, Lb. 5*/c to Wake AD, /D Stream, D:» D. — (D; — AD,), Lb. Cp = D pl? 
Bicusped 
0.678 () 0172 0.0180 0.673 0.007 0 O183 
0.678 0.0499 0.0518 0.653 0.009 0 O1SS 
0.678 0 0695 0.075 0). 6438 0.012 0 OSS 
0.859 0.0168 0.0180 0.854 0.009 QO O1LSS 
0.859 0.0496 0.0533 0.829 0.012 0 O83 
0.859 4}. 0700 0.0752 0.814 0.015 0.0183 
1. 12 0.0175 0.0188 1.110 0.009 0 0185 
1.120 0.0496 0.0515 1.080 0.015 0 0183 
1.115 0.0715 0.0736 1.055 0.019 0.0185 
NACA 16-006 
0.699 0.0172 0.0180 0.701 0.022 ISO 
0.700 0.0499 0.0518 0.690 0.025 0.0189 
0.700 0.0695 0.075 0.675 0.023 0. O1S9 
0.877 0.0168 0.0180 0.893 0.030 0.0187 
0.893 0.0496 0.0533 0 877 0.028 0.0190 
0.898 0.0700 0.0752 0.867 0.032 0.019] 
1.152 0.0175 0.0188 1.107 0.085 0 0189 
1.157 0.0496 0.0515 1.136 0.037 0.0189 
1.162 0.0715 0.0736 1.115 0.036 0.0191 


As the boundary 


layer is swept away, it must be replaced by air of higher 


energy, with a resulting increase in shear forces and 


drag. 
affected by the secondary velocities. 


The form of the boundary layer may also be 
These possi- 


bilities need more detailed experimental investigation. 
Although the experimental survey was not complete 
and was confined to only one condition of test, the fair 


agreement between the measured and theoretical ener- 


gies in the secondary flow is promising. 


It suggests 


that the theory enables secondary velocities to be pre- 
dicted with satisfactory accuracy. 
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11) Conclusion 


The theoretical methods used in this paper are easily 
applied to determine the relative magnitudes of second- 
ary flows. The assumption that the flow is inviscid 
does not necessarily invalidate the comparisons made 
possible by the theory, for shearing stresses are likely 
to reduce rather than augment secondary velocities. 
The main conclusions of applying the theory to struts 
and airfoils relate to the importance of nose shape, to 
the possibility of designing profiles to minimize sec 
ondary effects, and to the nature of the variation of 
the energy in the secondary flow with strut and bound- 
ary layer thicknesses for the particular case studied. 
A biconvex profile gives a good practical shape for 
minimizing although the 
slightly less practical bicusped shape is, according to 


secondary disturbances, 


theory, even better. 


The calculations for only one lifting airfoil shape have 
not been extended to more than a determination of the 
values of &,. The existence of a vortex sheet trailing 
off the airfoil is indicated. 
along the span outside the approaching boundary layer, 


This vortex sheet extends 


but it is not clear whether its strength is related di- 
rectly to the change of circulation along the span of the 


airfoil. 


The experimental results confirm the major effects 
predicted by the theory, particularly regarding the 
influence of nose shape on the secondary disturbances. 
It seems likely that the changes in level of river beds 
which occur around obstructions, and the depressions 
found in snow at the windward side of a tree or tele- 
graph pole are caused by the scouring effect of the 
secondary flow. The disturbances of the boundary 
layer by small excrescences on the wings of aircraft 
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are similarly due to the strong secondary vorticity gen- REFERENCES 


erated by the presence of a stagnation point. eit Wh: andi ean. WS. Tete Lane Theory 1 


The wind tunnel experiments by Ling confirm the Wing in Non-Uniform Flow, Quart. Appl. Math., Vol. III, p. 1, 


prediction that profiles can be designed to minimize April, 1945 
secondary flows. The drags due to them are small and * Squire, H. B., and Winter, K. G., The Secondary Flow 
difficult to measure accurately. The most important : te: es £ 
. Aeronautical Sciences, Vol. 18, p. 271, April, 1951 


; 


i i ' ; >} Hawthorne, W..R., Secondary Circulation in Fluid Ff 
section may therefore not necessarily be to increase the proce. Roy. Soc., Ser. A, Vol. 206, p. 374, 1951 


effects of the secondary flow at a wing-fuselage inter- 


drag of the wing, but to cause a harmful disturbance in ‘ Kronauer, R. E., Secondary Flows in Fluid Dynamics, Ph.D 
the boundary layer on the fuselage downstream of the Thesis, Harvard University, March, 1951; Proc. Ist National 
Congress of App. Mech., pp. 747-756, 1952 

5 Milne-Thomson, L. M., Theoretical Hydrodynamics, 2nd Ed 
Maemillan, 1950, p. 167 

6 Wilson, W. H., Shear Flow in Boundary Layers, Thesis for 
degree of Aero. Engineer, Mass. Inst. of Technology, 1951 


wing, to sweep and collect the boundary layer on the 
wing into regions of locally stalled flow, and to change 
the lift on the wing. An effect of the secondary flows 
on the boundary layer is suggested by the considerable 


discrepancy between the measured drags and the i Ling, A. T., Study of Secondary Flow Effect on Single Sy 
measured or calculated energy in the secondary flow. metrical Airfoils, M.Sc. Thesis, Mass. Inst. of Technology, 1952 
Appendix 


EVALUATION OF DEFINITE INTEGRALS FOR 61° (BICONVEX PROFILE) 


The expression for 61 for the biconvex profile is given in Eq. (7.9) 


i: aU ¢€ ™ n* sinh? V/n ; 
| a i-— —— d (coth Y/n) + 
dz #- Le 16 sinht Y/2 
: n* ! cosh? V/n af 
1 — }d (tanh V/n 7.9) 
J 0 16 [cosh VY — cos n(2/2) |? 
Writing x = Y/2 and y = 2/n the first integral becomes 
” 1 sinh? yx | | *™ sinh? yx 
t= d (coth yx) = -—])—- ——— 4x A.1) 
‘ y* sinh! x 76 y* Ji sinh? x 
where 6 > Osothat coth yi—>1/y5. Consider f° (e°’*/sinh‘ z) dz, where s = x + iy taken round the rectangular 
contour shown in Fig. 16(a) with the corners of the rectangle ABCD at R,71; —R, wi; —R, —7i; R, —7i; re- 
spectively. The contour is indented at 0, +77 where the radius of the semicircles ¢; and ¢» is 6. 


Now | <  < 2and hence 2 > y > 1 sothat along AD as R—> ~ 


272 ° { 2 1)R 
e sinh! s ~e — () 


Along BC 


dae 


e~’* sinhts > e — 0) 


Along AB, but omitting contour ¢, 


. - z } e- x + m1 t e~? t ] 
; dz = — = dx — - dx (A.2) 
n sinh? z Js sinh? x Fa sinh? x 


Along CD, but omitting contour ¢» 


eu ev . e x rt . en x 7t 
- dz = aaa dx + ae dx (A.3) 
sinh? g J sinh? x Jt sinh? x 


For the contour c, write z = mi + ¢ where ¢, = 6. Then the integral along this contour becomes 


ee "ig .; 
ihe : d¢ 
a sinh* ¢ 


and this may be expanded to give 


Cascade of Aerofoils in a Non-Uniform Stream, Journal of the 
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2 | — 2\1 { 
=¢ — + 2 (2? — — — — y(y? — 1) mi Ad 
> 6 3 6 3 
For the contour c write z = — ai + ¢ and the integral becomes 
— ev - 
=e : —— 
i sinh? ¢ 
2 | —  2\ I { : 
=e — —2{27° - emia. % aie (A.5 
> Oo } 0 
The integral round the contour becomes from Eqs. (A.3), (A.4), and (A.5) 
aa *“ cosh 2yx 7 2 | 1\ 1 
—41 sin 2y7 ae as + 2% an 20) - — + stir - —}-—i- 
J5 sinh? x 3 68 3/ 6 
| ; ee 
> cos 2y7 v(y* — 1)at = 2mt X residue (A.6 
3 
Now by inspection of Eq. (A.4) there is a pole of order 4 at the origin, Fig. 16(a), and the residue is (4/3)y(7*° — 
1). Hence Eq. (A.6) becomes 
-_ 2 sinh? yx + 1 — 2 |] | l 
—4i sin 2y7 - dx + 2% sin 2yr | —- — + 4[7’- “ie Ths 
J 5 sinh* x 3 6 3/ 6 
_ 4 
2at— ly" — 1) (1 + cos 2ye ALi 
o 
Now 
; dx l - 2? l l | 0 ‘ 
: —- = | Cot a2 — — coll’ x PS es a _ : —., .) 
Js sinh *x 3 3 6 3 \6 } 
2 l 1 | ; 
= — > —— : | —— @ 
3 0 > Oo 
Zz z 1 1 | 
= — : + (A Q 
) > O + 0 
omitting terms in 6 and powers of 6. Whence Eq. (A.7) becomes 
sinh? yx . 9 i l Tt 
: ax = — — _ y\¥ _ l cot Tr A 9 
‘ sinh? x 6 3 3 
which, when substituted in Eq. (A.1) gives for the first integral as 6 > 0 
l ge (7° — 1] n| | n wiz n 2a 2 
+ cot yr — 1 = + _ cot - (A.10 
33 38 2L3 \2 ae 64 - 6% 
Writing x = Y and \ = 1/n the second integral in Eq. (7.9) becomes 
. l cosh* Av ] : cosh® Ax : 
= ; d(tanh Ax = dx A.1] 
0 1.4 [cosh x — cos (/2A) | tr? JO cosh x — cos (2/2A 
ular contour ABCD shown in Fig. 16(b The 


Consider f e~’ dz/[cosh z + cos (7/2) ]* taken round a rectang 
corners of this contour have the same coordinates as in Fig. 16(a 


l along ADasR—> o 


Since 31/2 <A < 
2re 
e 

—> @ 
[cosh 2+ COS (fF 2X) ]? 


and there are no indentations 


—> () 
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Along BC 


2Az 


e 2(1 + A)R 


_— — 0) 
[cosh z + cos (2/2d) }? 


Along AB 


eo e* LT wt 
———— d= — ff 5 e“"z 
te [cosh z+ cos (#/2A) |? 0 [cos (7/2) — cosh x]? 
. 2A(x 4 1 
e 
/ ae (AT?) 
Jo  |cos (27/2) — cosh x]: 


Along CD 


* + 2Az . 2A(2 t 4 x ri 
J ae | a / a 
= - ax + dx (A.13) 
[cosh z + cos (7/2A) |}? 0 [cos (7/2) — cosh x/° Jo [eos (7/2) — cosh x]? 


Hence the integral round the contour ABCD becomes 


et ~ cosh 2Ax ; : f 
— 41 sin 2\7r - dx = 27ri X sum of residues 
J0 [cosh x — cos (2/2d)]? 


Whence 


5 cosh? Ax w/4 7 l dx 
: : = — X residues + — (A.14) 
0 [cosh x — cos (1/2A) ]? sin 2\7 2 Jo [cosh x — cos (#/2A) |? 
The value of the integral on the right-hand side may be determined by substitution of y = e* and reduction to the 
form 
] fy dy a{1 — (1/2A)] cot (4/2A) + 1 5 
<< an = (A.15) 
2 l = 2y cos (9/2) + 1}? 2 sin? (m/2A) 


the result being readily obtainable from tables. t 


If the poles of e’/[cosh z + cos (#/2A)]° occur at s = a, then since cosh a = —cos (1/2) poles occur inside 
the contour at 


a = mi{l — (1/2A)] and a = —7ill — (1/2d)] 


To find the residues put ¢ = s — a then 


rz 2r\a 2Xdf 
_ > 


e~* (cosh z — cosh a)~? = e € (cosh a (cosh ¢ — 1) + sinh a sinh ¢)~? 


giving on expansion 


I 


+.. .) cosech? a 


i. é 


a ( l (2\ — coth a) 
” i + rs 


Hence the poles are of order 2 and the residues are 


wi(2h — 1) Fe . ig 2 /. /é wi(1 —2d) [ « oa v T 
— e”™ [2\ — i cot (2/2d)] cosec? (r/2A) and — e™ ‘ (2 + 2 cot 5 Jeoseet “ 
Zi 2r 


Hence the sum of the residues is 


[4X cos 2Am + 2 sin 2m cot (2/2A)] cosec? (x/2d) (A.16) 
Substituting Eqs. (A.15) and (A.16) in Eq. (A.14) 
ao cosh? Ax dx me eres ie ie ae -_ a 
| alae = osteo? ae” cosec* (m/2d) E —T (2 cot 27 + nH cot (x 2») (A.17 
Eq. (A.11) then becomes 


nm V2 n \* , un 27 2r nT nn\ |) 
— _ cosec? l— cot + cot (A.18) 
2'yn 2 2 n n = > f 


+ Dwight, H. B., Tables of Integrals and Other Mathematical Data, Macmillan, 1947. 
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Supersonic Flap Lift Effectiveness for Some 
General Plan Forms 


RICHARD C. FROST* 
Convair-Fort Worth Division of General Dynamics Corporation 


SUMMARY 


rhe application of the concept of reverse flow to the solution of 
the supersonic flap lift effectiveness is presented hose solutions 
which could be obtained by the use of conical pressure fie'ds are 
given including solutions for both the case of a supersonic trailing 
edge and a subsonic trailing edge. The value of these solutions 
lies in the fact that a good correlation with test data can be ob- 
tained provided a theoretical correction is made for the profile 


ngle at the trailing edge. 


SYMBOLS 


= aray through the origin of the coordinate system and 
equal to By/x 
= the value of 8 tan e, where ¢ is the angle between the 
free stream direction and the ray ‘“‘a’’ which passes 
through the juncture of the hinge 'ine and the inboard 
flap edge. 
1 = aspect ratio of the wing 
= wing span 
= s/ b/2) 
i = pressure coefficient 
7 = Mach Number 
= 6 an €zs 
= Stan ey! 
Wy = § tan er; 
span of flap on one wing semispan 


= wing area 


,¥ = coordinate axes (s e Fig. 1) 
a = angle of attack 
3 = iP 1 


= flap deflection measured normal to the hinge line 
€ = 90° —A 


A = sweep angle 


S ibscripts 

HTL refers to’ the hinge line 

LE refers to the leading edge 

TE refers to the trailing edge 

1 refers to wing in direct flow 

r refers to wing in reverse flow 

(refer to upper limit of integration 

10 refer to various regions of conical pressure distribution on 


the wing 


INTRODUCTION 


A GENERAL APPROACH to the solution of the super- 

sonic lift effectiveness for a trailing edge flap on 

an arbitrary wing plan form has been presented by 

Frick,! although he does not carry out the integrations. 

Goin’ has presented the flap characteristics in equation 

and chart form for the case with supersonic edges. 
Received November 23, 19538 


Senior Acrodynamics Engineer. 
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Frick’s approach involves the concept of a swept line 
source along the hinge line with corrections made at the 
edges of the plan form to cancel the lift off the wing as 
needed. The method yields closed form solutions ex 
cept when the trailing edge becomes subsonic. Here 
the corrections at the plan-form edges affect other edges 
of the plan form and consequently, the corrections must 
be corrected ad infinitum to obtain a complete solution. 
An application of the reverse flow concept overcomes 
this difficulty and a closed form solution is readily ob- 
tained. In addition, it is possible to build up solutions 
for the case with a supersonic trailing edge, as well as a 
subsonic trailing edge, by a straightforward application 
of conical solutions. It is the intent of this paper to 
build up solutions for flap lift effectiveness for both the 
cases of subsonic and supersonic trailing edges by an 
application of the reverse flow concept and known or 


easily obtained conical solutions. 


REVERSE FLOW CONCEPT 


The reverse flow concept has been quite extensively 
discussed by Heaslet and Spreiter? and the application 
to a wing with flaps is indicated. From Eq. (61) of their 
paper the equation to be used here is presented: 


| ag i A ), 
i = | ( - )as, 
0 cos Ay; sy a, 


In coefficient form this becomes: 


where: 

C,,/a@,) = pressure coefficient per unit angle of at 
tack of the flat plate wing (flap unde- 
flected) flying in reverse flow 

Sy = portion of the wing actually covered by 


the flap 


For wings flying in reverse flow the application of 
available conical solutions for (C,/a@) is all that is re- 
quired to include a very broad range of wing-flap con- 
figurations. Thus, the flap lift effectiveness may be 
readily obtained for the wing-flap configuration flying 
in direct flow by an integration, over the portion of the 
wing covered by the flap, of the pressure coefficient per 
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unit angle of attack of the wing alone flying in reverse 
flow at the same Mach Number. 


CONICAL PRESSURE COEFFICIENTS 


The pressure coefficients per unit angle of attack that 
will be required for this discussion are located in the ten 
regions as indicated in Fig. 1. Their values are as 


follows: 


c a= 


mM, 
cos! 
— | 


12,/BV m,- — | 


fi, +a— ome) 
m, — a 


= 
a/s2 7r3V m," 


201 —-Vl1l— 
, = \ 
E(k) + nkim(— m2, k)! fmm, 
' ‘ ' 
ar /2 
E(k) = | V1 — k* sin? 6 dé 
7 
m/e dé 
T)(— my", R agent aa 
Jo (1 — m," sin? @) V1 k? sin? 6 
and 
ki = V1 — k?: Mm = V1 — n? 
2V (1 — m*)(1 — mm, 
R= : 
l1— mm, + V (1 — m?*)(1 — mm; 
m+ m, 
i = 
l+mm,+ V (1 — m*)(1 — m,*) 


(C,/a)5 and (C, 
Haves, Browne, and Lew.* 


a@)j Were obtained from a paper by 
The remaining values were 
derived following a method outlined by Hayes and Lin- 
stone.* 


FLAP LirT EFFECTIVENESS SOLUTIONS 


In the analysis of flap effectiveness, only those flaps 
extending from the wing tip inboard will be considered. 
The effectiveness of a flap located arbitrarily along the 
semispan of the wing may be easily determined by sub- 
tracting the effectiveness of a flap extended to the out- 
board edge of the desired flap from the effectiveness of a 
flap extended to the inboard edge of the desired flap. 
The fact that this is correct is very easily deduced from 
Eq. (1) where the integration is only carried out over the 
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where /(V 1 — m,’*) is the complete elliptic integral of 
the second kind with modulus V 1 — m,”. 
CaF ay = (oF Q)} 
ce tm, 
= : x 
aly r3B Vm," — 1 
(1 + m)(m, + a) — 2(m + am,) 
COS 1 
(1 — m)(m, — a) 


(=) | ( +m — ~*) 
= cos ~* 
a/sy TO l—m 
(2), - S74 
a/; a re] % m,—a 


where: 


m,*)\(1—- Vil — m- 
— (1 — V1 —m,) (1 — V1 — m?)! 
is of the following form 


actual flapped area. Eq. (1 


since we are dealing with conical pressures: 
C; l 


Lé ze ' a 
cos Ay, S | Js; @ 


This expression may be put into the following two 


(a) dx dy 


forms suitable for the integrations required by means 
of a transformation of the variables of integration using 
the appropriate Jacobian as indicated; also a factor of 
2 may be included if only one semispan is considered: 


Cis 2 ft: psa. x9 
= — (a) | J dx da 
cos Any S eo a - F a x, @ 
where: 
I(x, y)/(x, a)]| = x/6 
Cr; 2 fade. Xu 
= — (a da (2) 
cos Ay) Ila, © 28 
Similarly, 
Cis 2 fm fx C, cae 
= f | {a J dy da 
cos Ay) Sc i a y, a 
where: 
J [(x, y)/(y, a)]) = By/a 
Cis 2 fm ( By," . 
= (a) da (9) 
cos Avi p » a 2a 
where: 
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eS a SWEPT TAPERED PLANFORMS 
Since By = ax,, | - > Y 7 = 


Lg 3 2 


om) Ros foe i. fe 
5 \ SOS wer 


mo) y,= 5 SUPERSONIC TRAILING EDGES 


For the swept tapered plan forms a) = O and for the | 
‘ 5 . 
If the flap hinge line 





pointed tip plan forms ad) = m. ! 4 
intersects the wing tip at the trailing edge aj = my. “< 
: 5 ; es Pee M ss 
Where the functional dependence of C,/a@ on “‘a \ 
changes in the interval of integration in either Eq. (2 “a 
or Eq. (3) the integral is broken up into two integrals SUBSONIC TRAILING EDGES 





with the upper limit on the first integral and the lower 


limit on the second integral set equal to one (1). The POINTED TIP PLANFORMS 


complete solution for the flap lift effectiveness is then 


the sum of Eqs. (2) and (3). 
\9 ff 


Swept TAPERED PLAN FORMS ‘ \ 
, \ 
The solutions for the flap lift effectiveness are now 


readily obtained. Each solution has natural limits SUPERSONIC TRAILING EDGES 
based on the position of the Mach lines and the sweep of 


the hinge line and the trailing edge. The solutions are 





conveniently listed according to their limits as follows: 





orm 
I) Supersonic Swept Trailing Edges 
(A) a, > 1.—The Mach line from the origin of the 
a er nul SUBSONIC TRAILING EDGES 
coordinate system intersects the flap hinge line. : ; 
: : Fic. 1. The general plan forms in reverse flow and the ten re 
(1) Mp, ~ M;,: gions of conical pressure distribution 
-wO 
ans 
ing Cy m,(b,/b)? ‘m, — a,\* m, {My — ay m,(m | 
of eS = 
i A cos Ay i ae Aft oe 9 ay mM, ay m,(m, — 1 
re 
(14) m, = m,: 


Cr (m, — a,)(b,/b)? 2 m, — a, | 


A cos Ay; a: Wat — 3 
(B) a, < 1—-The Mach cone from the origin intersects the flap inboard edge 


(2) m, ~ m,: 
2) (a) m, > 1. 


F 
Cis (b,/b)* 2 m,(m, — 1)(1 — a, My (™ - wy ; (™ a ei 
= < _ cos ?* 
A cos Ant oon? 49 \ a; (m, — m) a; m,—-a 
mM, My, — Q,\" m(m, — | mM, + a, — 2m,a 
( coe * ( ) 
(mM, — Mm, a, ni,(m, — 1 My, — Qa; 


(b) m, < 1. 


/ 
Crs (b,/b)° Imm, — 1)(1 — a, My, mM, — a,\° fm, + a; — 2m, 
= Z ‘am cos 
| a, (m, — My, a, mM, — a, 


:) A cos Ay; aVm — 


mM, mM, — a,\* mm, — | my, + a, — 2m,a 
cosh! . 
(mM, — Mm) a; m,(1 — m, mM, — a 
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(44) m, = m,: 


II) Trailing Edge Normal to the Free Stream 


For this case m1, = 
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II) Subsonic Trailing Edge 


While the integrations are more straightforward and 


(A) a, 1.—(1) m, ¥ m;,: ee é; ; aie 
” h some saving in labor is realized by using the principle of 
Crs (" 2) (m, — a,)? | reverse flow for the problems with a supersonic trailing 
mM, — “ee . as . oF 
.dge ‘ase with ; tr: s edve t] 
A cos Ant a, Viet. ~ DJ edge, it is lor this ¢ Ase with 1 subsonic trailing edge that 
a substantial saving in labor is realized. In addition, a 
(44) m, =m, = « closed form solution is obtained which is not feasible by 
’ direct methods. 
C; F) re ha, — ] 
- (b; 0) 
. . . ») 2 
A cos Ay; 2Q, 1; m, ~m 
(B) a, < 1.—(t) m, #~ m,: (a) m,> 1 
; C ; (hb, h  H21 — WV | — H, 
C1; (b,/b)? = x 
= > | 24, Val — a,) A cos Ay 2a,(m, — m,) ;E (V1 — m,*)} 
A cos Ay "a," 
i ee 9 ! oe: a 
- (m, — a, = ma, — 2m,a, .|m, — 
My, COS ( — i = 
7 mM, ay 
V m,(m, — 1 " 
es m,(m;, — a,)* mm, + m,a, — 2m,a 
, (Mat a — ant) cosh-! 
cos” 9 / a ; MAM, — Ay 
ie gs ) 2a, V m),(m, My, ( 
(6) m, < 1 (B) m =m 
Cr; (b,/b)? ; ; 
- = 24, Vail — a,) + . 4(b,/b)2 ¥2 — 
2 t )7/D) 2.1—-V1l—~m, jm, 
A cos Ayr 1 Crs a f \ I ¢ 
(m, — a,)? A cos Ani 3m, | EX Vil—m/ } \ a 
m, cos! (1 — 2a,) — 4 
Vm~,(1 — my, 
‘Mn + a, — 2mpa, POINTED Tip PLAN FORMS 
cosh7! 
Mn — Q, ee ; : ’ 
Special consideration must be given to pointed tip 
(11) Mm, = m, = & plan forms because of the influence in the tip region of a 
, ; | supersonic or a subsonic leading edge. 
C15 (by b)? | (2a, + 1) jl — a, 
A cos Ay 7 T a; y ay ~~ : 
I) Supersonic Swept Trailing Edge 
(4a, — 1) 
cos~! (1 2a : : 
2a" (A) a, > 1.—(t) m, + m,: (a)m, ¥ a; m> 1. 
s a ! ae { ] - : ( [ (on » | 
C1;/A cos Anz = [(b,/b)?m,/a,V m2 — 1] [) [(mm — a,)(a, — m)|/aitn — my + [Om — a,)/m, 
(Ob) m = Ge my < 1. 
| 
Cr m,(b,/b)? (™ — =) My (“" — a ) (m, — m)(m, — 1 
A cos Aj; (ee ee a mM) a / (772, — m)(m, — 1 


(¢ Mp = Q;y. 
C;,/A cos Agr = [(6,/8 
(17) my m,: (a) m> 1. 


Crs (m, — 


A cos Ant a, ~ 


{ 
Vm, — 1] $[(m, — m,)/m, 
, ' 
a, (by bh ss nuim, — a,) 
m2 — 1 a,(m, — m) 
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m< | 
( n — ad, b, D } — 2,)(n + pi) 2m 
) 
1 cos A a t/a? — 1 a 2a,(m, — m)(n 
; “ (2), FH ( /} = a }} > l 


Wl ~ }} l 
Crs (b,/b)? m1) ” re 1 + m)(m, +a 2(m + aym 
e COS 
A COs Any rN m,- — | Hl — es a (1 — m)(n - a 
| ) 
Ho, m)(m, — 1 l + m)a, — 2n ” 
COS - 
m (l — ma (m 
l(m, — m)(m, — + m)(m, +a 
cosh 
(mM, — m)(1 1) | mi 
(¢ mM, =a 
( (b,/b m,— m (1 + m)(m, + m, 2(m + my, 
) cos 4 
1 cos Ay we A n 1 — m)(m, — My) 
lim mi) (mM - | fb ale 2 
COS 
n l 
(12) my, = mM: 
Ce b,/b m,(m, — a ha — a,)(a, — Hn 
A cos Aur T V m,- -_ l a,’ Vin, =— 3 } nos l 
| ) 9 
(m, — m)(m, — | Jl + ma, — 2m m,(m, — a,) | m, +a m, + 
cos™ T 2 - T 
m (1 — m)a a,(m, — | ' ma 2a 
(1 + m)(m, + a 
cos~! 
(1 — m)( 
II} Trailing Edge Normal to the Free Stream (m, = 3) m, wm. = wo: (a) mm > 1. 
(A) a, > 1.—(1) m, ~ m,: (a) m ~ ay m> 1. : 
( h, s 
— (* 4 j 2a,m, — mm, — 4," 1 cos A | ar ca 
A cos Ay a, mM, — M1 
h) m | 
mM, A~ a,~m l 
? . S } } 
C. (* 4 aoe ( (~~ ) ha 
My — - 
A cos Ayy a, 
A cos Ani ay V (m7, — m)(m, — 1 
m1 a These last two results agree with the results obtained 
h fe 
by Tucker and Nelson. 
® (b, a) : 
A cos Ay; M1) (B) a | 1) m, #Nn (a) n 
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Cis (b,/b)* | na, 1+ m — 2a, | 1+ mj)a, — 2m 
= : COs I COS == 
A cos Ay; T ay I m Vm (1 — m)a, 
(My — az) —f(1 + m) (mp, + a,) — 2(m + aym | 
COS ) | 
a," V (nm, — m)(m, — 3 (i = 90h, — a, 4 
(b) My aE Mm, < 1. 
Cz (b,/b)* | my, 1 +m — 2a,) | 1+ mja, — 2m 
~ , cos ) t cos — 
A cos Ant T ad; l1—m Vm 1 — mya, / 
(1, — a,)* (1 + m)(m, + a, 2(m + a.m, | 
cosh , 
a,” V (m, — m)\(1 — m (1 — m)(m, — a J 
(C) My, = Ay. 
Cy (b,/b)? | f1+m — 2m, | 1 + m)m, —2m 
= COS ‘ + COS 
A cos Ayrz a My l—m v/ mn (1 — m)m, 
(12) my, m, 
Cx, (b,/b)? | V (a, — m)(1 — a, l (1 + m)a, — 2m 
; r cos! 4. 
A cos Ay; T a; Vm (1 — moja, 
la, —1—m 1 — 2a, 4 m\] 
cos" )| 
a l1—m | 
(III) Subsonic Trailing Edges 
(A) m, ~ m,.—(1) m, ¥ a: 
Crs ; | m,(m, — a,) (m,— m  f{ma, + ma, — 2mm, 
= F(b,/b r V (a, — m)(m,— a,) + cos - 
A cos Ay; a,"(m, — My) 2V mm, (mm, — m)a, 
(m, — m)(m), — a,)*m, — f(m, + a,)(m, + m) — 2(m,a, + mm 
cosh 
2a,7(m, — My) V (my, — m)(m, — Mm, M, — M)\\My — A, Za 
_ 
(4t) ity = a): 
Cis ; | V (Mm, — m)(m, — my) (m1, —m m,M, + nM, — 2mm, 
= F(b, b) - cos ~ 
A cos Any My, 2V mn, (m, — m)m,, 
(B) my = m,, 
Cr; ; 5m,a, — 3ma, — 2mm, 
- = P(6,/b sess (a, — m)(m, — a,) + 
A cos An; 3a,°(m, — m) 
m,—m {mas + Ma, — 2mm 
cos” 
» . mim, (Mm, — m)a, 


DISCUSSION 


The preceding expressions are strictly valid only if 
the Mach line from the apex of the flap in direct flow 
does not cross the opposite semispan. The region of 
validity is given by a, > (b,/b)/(2 — b,/b). However, 
any carryover of loading to the opposite semispan ac- 
counts for only a very small contribution to the flap 
effectiveness if the affected area is not too large and 
consequently good approximations to the flap effective- 
ness for such cases is still obtained by using the pre- 
ceding values. 

The type of flap considered here is operating in a 
region of the wing where the local slopes of the surface 


have a marked effect on the loading. The discrepancy 
between test data and the results of the preceding 
analysis can be attributed largely to the profile angle 
at the trailing edge of the flap. Correction factors have 
been obtained by comparing the linearized theory two- 
dimensional pressure coefficients with the pressure co- 
efficients obtained by using the Busemann higher order 
approximation to account for thickness.* 7 This form 
of correction to the linearized theory results has been 
modified and extended empirically*® giving a very good 
correlation of test data. 


Continued on page 620 
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Optimum Design of Indeterminate Frames 


WILLIAM B. READEY* 
McDonnell Aircraft Corporation 


SUMMARY 


rhe paper presents methods for obtaining the optimum sec 
tion property distributions, for indeterminate frames and 
beams, directly from the design conditions. An optimum sec 
tion property distribution is assumed to be one which results in 
constant extreme fiber strains, where the tension strain is not 
nec ssurily assumed to be equal to the compression strain As a 
consequence of the constant bending strains, the change in 
geometry, due to a given strain, is determined when the inflec 
tion points are located 

With the inner and outer frame contours assumed to be known, 
i method is presented for determining the optimum inflection 
point locations, and consequently the section property distribu 
tion, for a frame (either symmetrical or unsymmetrical) loaded 
either symmetrically or unsymmetrically; and also, a method is 
presented for determining the inflection point locations, which 
result in optimum design, for a symmetrical frame loaded sym 
metrically by any known or unknown system of loads which 
yields exactly four inflection points. With the above qualifica 
tions, the optimum inflection point locations are shown to be in 
depe ndent of the loading condition In general, the solution is 
whieved by an iteration process, but a direct solution is pre 


sented for certain beam problems. 


INTRODUCTION 


SUALLY AN INDETERMINATE structure is designed 
U by arbitrary methods and then checked for 
strength, with no assurance that the resulting struc 
ture is efficient. This paper presents a method for ob- 
taining the optimum section property distribution for 
indeterminate frames directly from the design condi- 
tions. 

With the inner and outer frame contours assumed to 
be known, the following types of problems are con- 
sidered . 

Type I.—-Determine the inflection point locations 
which will result in optimum design for a symmetrical 
frame loaded symmetrically by any (known or un 
known) system of loads which yields exactly four in- 
flection points. It will be shown later that with the 
above qualifications, the optimum inflection point loca 
tions are independent of the loading condition. 

Type I7.—Determine the optinium section property 
distribution for a frame (either symmetrical or un- 
svmmetrical) loaded either symmetrically or unsym- 


metrically. 


DESCRIPTION 


Fig. 1 shows a frame and the same frame as it is 
idealized to facilitate analysis. The idealized frame is 
Received December 2, 1952 Revised and received January 
l?, 1954 
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assumed cut and fixed at one end. Unknown loads 
My, Pox, and Pyz are assumed to act at the free end. 
Any stress condition which existed in the frame before 
it was cut can be duplicated if these unknown loads are 


chosen correctly. 


DEFINITIONS, SYMBOLS, AND CONVENTIONS 


Optimum Design. For the purposes of this analysis, 


“optimum design”’ is defined as a design which results 
in constant strains in the extreme fibers of the compres 
sion and tension caps, but where the tension strain is 


not necessarily equal to the compression strain. 


H = frame beam depth 
= inflection point number 

V7 moment positive for compression in 
outer fiber 

VW = static moment in the idealized frame 

VM, = Pox = P 0) 

Wy, Pox, and Pyz = unknown loads 

S = clockwise distance from free end along 
assumed centroidal axis 

\ = horizontal distance from free end (posi 
tive to right 

Z = vertical distance from free end (positive 
up 

7 = total strain divided by the absolute value 
of the total strain 

€ = strain (tension strains positive 

Eo, € = strains in the outer and inner fibers, re 
spectively 

er = & — € = total strain 

Ay, Ax, Az = parameters (rotation or translation pet 


unit total strain 


LIMITATIONS OF METHOD 


Because only the bending deflections are considered, 
the usual limitation exists, that the ratio of // to the 
overall frame dimensions should be small. 

The slope of the stress-strain curve is assumed to be 


positive. 








: ncimianitiaba 
Poz 
' 
————_|] —__ » Po, —= 
Mo 
Fic. 1 Typical idealization of frame to facilitate analysis 
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ANALYSIS 
The total moment at any section is 


My = M' + My — ZPox + XPoz (1) 


If the idealized frame, shown in Fig. 1, is to be a 
true representation of the actual frame, the free end 
should not rotate or translate relative to the fixed end. 
The strain, ey, must be constant between adjacent in- 
flection points, 7 and 7 + 1, for optimum design; there- 
fore the rotation of the free end is 


d§ N “i+l d§ 


, = 
I] 

where ¢; is the absolute value of er, and er = yer, 

viz., y = +1.0 and has the same sign as J/7p. 


Eqs. (3) and (4) give the translation of the free end 
in the Y and Z directions, respectively. 


_ vs \ ~~ +1 nae. 2 
$ €r dS = |e >, / y is = 
IT jul Ji I] 


erAx (3) 


f x 
€r aS = 
: W 


Thus, if Eqs. (2), (3), and (4) are set equal to zero, 


> J = ee \ 
ar d. = "| { ) 
€7 = J; Vy €7 Z 


Eqs. (5), (6), and (7) result. 


Aw = 0 (5) 
Ax == §) (6) 
Az = O (7) 


APPLICATION 


Because the absolute value of the total strain, e7, 
is constant, the extreme fiber of the tension or compres- 
sion cap, of an optimum structure, goes plastic at the 
same rate everywhere, consequently, similar cap sec- 
tions retain their relative stiffnesses in the plastic range. 
The above equations are therefore not necessarily re- 
stricted to the elastic range. 


Analysis of Beams 


The analysis of frames, developed above, is also ap- 
plicable to many beam problems. The simple beam 
problems analyzed below are presented to illustrate that 
the assumption of constant bending strain determines 
the geometry between inflection points, even before the 
required section properties are calculated. 

PROBLEM: A constant depth beam of optimum de- 
sign, Z in. long, and fixed at both ends, is to be de- 
signed. If the beam loading is such that it produces 
two inflection points, where are these inflection points 
located? 

Solution: The beam is assumed to extend from Z = 


0 to Z = L, and the following equations result. 
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z dS 2. [- —Z 1S = 0.527 
a a a Wy’ 


The inflection points are designated as Z; and Z», and it 
is temporarily assumed that the moment is positive at 


Z = 0. 


Aw = 0 = (2 IT)(Z, - Zo + (Z 2)) 
2 L? 
Avy = 0= ( 0.52,? + 0.52 0.5 -) 
IT 2 
The solution to these equations is: Z,; = O0.25L, 
Z» = 0.75L. Obviously, the same result is obtained 
if the moment at Z = 0 is assumed negative. 


PROBLEM: A constant depth beam, LZ in. long, of 
optimum design with two equal spans, is to be designed 
with one end support fixed and the other two supports 
pinned. If the beam loading is such as to produce one 
inflection point in each span, where are these points 
located ? 

Solution: The beam is assumed to extend from Z = 
0 to Z = L, with the fixed end located at Z = L. 
The unknowns are Pox, at Z = Oand Poy, atZ = 0.5L. 
The requirement in this problem is that neither the 


0.5L 


therefore 


beam section at Z = O nor the section at Z = 
translate relative to the fixed end at Z = L; 
Eq. (6) must be satisfied from Z = 0 to L and from Z = 
0.5L to L. 


equations. 


[ yZ : [ y¥(Z — 0.5L) ' 
4 dS _ ‘i ds = J 
: HT Jost H 


rz 
/ (—Z/IH) dS = —(0.5Z?/H) 


J 0 


*Z (Z — 0.5L) 
/ ‘i dS = 
05 L TI 


— 0.52? + 0.5LZ — 0.125L? 
1 


This requirement yields the following 


« 


The inflection points are designated as Z; and Z» and, 
as in the previous problem, the moment at Z = 0 is 


assumed to be positive. 
2/H[—0.52,? + 0.52.2 — (0.5L2/2)] = 0 


») 


2 we 0.125L2 
; (0.52: — 0.5LZ. + 0.125L? — ) = 0 


9 
The solution to these equations is Z; = 0.47S8Z and 
Ze = O.8541. 

Whenever the beam depth, //, is constant, as in the 
previous problems, the beam bends in a series of circu- 
lar arcs with equal radii; therefore, the solution can be 
verified easily by geometric considerations 

Type I Prosltem. (Symmetrically Loaded Sym- 
metrical Frame with Four Inflection Points.) 
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OPTIMUM DE 


TABLE 1 TABLE 2 
ilues from Fig. 3 Values from Fig. 3 
“dS j Z ds 4 Z 
S _ ds S — dS 
H H H H 
7 3 —1 35 3.3 2 5 — ().79 
7 7 — 14.08 7.5 7.5 —16.71 
10 —2t 29 10.0 10.0 — 31.83 


Because of symmetry, only one side of the frame need 
be considered. According to the statement of the 
problem, there are two inflection points on each side. 
Eq. (7) is satisfied because of symmetry and Eqs. 
5) and (6) do not involve the loading condition; there- 
fore, the inflection point locations are independent of 
loading condition because any inflection point locations 
which satisfy Eqs. (5) and (6) can be achieved by some 
values of the unknown loads JA/p and Pox. 

The solution to Type I Problems is best illustrated by 
asimple example. PROBLEM: Determine the inflection 
points which yield optimum design for a symmetrical 
loading condition imposed upon the frame shown in 
Fig. 2. 

If inflection points are assumed at S = 3 and 7, the 
values from Fig. 3 are given in Table 1. 

If \/; is assumed to be positive at S = 0 (either sign 
can be assumed) : 

Aw = 213 — 7 + (10/2)] is positive. 1/) should be 
reduced, because \,;/.)/) is positive, as explained in the 
next section. 

ly = 2/—1.35 + 1408 — 
P,x should be increased. Both an increase in Pyx anda 
decrease in .J/) result in a decrease in .\/7; therefore, as 
can be seen from Fig. 4, the inflection points should 


(31.83/2)] is negative. 


SIGN OF INDETERMINATE 


FRAMES O17 





—)|—e 





A 
Mo 


Circular frame loaded symmetrically by a loading that 


Fic. 2. 
produces exactly four inflection points 
] } 


event that Ay, and Ax indicate that the inflection points 
should be moved in opposite directions, the integral 
curves (Fig. 3) should be consulted to determine which 
inflection point has the greater effect upon a particular 
integral. 

If new inflection points are assumed at 2.5 and 7.5, 
the values from Fig. 3 are shown in Table 2. 


Aw = 212.5 - daa + (10/2 } = U 


2(2. 
Ax = 2[-—0.79 + 16.71 — (31.83/2) ] 


The inflection points for any symmetrical load which 
produces exactly four inflection points are therefore 
located at S = 2.5, 7.5, 12.5, and 17.5. 

Type II Prosptem. (Any Frame Designed for Any 
Loading. ) 

An outline of the procedure used in solving Type II 


Problems is presented below. 


Values of 1/17, —(Z/H), X/H, and Al’/H7 are calcu- 








move in the direction of the positive moments. Inthe lated for various values of S. 
(¢ 
" 0 5 10 
lor T 1 
5 r “23 fF 
L J -50 
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Fic. 3. 





Geometrical values for frame shown in Fig 
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Mr 





: rg 


\ “4 S 





Mr curve showing inflection point shift due to decrease 
in Mr. 


Curvesot (dS/f1), | (—Z/H)dS, and 
0 0 


| (X/H) dS vs. S are plotted. 
0 


Fic. 4. 


Values of 1p, Pox, and Poz are estimated, and values 
of \/,/IT from Eq. (1) are then calculated at the points 
where 1///, —(Z//1), etc., have been calculated. 

Signs of Ay, Ax, and Az are determined after in- 
flection points have been located by plotting .1/,/J/1/. 
As can be seen from Eqs. (1), (2), (3), and (4), the func- 
tions OAyz/O.\Wo, OAx/Opox, and OdAz/Poz are positive 
in the vicinity of \ = 0; therefore, the value of the un- 
known load JJ», Pox, or Poz should be increased if the 
value of the corresponding lambda is negative and de- 
creased if the corresponding lambda is positive. This 
process should be repeated until all lambdas become 
At this point the values of A/o, Pox, and Poz 
Because 


Zero. 
satisfy the conditions for optimum design. 
the moments and the arbitrarily chosen constant strains, 
and consequently the stresses, are known, the required 
section properties can be determined. 

Each: lambda is a function of all three variables J», 
Pox, and Poz; 
supplement the above method. 
case, a frame has exactly four inflection points and the 


therefore, it is sometimes advisable to 
If, as is usually the 


moment curve slope is large at one inflection point, this 
inflection point is usually not sensitive to changes in the 
variables and can temporarily be assumed to be inde- 
pendent of the variables. The locations of the other 
three inflection points which satisfy Eqs. (5), (6), and 
(7) can be determined similarly to the method described 
in Type I Problem. The values of the three unknown 
end loads which yield these inflection point locations 
can then be determined by the solution of three simul- 
taneous equations. 

If a frame and the loading imposed upon it are such 
that both J/’ and the Z coordinate of the assumed 
centroidal axis curve can be expressed as known func- 
tions of XY, and if the number of inflection points can be 
assumed accurately, a direct solution is usually pos- 
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sible. The \ equations of the moments at the \ in 
flection points are, 


My, = 0 = M’'(X,) + Mo — Pox:Z(Xi) + Pow X 


(S 


If Eqs. (5), (6), and (7) are added to Eqs. (8), the re 
sult is V + 3 equations in V + 3 unknowns. 
Although the iteration method of solution is not so 
precisely defined as is the above method, it is applicable 
to almost all problems, and is usually relatively simple. 
To illustrate the basic technique involved in the itera- 
This 
problem also illustrates a typical discontinuity problem 


tion method, a simple problem is solved below. 


that can arise in the solution of a frame with straight 
sides. 

PROBLEM: (a) Determine the unknown end loads 
Pox, Poz, and Mo which will result in optimum design 
for the constant beam depth frame, shown in Fig. 5. 

(b) If it is assumed that additional loads of P kips 
magnitude are to be applied outward at the upper left 
and lower right corners, along a line connecting these 
corners, what should be the area of an aluminum ten- 
sion tie connecting the corners? The area of the ten- 
sion tie is to be so chosen that it will take P kips, and 
the stresses in the original frame will be unaltered by 
the addition of the P loads and the tension tie. 

Solution: (a) It is helpful to note that, for optimum 
design, the area between the \/;/// curve and _ the 


My, = 
¢ (yM7/H) dS is alwaysa minimum. This is verified 


0 line is always a minimum, viz., @ 


by taking the following partials: 





Z 
+_26-—_— 





Kips \ 
24 \ 














Fic. 5, Frame loaded unsymmetrically. 
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TABLE 3 


Functions for Fig. 5 Frame 


f° se F ol 
H / H uf H d 
S 70 ‘s 30S — 0.58? — 450 —_ 
70 ra — S00 
=| S 100 S — 40S + 2,000 O58? 
100 100 —? OOO 
) S 140) S 0.5.8? — 1408S + 7,000 
140 140 —? S00 
Od $| O.\/, I O7 | dS 
Ww JL ow "OM WH 
$ OM, dS 
Y =Ay=0 (9) 
j OM, LH 
0¢/O0Pox = Ax = O (10 
0¢/OPoz = rz = O (11) 


From Fig. 5, it is almost obvious that Pox is positive. 
It appears likely that when ¢ is a minimum, P»z is 
negative. Accordingly, if Poy and Poz are assumed to be 
5 and —1, respectively, a plot of 17,// indicates that 
Vy) must be about 20 in order to make ¢ a minimum. 
The required functions are listed in Table 3 for conven- 
The values of the lambdas resulting from the 
The 
large positive value of Ax indicates that Pox 1s too 


jence. 
first trial are indicated in the first row of Table 4. 
large. For the second trial, Pox is reduced to 4, and 
M is taken as 2 in an attempt to locate the J; = 0 
line approximately along a line dividing the 1/7/H 
curve into two equal areas. The resulting large posi- 
value of Ay indicates that J) should 
smaller. J) and Py, are taken as an infinitesimal 
larger than 0 and 4, respectively, while Poz is left un- 
The infinitesimal is added 


tive be even 


changed for the third trial. 
to avoid the discontinuity in inflection point locations 
at Vf, = Oand Pyy = 4. In trial four, Poz is increased 
to —0.5 because of the large negative value of dz 
in trial three. In trial five, only Poz is changed to 
avoid crossing the discontinuity. The result of trial 
five leaves no alternative but to reduce Pox. When 
this is done and lJ, is reduced to keep ¢ a minimum, 
\y becomes negative. It therefore appears that Pox 
should lie between 4 minus and 4 plus an infinitesimal, 


H / * ds 
0 H V \ —Z 


()} ‘) () 


—() 5.8? () » § 
— 450) ) - 30) 
308 + 450 S 54) 0 — *3() Sf os & 
S 94) 10S — 540 
— 1,650 160 30 1) 
LOOS + 2.900 160 5 — 109 1() 
—2,100 160 () 1() 
—2,100 HHO—AS 0 S 4() 


—?) 100) 0) 


viewpoint, it can be seen that it is unnecessary to make 
Ay = AM, O and Py, 
values, the frame member from S 


! because, for these 
100 to S 140 


becomes a tension member with zero moment of inertia, 


0 when 


thus, the deflection of the free end in the Y direction, 
relative to the fixed end of the idealized frame, is unim 


portant. 


(b) The strain in the tension tie 1s, 
(30 + Ay.jer')*> + (40 — Az. €7 
e = = 1.0 
20 


and the load is P = eAE, therefore, the required area is 


A, id a ef 
The lambdas, Ax., and Az.,, which are the values of the 
lambdas integrated between S 0 and S$ 70, vary 
inversely with //. It can be seen that the relative de 
flections, Ay. €r and Az, er, between the ends of 


the tension tie, become excessive for small values of //; 
therefore, in frames with small beam depths, it is 
necessary to consider the changes in moments due to 


the geometric changes resulting from large strains. 


CONCLUSIONS 


If a symmetrical frame is loaded symmetrically by a 
loading which produces exactly four inflection points, 
the locations of the inflection points which will result in 
optimum design are independent of the loading and can 
be located before the exact loading is known. 

A constant height beam of optimum design, with two 
inflection points, has these inflection points located at 
0.25L and 0.75L if the beam is built in at both ends. 

A constant height beam of optimum design, ZL in. 


viz., Pox = 4. Looking at this result from another long and built in at one end and simply supported at the 
TABLE 4 
Iteration Values 
Trial : 
No V Pox P S S p S5 O.5HXy 0. 5HXx 0 5H Xz 
F l 20) 5 —_] 10) 68 80 120 1.9 L670 = $0) 
7, 2 2 { = 38 64.7 + 43.3 
2 0-4 44 a f 7.5 65.0 100 140 +2.5 —14 — 225 
) j (+ $+ —0.5 3.7 67.5 100 136.2 0) +81 —37 
5 O+ {+ —().4 33 68 100 135 0 +105 0 
6 Oo— }— —(0.4 33 68 105 140 0) —70 0) 
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center and other end, has its inflection points located 
at O.ATSL and 0.8542 if the loading is such as to pro- 
duce exactly one inflection point in each span. 

The inflection points of a circular ring of optimum 
design, with constant beam depth, and loaded sym- 
metrically so as to produce four inflection points, are 


located at 90° intervals around the ring. 


The method presented in this paper is so basic that 
it should lend itself to many types of problems not fully 
developed in this paper. A few obvious extensions are 
listed below. 

If a frame has a tie connecting one side to the other, 
and the sign of the load in the member is known, any 
desired load can be assumed acting in this tie; and, be- 
cause the strains, and consequently the deflections, of 
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le- 
termined, the tie area which will result in the desired 


the frame are known even before the cap areas are¢ 
load can be calculated. Similar reasoning can be used 
to extend the methods to consider the change in geom- 
etry resulting from applied loads on thin frames. 
Inasmuch as J/,//7, which is approximately propor 
tional to cap load, can be calculated for any loading 
condition when only the frame contour and // are 
known, it is easy to plot approximate curves of frame 
weight vs. frame beam height, //, and thus determine 
the optimum // before the frame is actually designed 
This method might also lend itself to preliminary 
weight estimates inasmuch as it yields the minimum 
weight necessary to carry any given loading, and thus 
is also a measure of the structural efficiency of the final 


design. 


Prediction of the Characteristics of Two-Dimensional Aiirfoils (Continued from page 587 


whose location, in the low transonic regime at least, 
may well be determined largely by the boundary layer. 


REFERENCES 


1 


' Spence, D. A., The Growth of Turbulent Boundary Layers, 
ARC 14, 162, July, 1951. 

2 Spence, D. A., Growth of the Turbulent Wake Close Behind 
an Aerofoil at Incidence, ARC Current Paper No. 125, May, 1952. 

3 Squire, H. B., and Young, A. D., The Calculation of the 
Profile Drag of Aerofoils, R. & M. 1838, November, 1937. 

4 Bryant, L. W., Williams, D. H., and Taylor, G. I., The Flow 
of Air Around an Aerofoil of Finite Span, Phil. Trans. Roy. Soc., 
A, 225, 1925. 

5 Preston, J. H., The Effect of the Boundary Layer and Wake 
in the Flow Past a Symmetrical Aerofoil at Zero Incidence, July, 
1945, R. & M. 2107. 

6 Temple, G., Vorticity Transport and the Theory of the Wake, 
ARC 7118, 19438. 

7 Howarth, L., The Lift Coefficient of a Thin Elliptic Cylinder, 
P.R.S.(A) 1949, pp. 558-586, 1935 


’ Piercy, N. A. V., Preston, J. H., and Whitehead, L. G., 
The Approximate Prediction of Skin Friction and Lift, Phil. Mag., 
XXVI, p. 791, 1938 

’ Preston, J. H., The Approximate 
Symmetrical Aerofoils Taking Account of the Boundary Layer, 
with A pplication to Control Problems, R. & M. 1996, May, 1943 

Preston, J. H., The Calculation of Lift Taking Account of the 
Boundary Layer, R. & M. 2725, November, 1949 

'! Theodorsen, T., and Garrick, I. E., General Potential Theory 
of Arbitrary Wing Sections, NACA Report 452, 1933. 

2 Preston, J. H., and Sweeting, N. E., The Experimental D 
termination of the Boundary Layer and Wake Characteristics 
Simple Joukowski Aerofoil, with Particular Refe 
ing Edge Region, R. & M. 1998, March, 1943. 

13 Preston, J. H., and Sweeting, N. E., The Experimental D 
termination of the Boundary Layer and Wake Characteristic 
Piercy 1240 Aecrofoil, with Particular Reference to the Trailing Edge 
Region, R. & M. 2013, February, 1945. 

14 Pinkerton, R. M., Calculated and Measured Pressure 
tributions Over the Midspan Section of the NACA 4412 
NACA Report 563, 1936. 

% Thwaites, B., Approximate 
Boundary Layer, Aero. Quart. I, pp. 245-280, 1949 


Calculation of the Lift of 


rence to the Tra 


iS 
’ 


Calculation of the Laminar 


Supersonic Flap Lift Effectiveness for Some General Plan Forms (Continued from page 614 


Rolling moment and pitching moment effectiveness 
can be derived by a similar application of reverse flow 


concepts.” 
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The Effect of Tip Removal on the Natural 
Vibrations of Uniform Cantilevered 
Triangular Plates’ 


P. N. GUSTAFSON, W. F. STOKEY, anpb C. 


F. ZOROWSKI 


Carnegie Institute of Technology 


ABSTRACT 


Experimental results are obtained for the lowest six natural 
frequencies and the associated nodal lines of trapezoidal plates of 
uniform thickness clamped on one edge. Three series of plates 
are investigated. Each series consists of six plates and is de- 
veloped from a triangular plate by progressively cutting off por- 
tions of its tip parallel to its clamped edge, giving a clipped wing 
Two of the series originate from triangular plates of 
A triangu- 


effect 
delta-wing plan form with aspect ratios of 2 and 4. 
lar plate with a mean chordline sweepback angle of 45° is used 
Results are presented in the 


A method of treating 


as the basis for the third series 
form of graphs, tables, and photographs. 
plates of nonisotropic material is described. The use of inter- 
polation to extend the results to a broader field of plate shapes is 


also discussed. 


NOMENCLATURE 


C = root chord, in. 
] = half span of complete triangular plate, in. 
Al = amount of half span cut from tip, in. 


= natural frequency, cycles per sec. 
= natural frequency of a complete triangular plate, cycles 


per sec. 
® = 2xf = natural angular frequency, 1/sec. 
w = 2rfo 
wp = dimensionless angular frequency 
h = plate thickness, in. 
E = modulus of elasticity, lbs. per sq.in. 
p = mass density of plate material, lb. sec.?/in.‘ 
= Poisson's ratio 
D = flexural rigidity = Eh°/12(1 — yu?) 
a = scale factor 


INTRODUCTION 


AS AN EXTENSION to the experimental investigations 
that have been made of the natural vibrations of 
cantilevered triangular plates! this study was under- 
taken to explore the effect of cutting off portions of a 
triangular plate parallel to its clamped edge. The 
vibrational characteristics of the trapezoidal plates 
formed by such cuts appear to be of as much impor- 
tance as those of triangular plates, since uniform flat 
plates of both plan forms approximate basic wing and 
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tail surfaces on many guided missiles and some high 
speed aircraft. 

The variations in plan form of trapezoidal plates are 
seemingly endless. Yet, if one is willing to confine at- 
tention to the case of parallel root and tip chords, 
any shape may be regarded as a function of three con- 
venient parameters: aspect ratio of the plate (ex- 
tended to form a complete triangle), sweepback angle 
of the mean chord line, and ratio of the span removed 
from the triangle to the original span of the triangle. 
If, for any values assigned to two of these parameters, 
it were desired to examine five different cases of the re- 
maining parameter, then 125 plate shapes would be re- 
quired. Such a program while not impractical, seems 
overly ambitious as long as the application of flat 
plate vibrations to aircraft and missiles is still in an ex- 
ploratory stage, and while considerable attention is being 
given to finding analytical solutions. 

This report covers eighteen trapezoidal shapes de- 
rived from only three triangles. In a later section it 
will be shown how interpolation may be used to cover 
the field of delta wing shapes from aspect ratio 2 to 
aspect ratio 4, and for an aspect ratio of 4, interpolation 
may again be used to explore mean chord sweepback 
angles ranging from tan! 0.5 to tan! 1.0. 

The configurations and designations for the three 
series of plates investigated are given in Figs. 1, 2, 3, 
and 4. Series I and II are both developed from tri- 
angular plates of the delta wing shape (the trailing edge 
is a straight line from wing tip to wing tip). The root 
chord in Series I is equal to the total span while in 
Series II it is half of the total span. Series III origi- 
nates from a sweptback triangular plate with a root 
chord equal to half its span and a mean-chord sweep- 
back angle of 45 The parameter A//l, (see Fig. 1) 
varies over the same range for each series, from 0 to 
0.4. All three triangular plates from which the three 
series are developed have been investigated previously.’ 

The test procedure and equipment used for the ex- 
perimental determination of the natural frequencies 
and node lines was the same as that described in refer- 
ence 1. The plates for each series were cut progres- 
sively from the steel plates A-1, A-2, and S-5 used 1n the 
triangular plates. 


authors’ earlier investigations of 


Thus the material properties were kept fixed for a given 
») 


series. These properties are listed in Table 2. 
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| ree eS Fic. 2 Plate configurations for Series | 
Leterme 
Fr Plate nomenclature 
TABLE | 
Measured and Dimensionless Frequencies for Tip Removals of from 0 to 40 Per Cent of Span 
SERIES I 
Family C D 
Mode 2 3 kL 
Al/ko f wh f Wp f Wp f Wp f wh f Wh 
0 32.86 22.0 91 58.7 | 16, 110 | 181 119 | 283 186 | 348 228 
ol 34 22.7 93 60.0 | 179 120 |} 181 119 | 293 «#192 | 352 231 
02 38.5 25.8 97.6 63.0 | 184* 121 | 212% 142 | 302 198 | 362 238 
+2) 41.9 28.0 99.4 64.1 | 186* 122 | 235%* 157 | 30h 200 | 366 240 
3 48.3 32.3 | 103.4 66.8 | 190* 125 | 266%* 178 | 308 202 | 379 249 
035 53-7 35.9 | 107.4 69.4 | 196* 129 | 299% 200 | 314 206] 404 265 
4 60 40.1 | 112 72.3 | 202* 133 | 350** 234 | 312 205 | 436 286 
* Family D 
*%* Family C 
SERIES II 
Family B C D 
Mode 2 a 4 
Btflyo f wp f wp f wD f Wp f wp f wp 
0 34.5 5.87 | 136 23,8 | 190 32.4 | 325 56.1 }441 76.0/578 99.7 
ol 37 6.30] 142 24.8 | 198 33.7 | 335 57.8 | 482 83.1/583 101 
02 42 7.15 | 153 26.7 | 223 38.0 | 364 62.8 | 561 96.9/598 103 
025 46 7.83 | 157 27.4 | 243 41.4 | 37h 64.5 1596 103 |621 107 
3 50.5 8.60} 141 28.2 | 268 45.6 | 385 66.5|606 105 |660 114 
035 56 9.53 | 169 29.5 | 300 51.1 | 410 70.8 |629 109 | 695 120 
4 64, 10.9 | 175 30.6 | 339 57.7 | 434 74.9 |639 110 | 718 124 
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DISCUSSION OF RESULTS 


Series I 

Fig. 5 shows the photographs of the node lines for 
the first six natural modes of plates I-1, I-2, I-4, and 
I-6 arranged in order of increasing frequency from left 
to right. For comparison the node lines for the first 
six natural modes of the triangular plate from which the 
series was developed are also shown under the plate 
designation I-0. The frequencies corresponding to 
these modes for all five plates shown, as well as those 
not shown, are given in cycles per sec. and in a dimen- 
sionless form in Table 1. A discussion of the material 
properties used in computing these dimensionless fre- 
quencies will follow in a subsequent section. 

By choosing a given mode in Fig. 5 and observing the 
behavior of the node lines as the parameter A///, in- 
creases, it can be seen that the changes in mode shape 
from plate to plate follow definite patterns, thus estab- 
lishing a type of shape family. An exception to this 
case in Series I occurs in the 3rd modes of plates I-0 
and I-1 which belong to the shape family represented 
by the 4th modes of plates I-2, I-4, and I-6 and the 4th 
modes of plates I-O0 and I-1 which belong to the same 
family as the third modes of I-2, I-4, and I-6. In 
each mode the frequency increases as more of the 


tip is removed. This behavior is true for the other two 


Fic. 4 Plate configurations for Series III 

series also and is to be expected, since an increase in 
A///, brings about a decrease in the mass of the plate. 
The fundamental mode of all five plates can be gener 
ally classified as a spanwise bending mode since only a 
little bending occurs in the chordwise direction. Simi 
larly, the second mode can be classified roughly as a 
torsional mode, again remembering that some spanwise 
bending does occur. It is of little use to attempt a 
sunple classification of the remaining modes since they 
are composed of varying amounts of chordwise and 
spanwise bending. 

A clear picture of how the amount cut from the tip 
affects the frequency of the resulting plate is given by 
the graphs of Fig. 8. Here percentage change in fre 
quency, (w — w)/a), is plotted against the parameter 
Al//, for all six modes of all six plates. Each curve 
corresponds to a given shape family so that the per- 
centage increase in frequency is based on that mode 
from which the computed mode evolved. Families A, 
B, E, and F are composed of the Ist, 2nd, 5th, and 6th 
modes, respectively, of all six plates. Family C con- 
sists of the 3rd modes of plates I-0 and I-1 and the 4th 
modes of I-2, I-4, and I-6 while the 4th modes of I-0 
and I-1 along with the 3rd modes of I-2, I-4, and I-6 
make up family D. 

The curves show that for a given increase in A///, there 
is a greater percentage increase in frequency in families 









































TABLE 1 (Continued 
SERIES III 
Family A B C D E F 
Mode 1 2 3 _k 5 rc 
At/Ly f wp f Wp f Wp f wp f wp f wh 
) 26.3 4.50 | 101 17.6 | 171 29.3 | 259 44.7 | 346 59.7 | 522 90.0 
ft 27.9 eT} 110 19.2 | 1a 32.51/274 47631376 64.9] 525 90.5 
Ae 3165 5.39 122 ahe3 198 33.9] 289 49.8} 438 75.6} 542 93.5 
025 34.8 5.96 | 130 22.7 215 36.8} 300 51.7] 476 82.1 | 567 97.8 
03 38.5 6.59 136 23.8 24,3 41.6] 312 53.8} 505 87.1 | 623 107 
035 L4e9 7.69 143 25.0 277 47.4 | 327 56.4] 540 93.2] 674 116 
4 51.7 8.85} 151 26.4 | 314 53.8 | 347 59.91 573 98.9} 699 121 












































624 JOURNAL OF THE AERONAUTICAL SCIENCES—SEPTEMBER, 1954 
TABLE 2 
Plate Properties 
Property Experimentally Determined y/ 
Spanwige E Chordwiseg E fe 3 h (assumed) 
Plate (#/in.“) (#/in.“) (#/in.7) (in. ) 
6 6 
Series - I 29.3 x 10 31.5 x 10 281 ~0622 29 
; 6 6 
Series - II 2942 x ID 2706 x 10 28) 0613 029 
: 6 6 
Series - III 29.0 x 10 2746 x 10 281 0613 29 
10" x 10" 6 6 
29.3 x 10 31.5 x20 282 0627 029 
square 
gu " 6 
18" x 28 29.8 x 10 . 282 0585 29 
beam 




















A and C than in any other family. The fact that Fam- 
ily A is primarily a spanwise bending mode and Family 
C has a large percentage of spanwise bending in it while 
the remaining four families lean heavily toward tor- 
sional modes suggests that removing portions of the 
tip has a greater effect on the frequency of spanwise 
bending modes than it has on chordwise bending, or 
“torsional’’ modes. This can be explained by the fact 
that the spanwise stiffness increases rapidly as the 
effective span is decreased while the torsional stiffness 
is affected only slightly by these changes. 


Series II 


The node lines for the first six modes of plates II-1, 
II-2, II-4, and II-6 along with those of II-0 representing 
the triangular plate from which Series II is developed 
are given in Fig. 6, arranged again in order of increas- 
ing frequency from left to right. The frequencies for 
all modes of all plates in Series II including II-0 are 
given in Table 1. It can be seen from the photographs 
that for this series each given mode makes up a distinct 
A general classification can be extended 


shape family. 
The fundamental and 


to three modes of this series. 
third mode are primarily spanwise 
while the second mode is quite clearly a_ torsional 


bending modes 


mode. 

The graphs of Fig. 9 show how the frequencies of the 
plates of this series are affected by tip removal. The 
increased effect on bending modes is again shown by 
the behavior of the curves of the shape Families A and C 
which correspond to the first and third modes. For 
A//l, in the range 0 to 0.2 Curve E behaves much like 
Aand B. A check of the node lines of the fifth mode in 
this range shows that the mode consists of a large per- 


centage of spanwise bending. From 0.2 to 0.4 Curve E 
behaves much more like B which represents the second 
mode. A study of the node lines of the fifth mode of 
plates II-4 and II-6 bears this out as it shows that a 
considerable amount of torsion is present along with the 


bending. 


Series III 

The node lines of Series III are presented in the 
photographs of Fig. 7 arranged in the same manner 
as the first two series. Table 1 and the graphs of Fig. 
10 show the effect of tip removal on the frequency of 
Again in this series all modes of a given or- 
The same general 


vibration. 
der belong to the same shape family. 
classifications can be applied to the first and second 
modes of this series that have been given the 
responding modes of the preceding series. The other 
modes are too complicated to be given simple classifi- 
This is borne out by the absence in Fig. 10 


cor- 


cations. 
of the distinct difference in behavior between bending 
and torsional and combined modes so apparent in Figs. 
S and 9. However, Curve C in Fig. 10 does seem to 
indicate that the third mode might behave as a bend- 
ing mode even though it may not appear so from its 


node lines. 


FREQUENCIES OF PLATES OF SIMILAR 
PLAN FORM 


The frequencies given in Table | can be used to pre- 
dict the frequency of a plate with a similar plan form 
but having different size, thickness, and material prop- 
erties by using Eq. (1) in which the primed quantities 
refer to the plate whose frequency is to be predicted 
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Node lines for Series I plates 
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Fic. 6 Node lines for Series II plates 
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Fic. 7. Node lines for Series III plates 
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Fic. 8. Effect of tip removal on frequency, Series I. 
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Fic. 9. Effect of tip removal on frequency, Series II. 


and the unprimed quantities refer to the plate whose 


frequency is known.! 


; l ae oo (l > gp) 
f = } l 
a? h LE p’ (1 — pw”) 


a is a scale factor relating all the dimensions of the two 
plates except the thickness. It is shown in the de- 
velopment of this expression in reference 1 that this 
conversion formula is exact only for plates having the 
same Poisson's ratio but it also provides a good ap- 
proximation for plates whose Poisson’s ratios are not 
much different. 

The use of Eq. (1) is made less cumbersome by re- 
arranging it and combining the quantities relating to the 
plate of known frequency into a dimensionless fre- 
quency. Introducing C’ = aC into Eq. (1) gives 


‘ Cr |= p (1 — pn?) f 
f! = ? 
CC" Le 6 GU — 2") 
or 


oH E' St [ ~ i oe 
“C1 p'(1 — pw) ] Uh E 5s Mes 


by E - ad ee @ 
h E 2rV 12 V D/phCt 


where w = 22f, D = [Eh3/12(1 — w*)]. Therefore Eq. 
(3) can be written as 


w’ = wp V D'/p'k'C'* (5 
where 
wp = w/V D/phC* = dimensionless frequency (6 


This form of the frequency is used in reference 3. 

Eq. (5) gives the unknown frequency of the new plate 
in terms of its properties and dimensions and the di- 
mensionless frequency, wp. It is seen that once wp 
is computed, the user of this transformation need no 
longer be concerned with the properties and dimensions 
of the original plate for which the frequency is known. 
This development also eliminates the scale factor a 
and introduces in its place a significant dimension of 
the plate, the root chord. The dimensionless fre- 
quencies tabulated in Table 1 were computed from Eq. 
(6). 


FREQUENCIES OF PLATES OF DIFFERENT PLAN FORM 


The graphs and tables of this report may be used to 
estimate the natural frequencies of a plate having a 
plan form within the range, or slightly outside the range, 
of the plan forms actually tested. This may be done 
by a simple interpolation or extrapolation. Confidence 
may be placed in such a procedure as long as the numer- 
ical values being interpolated between, or extrapolated 
from, are not far apart. A study of the effects of as- 
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EFFECT OF TIP REMOVAL 


TABLE 3 
Dimensionless Frequencies for Series I, 
Based on Span 


Mode wp for Al/lb = O 
l 5.50 
2 14.7 
3 27.5 
4 29.8 
5 46.5 
6 57.0 


pect ratio and sweepback on triangular plates' should 
also prove helpful in determining whether the fre- 
quency behavior is sufficiently linear in the interval be- 
ing considered. The method will be demonstrated by 
several examples. 


I) Interpolation mage to Plate of Truncated 
Delta-Wing Plan Form 
In Fig. 11(a) is shown a delta-wing plan form which is 
intermediate in shape between the plates of Series I 
and Series II. The aspect ratios of the complete tri- 
angular delta-wing shapes are as follows: 
Aspect Ratio = 2L/C = 2(20)/20 = 2 
Aspect Ratio = 2L/C = 2(20)/10 = 4 
Aspect Ratio = 2Z£/C = 2(20) 2.5 


Series I. 
Series IT. 
Plate (a). 


16 = 2.5 


The graphs of frequency vs. aspect ratio for complete 
triangular plates' reveal that in the range from 2 to 4 
the frequency of the first six natural modes are ap- 
proximately linear functions of aspect ratio and it was 
decided to use this as a basis for interpolation in the 
truncated plates as well. 

Furthermore it has been found that interpolation 
produced better results, on the whole, when the semi- 
span of the complete delta plate was used to compute 
dimensionless frequencies, rather than the chord. 
Noting from Eq. (6) that the dimensionless frequency 
varies as the square of the length, C, then the dimension- 
less frequencies of Table 1, Series I should be divided 
by 4 if the 10-in. semispan were used to compute 
wp. Series II dimensionless frequencies would be un- 
affected since the semispan and chord are equal. Table 
3 gives the dimensionless frequencies for the complete 
triangular plate (A///) = 0) of Series I. 

For a numerical example, consider the second mode 
of vibration (Family B) of the plate shape in Fig. 11(a). 
Al /) = 3.8/10 = 0.38. From the graph of Family B, 
Series I, Fig. 8, the frequency ratio is Aw/w = 0.22. 
From the graph of Family B, Series II, Fig. 9, the fre- 


quency ratio is Aw/w) = 0.27. Then wp, for Series I 
is 1.22 & 58.7/4 = 17.9 and wp for Series II is 1.27 X 
TABLE 4 
Comparison of Interpolation with Experimentally Observed 


Frequency, Plate (a 


Mode Family — A B D te E F 
Mode Number—> l 2 3 4 5 6 
Predicted frequency 55 118 241 308 372 465 
Measured frequency 52.2 116 225 308 370 447 
Per cent difference +5.4 4+1.7 +7.1 0 +0.5 +4.0 
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Fic. 10. Effect of tip removal on frequency, Series ITI 


23.8 = 30.2. Interpolating, we obtain w, for Plate (a) 
= 17.9 + (0.5/2.0) (12.3) = 21.0. 

Next the frequency of an actual steel plate is pre- 
dicted, assuming the elastic properties as given below. 


h = ().059 in. 

E = 29 X 10° Ibs. per sq.in. 

rm = 0.29 

D = Eh*/12(1 — yw?) = 560 Ib. in. 

(//2)* = 10* in.‘ 

p = (0.282/386) = 7.31 & 10~* Ib. sec.*/in.4 
Then 

o> Wp \ D phil 2 ‘= 39.4w p (7) 

Thus for the second mode of this plate, f = [(35.4 X 
21.0)/2r]| = 118 cycles per sec. The experimentally 


measured frequency of this plate was 116 cycles per sec. 
By this procedure of interpolating between Series I 
and II to find w,, and then using Eq. (7), the first six 
natural frequencies of this plate were calculated and are 
compared with experimental results in Table 4. 


II) Interpolation Applied to Plate of Truncated, 

Sweptback Plan Form 

In Fig. 11(b) is shown a plate which is intermediate 
in shape between the plates of Series II and Series ITI. 
The completed triangular wing has a total span of 20 
in. and root chord of 10 in. as do the plates of Series IT 
and III, so that all three plates have the same aspect 
ratio. The graphs of frequency vs. tangent of sweep- 
back angle for complete triangular plates' reveal that 
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in the range from 0.5 to 1.0 the frequencies of the first 
six natural modes are approximately linear functions of 
the tangent of the sweepback angle (based on mean 
chord line) and it was decided to use this as a basis for 
interpolation in the truncated plates as well. 

As an example, consider the third mode of Plate (b) 


in Fig. 11. 
Al/Ilg = (3.25/10) = 0.325 


From the graph of Family C, Series II, Fig. 9, the fre- 
quency ratio is Aww = 0.49. From the graph of 
Family C, Series III, Fig. 10, the frequency ratio is 
Aw/we = 0.51. 

Then w, for Series ITI is 1.49 K 32.4 = 
for Series III is 1.51 XK 29.3 = 44.2. Furthermore the 


sweepback tangents are 


18.2 and wp 


tan a = 0.5 
tan a = 1.0 
tan a = 0.7 


Series I. 
Series II. 
Plate (b). 


Then w, for Plate (b) = 48.2 — (0.2/0.5) K 4.0 = 46.6. 

Next the frequency of an actual steel plate is pre- 
dicted, assuming the elastic properties as given below. 
The modulus of elasticity is based on experiments 
made on plates of similar steel. 


h = 0.0613 in. 

E = 28.5 X 10° Ibs. per sq.in. 

pw == 0.29 

D = [Eh'/12(1 — w?)] = 598 Ib. in. 

C* = 10! in. 

p = (0.281/386) = 7.28 K 10~ Ib. sec.?/in.4 
Then 

Qo = Wp D phC# = 36.68 p (S) 

Thus for the third mode of this plate, f = [(36.6 X 
46.6)/2mr] = 272 cycles per sec. The experimentally 


measured frequency for this plate was 270 cycles per 
sec. A comparison between predicted and measured 
frequencies for the first six modes of this plate is given 


in Table 5. 


SCIENCES SEPTEMBER, 1954 

As a closing remark on interpolation, it seems evident 
from the photographs of Figs. 5, 6, and 7, that interpola- 
tion may be applied to shapes of natural modes as well 
as to the frequencies of vibration. 


THE EFFECT OF NONISOTROPY IN PLATE MATERIAL 


It is well known that the rolling process introduces a 
directional character into the grain structure of metals, 
the grains being longer in the direction of rolling. 
This orientation of grains is only partly relieved by 
annealing. It has been found that the ultimate ten- 
sile strength and ductility of steel may vary consider- 
ably, depending on the direction in which these prop- 
erties are measured in relation to the rolled plate.’ 
In attempting to check the accuracy of our techniques 
by comparison with work of Barton® and Dalley and 
Ripperger' we began to suspect a difference in the 
moduli of elasticity in the spanwise and chordwise di- 
rections for the steel plate being used. Careful static 
beam deflection tests revealed that the moduli in the 
two directions differed by as much as 7 per cent, one 
plate having 29.5 X 10° Ibs. per sq.in in the direction 
of rolling and 31.5 X 10° Ibs. per sq.in. perpendicular 
to the direction of rolling. 

In determining the modulus of elasticity, two beam 
specimens, '/, and '/: in. wide, were cut from both the 
spanwise and chordwise directions for each of the 
three original triangular plates. The spanwise direc- 
tion in Series I and the chordwise direction in Series II 
and III coincided with the direction of rolling. The 
material was hot rolled steel approximately 1/j¢ in. 
thick. 
to give a constant uniform cross section throughout. 


The edges of the beams were accurately ground 


They were then placed on specially constructed knife 
edges and deflections were measured at the center of a 
10-in. span, under the action of centrally applied loads. 
The load on each beam was gradually increased to a 
point near the limit of elastic behavior, and then 
gradually decreased to zero. Curves of load versus de- 
flection were then drawn to insure that the tests were 
conducted in the range in which deflection and load 
were linearly related. The modulus of elasticity was 
then computed from the compiled deflection and load 
data, using elementary beam theory. No significant 
difference was found between the moduli given by the 
'’, and '/,-in. beams and it may be concluded that 
specimens of both widths were amenable to simple 
beam theory. 


TABLE 5 
Interpolation with Experimentally 
Frequency, Plate (b) 


Comparison of Observed 


Mode Family—> A B ee D E F 
Mode Number—> l 2 3 } 5 6 
Predicted fre- 

quency 48.2 158 272 366 581 668 
Measured fre- 

quency 47.3 150 270 350 579 650 
Per cent dif- 

ference +1.9 +5.3 +0.7 +4.6 +0.3 42.8 
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Table 2 gives the results of these experiments in the 
form of an average spanwise and chordwise modulus of 
elasticity for each plate. Also included in this table 
are the measured densities of the plate material and the 
values of Poisson’s ratio used in computing the dimen- 
sionless frequencies. 

The variation of modulus of elasticity with direction 
in nonisotropic plates poses the question of what value 
of & should be used in predicting new frequencies and 
computing dimensionless frequencies to insure the 
greatest accuracy. Since bending in a given direction is 
related to the modulus of elasticity in that same direc- 
tion it appears that it would be proper when attempt- 
ing to predict, say, the frequency of a spanwise bend- 
ing mode to use the spanwise modulus of elasticity. 
The dimensionless frequency used for this purpose 
should be given the same consideration if it is obtained 
from a plate of nonisotropic material. This procedure 
is quite straightforward for modes consisting primarily 
of either spanwise or chordwise bending providing the 
moduli of elasticity of the material are known in 
these directions. However, modes in which bending in 
both directions are combined make it necessary to use 
some average of the two moduli depending on the mode 
shape. The spanwise /’s were used in computing the 
dimensionless frequencies of shape Families A and C in 
each series while the chordwise -’s were used for shape 
Families B. Due to the complicated shapes contained 
in Families D, E, and F, arithmetic averages of the 
spanwise and chordwise moduli were used in computing 
their dimensionless frequencies. 

Another effect of nonisotropy should be pointed out. 
The shape of a natural mode is unquestionably altered 
irom the shape that would exist if the plate were iso- 
tropic. We feel that even for appreciably different 
stiffness in the spanwise and chordwise directions, the 
effect on modes which are predominantly spanwise or 


chordwise bending would be quite small. For more 
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complicated modes the effect would be more 
nounced, and no doubt some of the node lines in Figs. 
5-7 are not quite the same as would be observed for a 


However, it is unlikely that the 


pro 


truly isotropic plate. 
7 per cent maximum difference in spanwise and chord- 
wise moduli has introduced serious differences in shape. 


CHECK ON ACCURACY OF EXPERIMENTS 


In order to establish some idea as to the degree of ac 
curacy of the experimentally determined plate fre 
quencies, and as a check on effectiveness of the clamp- 
ing, the first five natural frequencies of a 10- by 10-in. 
cantilevered square plate and of an 18- by 2!/»-in. beam 
were determined using the same test equipment and 
procedure used on the trapezoidal plates. A compari- 
son was then made between the observed frequencies 
The 


results of this investigation are presented in Table 6. 


and existing experimental and theoretical results. 


Before the dimensionless frequencies were computed, 
tests were run to determine the moduli of elasticity of 
the square plate and beam, which were cut from ' /;¢-in. 


commercial cold-rolled steel. These tests were con- 
ducted in a manner similar to those previously dis 
cussed with the exception that only the spanwise 
modulus was determined for the beam specimen since 
only spanwise bending occurred in the first five modes. 
The moduli determined for the square plate and beam 
are listed with those of the triangular plates in Table 2. 
The first and third modes of the square plate are 
primarily bending modes, the second and fourth con 
sist chiefly of torsion, and the fifth is a combination of 
both. Hence in computing the theoretical frequencies 
the spanwise modulus was used for the bending modes, 
the chordwise modulus for the torsional modes and an 
arithmetic average of both was used for the fifth mode. 

When the appropriate modulus is used, the dimen- 
sionless frequencies obtained for the square plate agree 


Continued on page 648 


























TABLE 6 
Comparison of Analytical and Experimental Frequencies for Square Plate and Beam 
. 
10" x 10" Plate 16" x 23" Beam 
Ww Author's | Computed Experimental Author's Comouted 
0 (experimental)! value 4 value (experimental) value % 
Mode value reference (3)i deviation | reference (4) value reference (§) devia» 
1 3.35 3.49 - 4.0 3237 3.46 3.52 ~}.9. 
2 8.53 8.55 — 0.2 8.26 21.54 22.4 -4.0 
3 20.90 21.44 = 2.5 20.55 60.85 61.70 - 1.4 
4 26.72 27.46 = 2.7 27.15 119.9 121.0 - 0.9 
5 30.61 31.17 —- 1.8 29.75 197.5 200.0 - 1.3 
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A Note on the Laminar Boundary Layer on a 
Circular Cylinder in Axial Incompressible 


Flow 


Howard R. Kelly 
U.S. Naval Ordnance Test Station, Inyokern, China Lake, Calif. 


February 10, 1954 


VERY INGENIOUS AND useful method of computing laminar 
A boundary layers on circular cylinders in incompressible 
flow has been published by Seban and Bond.! Among their 
published results they compute the ratio of the displacement 
thickness 6, on a cylinder to the corresponding thickness on a 
flat plate to be 


5:/dp. = 1 — 0.029F + 0.0548? . .. (1) 


to second order in the transformed axial coordinate 


— = (4/a)V vx/uo (2 


In a later paper, Eckert? mentioned that Seban and Bond ob- 
tained ‘‘no clear results because of the limited number of terms 
in the series expansion.”’ It is felt by the present writer that the 
series in £ should converge rapidly for —< 1.0, but that the re- 
sults in Eq. (1) are in error for another reason, as will be shown 

The authors of reference 1 define a displacement thickness by 


. 
i= | (1 — u/uo) (r/a) dr ( 
a 


This definition accounts for their use of cylindrical coordinates 
and for finite thickness of boundary layer in the right-hand term, 
but does not account for them in the left-hand term. In order to 
account for the mass-flow defect in a way comparable to flat-plate 
theory, one should obtain Eq. (3) by the relation 


a + 6* y 
uj(r/a)dr = (to — u)(r/a) dr (4) 
a a 


which leads to 


the equation 


ww 


Px 
6, = 6*(1 + 6* 2a) = | (1 — u/uo) (r/a) dr (5 
a 


where 6* is the true displacement thickness. Solving for 6* in 
terms of 6; from Eq. (5) 
6* = 6,(1 — 6,/2a + 6,?/2a?... ) (6 


The solution of Eq. (3) given in reference 1 may be written 


5, = 1.721 VW wx/uo (1 — Efi( )/1.721 — £fx( ©) /1.721 .. . ) (7) 
where fi(o) = 0.051, foo) = —0.093. Substituting this in 


Eq. (6), then using Eq. (2) and retaining terms only to £2, we ob- 


tain 
6* = 1.721 V wx/uo (1 — 0.245¢ + 0.160? . . . ) (8) 
so that, if 6*r.p, = 1.721 VY vx/to, 


6*/d*rp. = 1 — 0.245 + O.160E? .. 9 


Comparing Eq. (9) to Eq. (1), it is seen that the results given by 
Seban and Bond are in error by a factor of about 8 for the linear 
term in é and a factor of 3 for the quadratic term. As compared 
with Eq. (1), Eq. (9) shows results which are more nearly the same 
order of magnitude as those of Eckert for the turbulent boundary 
layer. 

Some further corrections are also indicated for reference | 
Due to a typographical error, Eq. (26) contains a term 3f,y’ 
The correct equations have been 


which should have been 3f "fs. 
integrated at the Naval Ordnance Test Station+ in connection 
with another problem, using the REAC, an analog computer that 
allows continuous integration with a high degree of accuracy. 
The Blasius function fy was duplicated quite accurately, while the 
solutions for f; and fz varied somewhat from those given by Seban 
and Bond, as follows: 


Seban and Bond NOTS REAC 


f,"(O) 0.704 0.696 
fe”(0) —0.095 —(). 160 
fil 2) 0.051 —0.019 
fol ~ ) —(). 093 —(). 100 


It is believed that the difference is probably due to the finite in 
terval size used in the former integration 
The above results will again change the equation for displace- 


ment thickness. Eq. (8) will become 


6* = 1.721 V wx/uo (1 — 0.204 + 0.14682... (10 
and Eq. (9) will be 
6*/6*r.p. = 1 — O.204E + 0.1468". . (11 


The equation given by Seban and Bond for skin friction 


cr = 0.664 / » xu 1 + 0.58—& — 0.071E (12) 


will become 


~9- 


cr = 0.664 VY v/xuo(1 + 0.525 — 0.120€?) (13 


€ 
so that 


, = wf 595: — Wye2 ( 
C/Ctpp = 1 0.525€ 0.120¢ (14 
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t The REAC integrations were performed by R. H. Olds of the 
Mathematics Division. 


634 





rh 


whi 


Col 


whi 


Th 


Th 


Li 


He 


een 
‘ion 
hat 


the 


yan 





READERS’ FORUM 635 


The Kernel Function of the Integral Equation It is convenient to express this function as a power series of 
Relating the Lift and Downwash the circular frequency w. Such a series is obtained in reference 1, 
Distributions of Oscillating Finite Wings in = Eq. (26), and the object of the present note is to show how the 
Subsonic Flow* result may be obtained in a more straightforward man 

ner 
N. Lancet We define the Laplace transform of f(£) to be F(p) where 
Department of Engineering, University of California, Los Angeles E 
March 31, 1954 F(t { f(t )e— PE di ° 
. 0 


I 4 RECENT PAPER Watkins, Runyan, and Woolston! have dis- Consider the function 
cussed the kernel function 


. Ko pt lA WA + B i g 
' li o* xo/ J f 1 
A(xo, Yo) = lin é s t) = d , 
- ; . ‘ Ar +r oh tT OH 


2 —> 0 02? 5 . 


Vg M(x? + B?yo? + Bs The Laplace transform of Eq. (3) with respect to & is 
dd l 
A? + By rT §*2* . 
dd 
Fip) = } 
* This work was carried out under the auspices of F.O.A., Scientific Re \2 + p n+ Mir? + ¢ 
search Project TA 01-101-3006 (OEEC 151 * 
+ The author wishes to express his gratitude to Prof. John W. Miles who 
gested the use of the Laplace Transform Method discussed herein where 7? = 8*(y)? 4 


he integral for F(p) can be evaluated in terms of elementary functions and we obtain 


" i > 9\!/s 
= - e p* p*r* 7 
F(p) = tanh™! j ' ' “ 
(p? + B*r 2 Pas M + 1)j (x? + 7° ~ *oy 


In order to obtain a series of K(x), yo) valid for small & we expand F(p) in powers of p~', thus 
B?r? l { 2 2 ly : : 2 ' 
F(p) = — log [((M + 1)j) (x0? + r?)/? — xo /2p] 4 (M+ 1) \(xo? + 7? - Xo} + 
2p> pp 2p? 


B?r? l l Birt 
; + B?r? (MU + 1)2) (xo? + 72)? — xof? + + (Kp 
WU + 1)} (xo? + r?2)/? — xo} 4p 2 20.7 + 1)?) (xo? + 7? Xo} 2 y 


his expression can be differentiated twice with respect to s and 2 equated to zero subsequently, with the result that 


0? F(p) B , Xo + Ss g? B2[xo + (1 — WM — M?)s 
= log [((.M + 1) (s — xe)/2p] - a + O(p-4 6 
Os? z= 0 p* PSyo" prs p's 
where s = (x)? + B*yo?) /? 
We take the inverse transform of Eq. (6), which involves only simple functions of p and obtain 
O*f(E) Bie? 1 3 ; B7E&(s — Xo) (xo + 5 BrE B7E2[x%o + (1 — M — M?)s : m 
= — 7 — log = — of ih + OE r 
Os? z 0 2 2 21 — M) So? § 2s . 
Collecting together like powers of w and multiplying by e~ ‘*** VY) gives 


Ae er re, J _ Xo + Ss we w? w(s — X z w? a ] im l _* w*\| 
K(xo, Yo) = € + log Y + M + 0 8) 
: ( S¥o? Vs 2V? 2vVal — M 2V2 2 2 8? 5 B°} | 


which is the same as Eq. (26) of reference 1 


The series for K(xo, yo) when .f — 1 can be obtained similarly. We write 17 = 1 — e€ and find 
/¥ : lim 0? x¢ elf 2e} [A — (1 eid? + ely + 2? 
e\texX0/") K(x, Yo)| M=1 = s—>0 1d 3 
: : ds? 1A? + 2e(yo? + 2%)} : 
«>0 ? eet ad 
Transforming again with respect to £, which is jw/ V in this case, we have 
$ ” 2e dd 
F( p) = ; ; ; ; ; (10) 
_o [Zep — A H+ (1 — EA? + Ze yo? + 27) AA? + Ze o*® + 27)} 
The integration with respect to \ can again be performed and if «— 0 we obtain 
: : 1 (et + a 4 
Lim F(p) = —- Z - tanh” |; “ . ; 11) 
«0 ie? -— 3 + Ss") "" iP + |(y¥o* + 3*)/Xo} 5 


The series expansion of Eq. (11) is 


: : ie a 1 y° + 3% 1 yo? + 2? 1 1 (yo? + 22)? 
Lim F(p) = _ log + xo + + 3 Xo? + Yo* + 2? — — + O(p-*) (12 
«> 0 4p3 2p 2pxo 4p? Xo 8p \2 2x0? 


Hence 
. OF(p) Te ] 1 1 2xop 1 1 Yo? 
a = — _— - log + + -3- + O(p~é 
«0 2? 2=0 2p* p/ yo? 2p? Yo? 2xop? 4p3 Xo? 
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Thus 
K saxe/V) 3 2 lw w? , wy" w* iw V9? ; ) 
(Xo Y)|M=1=e oo ae r, 108 > = —ae es - + O(w*) > 1 
ty? Vxo 2V? 2Vx 2V? 2 2x9? { 
iresult that could have been obtained directly from Eq. (8) provided due consideration is given to the order of terms neglected 2 
The value of K(xo, yo) may also be obtained for the case when w —> Copson? shows that if -' 
*) 2 
le J ere” WE) a 14) 2 
a 
then, provided @'(£) # 0 anywhere within the range a < & < 3b, 
y( h jet? b "7 a etvela 
[= + O(w~? 15 
lwo | b Twp (dad) 
The formula, Eq. (15), may be applied to Eq. (1) and we derive the expression 
( [4 My. ‘7 [ x i } 
; _ OF Jyge Lb Vo Ve . Vee LV Vat” Vas diet 
K(xo, Yo) = lim . ~ lim + O(w-2 
s—> 0 Os" | iw (xy? + r?) — Mxo] \— tw[(A? + r2)/7 + Nr 4 
M 2x0 \ 
; ( ial 1x 1 ; ‘| / Eq 
° . < | 3°¢ Y VB > 
= lim : é' oe a 
2—> 00s" | tw |(X9? + 1?) - Mx \ re] 
Hence bv 
w 1 
G )( Mix — 2 ) a) 
, MBre het 
K(xo, Yo) = — - + O(a ~! (16 
SS — Mxy su 
which is the value of K(xo, yo) valid for large w. The term displayed in Eq. (16) is sufficient for practical applications but it should be the 
noted that later terms will contain (s — \/xy)~’ where r > 1 _ 
REFERENCES ( 
1 Watkins, C. E., Runyan, H. L., and Woolston, D. S., On the Kernel Func- 2 Copson, E. T., The Asymptotic Expansion of a Function Defined by a Def 
tion of the Integral Equation Relating the Lift and Downwash Distributions of inite Integral or Contour Integral, Admiralty Computing Service, S.R.F. in} 
Oscillating Finite Wings in Subsonic Flow, NACA, TN 3131, January, 1954 ACS 106 , 
nel 
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? ) 
the 
A Note Concerning the Effect of Wing Twist _ by the relation for 
on the Response of a Wing to a Sharp-Edged eit = he) —ee (1 
Gust ww 
where yw is the mass parameter defined in conventional symbols 
: by 
Franklin W. Diederich Re 
am 
California Institute of Technology, Pasadena, Calif. (on leave ia 8. - , 
from Langley Aeronautical Laboratory, NACA) SC 2 
April 27, 1954 ; ' : 
and is equal to 4 times the mass parameter A used in reference 1; 
‘tt | ; _ the quantity gp is the dynamic pressure at divergence, which for ‘yt 
: SCENT PAPE > acceleration response to a sharp-edge } , ‘ gps ; 
T° A RECENT mae ilk i ne aoe gee Se ss wy rs Bee sweptback Wings Is negative (hence fictitious), so that —qdp isa ‘ 
an < ane Ww one degree o -edom, vertical mo- ee ; a ; : , 
gust of an airp ‘an ith one en ; sete ; read cae me positive quantity. The stiffness parameter cR thus involves the an 
i ras calculated us insteady-lift functions for incompres- ; : Seite 
Gon, was calcul ex age nettles : ; ; sd A apes nei mass parameter as well as a true index of the wing flexibility 
. y > basis he assumption that win sxibility car on ; ‘ é : 
sible flow - the basis . one san jon that 8g oe ; ‘ si The effects of mass (or inertia) and stiffness can be separated by 
> take ace -onsidering an average twist which is . ; pers aes 
be t iken into ss mS nyse “ 8 E ot Be si . ‘ : plotting the results of reference 1 in the form used in Fig. 1 The 
é > al acceleration. is mption i . . Ae 
sete ; ss a sence eres i : ; be “ “ ce wing-twist alleviation factor y can then be seen to be substan 
~é é P 2 t 1buted by the tall and fuselage is not : . . . 
reasonable if the : ? — 7 ' A i! sce i . tially independent of the mass parameter and to depend primarily Tt 
arge > mi > wi including nacelles, tip tanks, etc.) . hag ue 
large, the mass 0 a ae ( rpamealiig _— me re 2 van on the stiffness parameter (q¢/—qp); furthermore, it is given | 
all, < » effective : -onsta ; selected properly ‘ ; 
is small, and the e ecti re spring const nt is selecte i roperly approximately by the relation 
Also, the neglect of pitching motion while including twist can be , 
justified if the wing is relatively flexible, the moment of inertia y = e4 
of the airplane in pitch is large, and the airplane has little or no 1 + (g/—qp) ne 
sweep. However, if the wing is highly swept back, the practical c. ss : (3) ( 
Nae : , ar 
significance of the results obtained by assuming that a sharp- = C Ma 
edged gust hits all points on the leading edge simultaneously is La 
somewhat uncertain. where Cy, is the lift-curve slope corrected for static aeroelastic 
Without attempting to delimit the conditions for which these rs F 4 me i eff ; i f ; | 
. aa . ; effects.*, Some of the dynamic effects ignored in reference 1 c: 
results are valid and significant, the purpose of this note is to re- as , lai tai - lati va “ = (3 re = — . “a sac 
; A oe ye included in the simple relation given by Eq. (: asi-steady ur 
cast the results of reference 1 in a form which exhibits the rela- ; i eff " re i . given by Eq. (9) ae rehes es 
é . dynamic effects (deflections due ‘onstant-g pull- are in- 
tive effects of the mass parameter and the stiffness parameter y dint agia \cemections due to constant-¢ pull-out) are in 
cluded in Cr ek 
somewhat more clearly. ay : 
The stiffness parameter cR used in reference 1 can, by reasoning The conclusion that the alleviation factor is relatively insensi- * 
similar to that of Appendix C of reference 1, be shown to be given tive to mass parameter and can be represented approximately by 
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Eq. (8) had been reached previously as a result of unpublished 
calculations performed for positive effective spring constants, 
representative of unswept and swept-forward wings (for which the 
alleviation”’ is actually an overstress) by the undersigned and 
by William H. Phillips (who included pitching motion) and thus 
ppears to be valid quite generally, within the assumptions in 
herent in this type of analysis. The maximum acceleration pur- 
suant to entry into a sharp-edged gust can thus be obtained within 
these approximations from a rigid-wing analysis, such as the 
one of reference 3, but using the lift-curve slope C7, which is cor- 


rected for static aeroelastic effects instead of the lift-curve slope 


a in the mass parameter yu 
Che effects of ‘‘twist’’ considered here (which is actually ar 


.¢ 


ingle-of-attack change due to structural deformation and is, 
hence, caused largely by wing bending in the case of a sweptback 
wing) represent only a part of the effects of wing flexibility on the 
acceleration response to a sharp-edged gust. However, used in 
conjunction with a simple analysis of the dynamic effects due to 
vertical deflection of the wing (such as the one of reference 4) 
they serve as a quick method of estimating gust loads in the cases 


for which the assumptions made are valid 
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wept Wings, NACA TN 2608, February, 1952 
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irbitrary Shape, NACA TN 2036, February, 1950 
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ng Flexibility, NACA TN 2763, August, 1952 


Airplane Including Wing Bendi 


Theoretical and Experimental Investigations 
on Swept and Delta Wings in Incompressible 
Flow* 


E. Truckenbrodt 

nstitute of Fluid Mechanics, Technical University of Braunschweig, 
Germany 

March 29, 1954 


r I ‘HIS IS A BRIEF SUMMARY of some systematic investigations on 
swept and delta wings in incompressible flow, which will be 
published in detail later.! 


* These investigations have been sponsored by the Air Research and De 


elopment Command, European Office at Brussels, Contract No. AF 61(514 
426) \ brief report was given by H. Schlichting at the IAS Meeting in 


New York on Jan., 26, 1954 
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Fic. 1. Plan forms of swept and delta wings (dimensions in 


millimeters), profile NACA 0012 


Qur investigations are an extension of those by W Jacobs,? 
who did systematic six-component measurements and pressure 
distribution measurements on swept wings with different angles 


ot sweep but constant aspect ratio, 1 = Dia Systematic six- 
component measurements on delta wings have been done earlier 
by Lange and Wacke 

We wanted to explore the influence of the aspect ratio on the 


verodynamic parameters of swept and delta wings in symmetrical 
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Fic. 3. Lift slope and aerodynamic center of delta wings. 
XG.c./C¢ = 0.493. 
flow. We carried out three-component measurements on five 


swept wings and five delta wings and also theoretical calculations. 
In Fig. 1 the plan forms of the ten wings are given. The five 
swept wings have constant chord and constant angle of sweep, 
¢ = 45°, with aspect ratios between A = 5.2andA = 1.2. The 
five delta wings all have a taper ratio of 1:8 and aspect ratios 
between A = 3.9and A = 0.8. All wings have the profile NACA 
0012 and have no twist. The measurements were done at a 
velocity of 40 m.p.s., thus the Reynolds Number, based on the 
mean chord, is 6.10%. 

The main results for the swept wings are presented in Fig. 2. 
The upper diagram gives the slope of the lift coefficient, and the 
lower diagram the distance between the aerodynamic center 
(A.C.) and the geometric center (G.C.) against the aspect ratio 
In Fig. 3 the same results are plotted for the delta wings. In the 


latter case the distance x, ,, is referred to the geometric mean 


+s 
Cm = (1 of cy) dy 
—§ 


Here S is the wing area and c(y) the local wing chord 
These experiments were supported by systematic theoretical 


chord, defined by 


investigations on the spanwise and chordwise lift distributions of 
this series of wings. Since the method of J. Weissinger‘ for 
calculating the spanwise lift distribution of swept wings is not 
sufficient for obtaining the aerodynamic center, we worked out an 


ICAL 
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“Extended Lifting Line Theory’’® and also a full “Lifting Surf 
Theory,’’® which yields a continuous chordwise and spanwise lift 
distribution. Our Lifting Surface Theory is applicable to 
arbitrary wing plan form with arbitrary twist distribution and 
arbitrary camber Multhop; 
According to our Lifting Surface Theory we have calculated the 
spanwise and chordwise lift distribution of all wings of Fig 
The main results of these theoretical calculations are included in 


It is similar to the method of H 


The agreement 
Ther 


Figs. 2 and 3 for comparison with experiments 
between theory and experiment is very good in all cases 
fore the new theory has full credit 
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Remarks on Compressible Flow Actuator Disc 
Theory 


Lawrence Talbot and E. V. Laitone 

Assistant Professor and Associate Professor, Respectively, University 
of California at Berkeley 

April 3, 1954 

ie A RECENT CONTRIBUTION to the Readers’ Forum, v. Krzywo- 
blocki! raised the objection that although the authors of 

reference 2 criticized several previous assumptions*® * made in 

treating the actuator disc in compressible flow, they did not them- 

selves present anything more than the exact incompressible flow 

relations. 

The discussion in reference 2 was intended to emphasize two 
main points: 

(1) Entropy increase through an actuator disc, which pre- 
viously has been ascribed to blade friction losses, can occur in a 
frictionless propeller and is in fact inherent in the actuator disc 
concept. 

(2) The first compressibility corrections to the axial interfer 
ence factor due to the trailing helical vortices and to the rotation 
of the slipstream are of the same order of magnitude as the first 
compressibility correction term associated with the axial motion 
These vortex and rotation corrections are of opposite sign and 
possibly may cancel each other to a first approximation. 

A solution for the one-dimensional actuator disc flow which is 
applicable to situations where the compressibility effects are 
important was presented in reference 5. In this note the as- 
sumption was made that the power input per unit area at the disc, 
P, could be expressed as the product of the disc loading (thrust 
per unit area, 7, and an average velocity through the disc. For 
this average velocity, the simple arithmetic average (V, + V2)/2 
of the two velocities immediately ahead of and behind the disc 


was chosen, so that P = 7(V; + V>2)/2, and it was found that the 
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Fic. 1. Calculations for Figs. 1-3 are based on flight Mach 
Number J, = 0.7 at an altitude of 40,000 ft. The horsepower 
per square foot of disc area, P/550, is defined for the different 
theories according to the relations P = [7(Vi + V2)| /2, Stewart; 





P =§7V;, incompressible; P = !/2 Vi (V3? — Vo?), Vogeley 


Correction: V2 above 800 line should read V3.) 
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states at 1 and 2 were related through the Rankine-Hugoniot 


equations for a plane shock wave. An explicit expression was 
obtained for S; — S,, the entropy gain through the disc, and 
thus for the density ratio po/p;, where the subscripts 0 and 3 refer 
respectively to conditions far upstream and downstream of the 
actuator disc. Since S. — S; > 0, p/p 1 for all finite disc 
loadings 

In Figs. 1, 2, and 3, results calculated by Stewart® using the 
theory of reference 5 are compared with the incompressible flow 
theory and with the results obtained by Vogeley, who assumed 
isentropic compression through the disc. It is seen that Stew- 
art's and Vogeley’s values coincide almost exactly for light disc 
loadings. (Some errors may have been made in reproducing 
Vogeley’s values, since they were read from rather small curves 
Also, it can be seen that the density ratio, pp/p;, does not differ 
appreciably from unity until moderately high disc loadings are 
reached. On the other hand, the incompressible flow solution 
has the wrong slope right from the start and is considerably in 
error even for light disc loadings 

In light of these results, it seems probable that any com- 
pressibility correction iteration series which has for its first term 
the incompressible flow actuator disc solution will be wrong unless 
MJ, is very small. This applies both to the approximation 6 
(z 4q 
V;, and q is the free stream dynamic pressure; and to the approxi 


made in reference 2, where 6; is defined by Vo(1 + 6;) = 


mation p; p2 used by v. Krzywoblocki 
The assumption of isentropic compression through the disc is 
However, calcula 


correct for light disc loadings, and .1/,; < 1 
tions based on (p2/pi)” = pe/p; are not much (if any) easier to 
perform than those based on the theory of reference 5, and the 
latter has the advantage of indicating when entropy gain through 
the disc becomes important. In this connection it is worth noting 
that for 1/; = 1 there is a limiting value of the power for which 
isentropic flow is possible.’ If ? is greater than this limiting 
value, losses must be introduced in the dise in order to obtain 


solutions. 
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On Jones’s Criterion for Thin Wings of 
Minimum Drag 


R. Sedney 
Ballistic Research Laboratories, Aberdeen Proving Ground, Md 


April 4, 1954 


T° REFERENCE 1 Jones has derived necessary conditions for 
minimum drag shapes of thin wings subject to conditions of 


given lift, given maximum drag, or given volume. In deriving 
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these conditions, the concept of the ‘‘combined flow field’’ was 


used. This consists of superimposing the flow fields of the for- 
ward and reversed motions. In this note, it is shown how these 
conditions may be derived using the standard methods of calculus 


theorem.?. The 


of variations and the general reverse flow 
method will be illustrated for the case of given lift. 
it is shown that the necessary condition actually yields a mini- 
mum for the drag. Linearized theory is assumed throughout. 
Let p and a@ denote the local lift and angie of attack distribu 


These are defined on the projection of 


In addition 


tions in forward flow 
the wing surface on a mean plane. This area is denoted by S 
As in reference 3, two 


In the first, the lift in re 


and all integrations are over this area 
types of reverse flow are considered. 

verse flow p = p, and @ is determined from linearized theory 
In the second, the angle of attack in reverse flow a = —a, and 
p is determined from linearized theory. The reverse flow theo 


rem states that 


S pa dS = S pads «i) 


and similarly for the second type. Since p and a@ are related 
linearly, p + eg must correspond to a + €8 where @ is the angle of 
attack distribution corresponding to the lift g. Defining the vari 
ation of p, 6p = eq, then the variation of a is da = «8. The re- 
verse flow theorem can, of course, be applied to 6p and éa. 
Consider now the variational problem of minimizing the drag, 


D, subject to condition of a given lift, L. Since 


D = S pe ds (2 
12 Sp ds (3 


the quantity D + AL, where A is an undetermined constant multi- 


tl 


plier, should be minimized according to the rules of the calculus of 
variations. Then 
6(D + AL) = S (pic + adp + dbp) dS (4 
Since Pp = p, the first term of the integrand can be written 
S Pia dS = S pec ds 
= S ap dS 


Thus Eq. (4) can be written 
(D+ AL) = f(a+ a+ dsp dS 
Since the variation of D + AL must be zero for arbitrary 6p, the 
necessary condition for minimum drag with given lift isa + @ = 
\ = Constant, which is Jones’s criterion. Multiplying the last 
equation by p and integrating determines A, so that 
a+a@=2D/L (5 
The criterion can be expressed in a different form involving the 


second type of reverse flow. Let the flat plate of unit angle of 


attack be denoted by ay,i.e., a; = 1. The lift can be written 
L= f[pdS = fa;pds 
= — fa;pds 
S pia ds 


Taking the first variation as in Eq. (4) 


yields 
6D +rL) = f'(p — pb — py ba dS 


and since 6a is arbitrary p — p = Ap; where now \ = —(2D/L). 
The question of what type of functions p and a are admissible 
in the variational process needs to be discussed. It is known 
that the reverse flow theorem, Eq. (1), does not hold if the lift has 
a leading edge singularity in forward flow, since the wing in re- 
verse flow does not satisfy the Kutta condition at the trailing 
edge. Therefore admissible functions p must be such that they 
do not have leading edge singularities. However there are indi- 
cations that Jones’s criterion is valid even for functions p which 


have these singularities. * 
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Using the concept of orthogonal loadings* 4 it can be shown 
that Jones's criterion is also sufficient, that is, it actually yields q 
minimum. Let a; and a be two angle-of-attack distributions 


and p; and p. the corresponding lifts. Define 


> > 
= 2 + rt 
l 2f (ap a2p 


is one-half the mutual interference drag of the tw 


ds 6 


(@1, Qe 


e., (ai, @ 


Then q is orthogonal to ae if (a;, a2) = 0. Sine 


. 
(ai, @ = JS pia dS 


= D, > 0 


distributions 


it is easy to show that Schwarz’s inequality holds 
(a1, a2)? < (an, a) (ae, a ‘ 
Now let ay be any @ such that 
. 
= == ) 
L JS pi ds ( 


Then multiplying Eq. (5) by po, integrating, and applying Eq. (1 
= 0 so that any a with zero total lift is orthogonal 
The con 


yields (a, ap) 
to the optimal a, i.e., the a which satisfies Eq. (5 
If a’ is any a such that 


* py! dS = ‘ pdS = L 
St St 


verse is also true 


then a — a’ = ay since 
S (p — p')dS =0 
Thus (a, a-—-a’)= O or (a, a’) = (a, @ Using Eq (7 
(a, a)? = (a, a’)? < (a, @ (a’, a’ 
(a, a < (a’a’) 


Thus it is shown that of all admissible functions with total lift 
L, the angle of attack distribution satisfying Eq. (5) vields the 
least drag. 
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Addendum—“‘Effect of Diffusion Fields on the 
Laminar Boundary Layer’’ 


John W. Smith 
Aerodynamicist, Bell Aircraft Corporation, Buffalo, N.Y 
April 2, 1954 


[ SHOULD BE NOTED that the various forms of the diffusion 
equation given in reference 1 are descriptive of the diffusion 
field at low concentration of foreign gas. The general diffusion 
equation, * valid for binary gas mixtures of arbitrary strength, can 
be derived from the kinetic theory of Chapman and Cowling.’ 
One can start from the equation of conservation for foreign gas 
molecules (species 1) in a steady state, viz., 
~ 


div( mc) 


- 
where , is the foreign particle density and ¢; the mean molecular 


* Kindly brought to the author's attention by Dr. Jerome L. Fox of the 
Aerophysics Department, North American Aviation, Inc 





it fe 


base 
chat 


COTIC 


she wn 
rields 1 


yuTtIONS 


1q. (1 
»gonal 


- con 


il lift 


s the 


she 


ion 
ion 


ion 


ig.° 
gis 


the 





READERS' 


a 

velocity for species 1 The mean mass velocity c, the mean 
ee _ th —_> 

peculiar velocities C; and C:, and the diffusion current J, for 


species | are defined as follows: 


> > -_ 
pcy = pili + prt 
- - om 

( = ( ( 

—_ ‘s 

J, = ny 


where p refers to the total density, and p, refer to the partial den- 
sities. It is important to observe that the mean mass velocity 


is in general different from the mean molecular velocity for the 


mixture, and that the quantities « and v occurring in the boundary 
- 


layer equations are precisely the components of « From the 


fundamental equation of Chapman and Cowling? 


—_ _ n? j AN oN ( 
( © = D grad grad 77* 
nyn | { 
it follows that 
os Np» AN on | 
J), = - dD, = Ny — grad 7 
Nop | 7 f 


and the desired form of the diffusion equation is obtained when 
this result is substituted into the original equation of conserva 
tion. The form of the diffusion equation occurring in reference 
1 is a special case of this, applicable with the restriction noted 
above 

rhis restriction is removed through the following changes in 
the equations of the general theory (Section 3, reference 1 

1) The left side of Eqs. (3.6) and (3.8) should read 

(0/Ox){un,) + (O/O¥) (vn 

(2) In Eq. (3.8) and in the definition of @, lef. Eq. (3.11) 

replace n by pe / Neop. 


3) In Eq. (3.16), replace g2~'f[ smi] ' bv [p27 pert 


(4) In Eq. (3.21), replace ¢2~'g5’ in F; by [¢2~'os] ’. 

A numerical calculation of the sample problem for J/, = 8, 
based upon the general diffusion equation, has led only to small 
changes in the concentration profile, which tend to increase the 
concentration of helium present in the boundary layer 
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Design of Total Pressure Probes for Minimum 
Interference with Measured Flow 


J. L. Livesey 

Lecturer in Mechanical Engineering, University of Manchester, 
England 

April 9, 1954 


—ypesasongeed IN THESE COLUMNS Kaye, Brown, and Dieckmann! 
gave results of traverses of supersonic and subsonic streams 
in a '/,” bore pipe using extremely small total pressure tubes de- 
signed to have minimum interference with the measured flow 
The results compared extremely well with theoretical predictions 
but, for example, in Fig. 5, p. 204, in the case of turbulent sub- 
sonic stream, it appears that the tube reads consistently high 
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throughout the radial traverse. This tends to indicate an error 
i.e., the effective center of the probe 


rhis effect has 


of the “Displacement Type,” 
being displaced toward the higher total pressure 
often been noticed before in incompressible flow (see reference | 
and unpublished work in the United Kingdom indicates that the 
error is a function of the solidity of the cross section of the prob« 
tip and of the aspect ratio of this same cross section. So far work 
on compressible sub- and supersonic flow has not given any con 
sistent indication of the error. It will be noticed (see Fig. 1 of the 
subject paper) that the probes used by Kaye, Brown, and Dieck 
mann are of high solidity in the tip cross section and of fairly low 
aspect ratio. The solidity could be considerably reduced, eve 
in probes so small, by careful hand stoning of a chamfer on tht 
tips of the probes and without doubt the displacement error would 
be considerably reduced 

An estimate from Fig. 5 of the subject paper indicates that 
displacement of roughly 0.005 in. (as accurately as can be esti 
mated from the printed curve) would correct the traverse to the 
theoretical curve, for the case of subsonic turbulent flow. The 
work of Young and Maas? and also unpublished work of the 
writer in the United Kingdom, indicate that the displacement 
error is of the order of 0.15—).18 of the leading dimension of the 
probe tip, i.e., in this case 0.0045—).0054 in. which gives some 
weight to the assumption that the error is of displacement type 

The traverses also show a general overestimation of velocity 
level near the wall which is probably due to the effect of the in- 
duced velocities due to the image vortices, in the wall, of the 
vortices cast off from the tip of the probe. The writer has also 
observed this effect previously 

The writer apologizes for giving such attention to detail in a 
note, the major purpose of which was to establish the existence 
of extended laminar bounttary layers im stepersonic tube flow, 
but these traverses are probably nearer the theoretical pre 
dictions than their authors realized at first and the writer was 
interested to observe these displacement and wall effects in which 
he has been interested for some time. 

The extreme accuracy of the probe manufacture and of the 
traversing arrangements indicated by the subject paper indicate 
to the writer that a detailed examination, which would be ex- 
tremely valuable, of the displacement effect for high speed flows, 
could be carried out very rapidly in this apparatus 
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An Approximate Method for Determining the 
Effects of Wing Twist on Airplane Gust Loads 


Frederic M. Hoblit 
Structures Engineer, California Division 
Corporation, Burbank, Calif 


April 15,1954 


Lockheed Aircraft 


I REFERENCE 1, YEH AND MARTINEK develop a procedure for 
determining the effect of wing flexibility on airplane gust 
loads and give computed curves which show the alleviation due 
to twist, as a function of non-dimensional mass and flexibility 
parameters 

The purpose of this note is to show that a “common sens 


approximate procedure gives results that agree very closely with 
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FIG.| TIME HISTORY FOR d=35,6, cR=I5Q, 1=.487 
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those obtained by the more precise theory and hence may be 
useful for cases not covered in reference 1. 

Reference 1 is confined to the study of an idealized model in 
which a massless rigid wing is attached to the fuselage rigidly with 
respect to bending and elastically with respect to twist. The 
fuselage is free to plunge but not to pitch. Transient aerodynamic 
effects are considered, although certain approximations are intro- 
duced to facilitate the numerical work. Only the case of a 
sharp-edge gust istreated. The results are given as time histories 
of the ‘‘acceleration ratio,” 


Load factor increment 


Load factor increment as computed by the 
gust-load formula, pmSUwWmar,/2W 


for various values of a dimensionless flexibility parameter, cR, 
and of the familiar mass parameter, \ = 2//pcSm. In these 


expressions, m is the slope of the lift coefficient curve for the 
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rigid wing, UL’ is airplane forward speed, w is gust velocity, and 
p, S,c, , and W have the usual meanings, given in reference | 
The ‘“‘common sense”’ approximate method is as follows: 
(1) Obtain an “‘effective’’ slope of the lift coefficient curve, 
m,, defined as rate of change with respect to fuselage angle of 
attack of the airplane steady state lift coefficient, with the 
effect of wing twist included. In terms of the parameters 
used in reference 1, this is given by m, = m/(1 + cR/4)). 

(2) Obtain Agsyective = 2M /pcSm:. 

(3) From Fig. 2 or 6 of reference 1 or from similar information 
available elsewhere (for example, references 2-10), obtain the 
acceleration ratio for an airplane with rigid wing having the ef- 
fective value of \ obtained in Step 2. 

(4) Multiply the result of Step 3 by m,/m to give the accelera- 
tion ratio. 

Step 4 is necessary so that the denominator of the acceleration 
ratio is based upon m (rigid airplane) instead of m;, (effective 
value taking account of twist). 

Results of the approximate method are compared with some 
of the curves of reference 1 in Figs. 1-4. It is seen that the agree- 
ment is very close, especially for the larger values of A (35.6 
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| 71.2) typical of current aircraft. The wing flexibility is indi- 


ited on the curves by values of the ratio r = m,/m as well as 
by the parameter cR, since r seems to have a more direct physical 
significance 

It would appear that the approximate method could be used 
profitably and with considerable confidence in determining the 
fflects of wing twist in cases not covered by reference 1 These 
would include cases characterized by 

1) Gradient, triangular, or other gust profile 

2) Flow not two-dimensional 

3) Wing twist distributed along span 

1) Wing taper. 

It should be emphasized that both here and in reference 1 the 
inertia forces resulting from the variation with time of the de 
formation of the structure have been neglected 

Incidentally, a typographical error has been noticed in refer- 
) 


ence 1. In Eq. 23 there should be a coefficient 2 in front of the 


first integral on the right-hand side. 
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Effect of Roll on Dynamic Instability of 
Symmetric Missiles 


C. H. Murphy 

Ballistic Research Laboratories, Aberdeen Proving Ground, Md 
April 15, 1954 

~~ INTEREST HAS BEEN EXPRESSED recently in the dy- 

namic stability of symmetric missiles in rolling flight.’ * 

This note attempts to extend the discussion by stating a slightly 
neater form of generalized stability conditions* 5 and describing 
certain experimental results on dynamic instability. The nomen- 
clature employed here is based on that of reference 6. Some of 
the aerodynamic coefficients which we will use are defined below 


ind the remainder are easily identified 


Cx, = normal force coefficient slope (algebraic sign is nega- 
tive due to sign convention of coordinate system 

Up = drag coefficient 

Cu, = OCw/O(qd/2V) (d is body diameter) 

Cu, = oC (ad /2V) 

Cy, = Magnus moment due to angle of attack + (1/2) X 


pV2Sd (pd/2V)a 
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We will make use of two stability factors: 
] The ballistic stability factor 
kh 4tk,-2 (pd/V 


2(pSd/m) ¢ 


where &, and k, are axial and transverse radii of gyration in diam- 
eters, p the axial angular velocity, and m mass 
(2 The generalized dynamic stability factor 
sA) = P/O 
where 
AA) = Cra + 2Cp + 1/2)k.-? (Cx + Cura) + (2m/pSd)xr 
P(\) = 2(Cna + Cp — (1/2 k, 2 Cs. + (2m /pSd)xr 


and A denotes a level of damping. It is well known that if the 
independent variable is taken to be Z, distance measured along 
the trajectory in diameters, and the assumption of constant di- 
mensionless spin pd/V is made, the pitching and yawing motion 
may be written in the form 


B+ ia = Kyl Mt Z 4 Kiel Met ide")Z 


where );, ¢;’ are real constants and K; are complex functions of 


initial conditions 


A generalized stability condition may now be stated.5  , will 
be greater than J if and only if 

(1) Q(A) <0 

(2) Either s > {1/[5(2 — 5) when Cm, and §(2 — §) have 


the same sign or Cy, 0 when signs are different. \ has the 


character of a minimum permissible level of damping. If \ is zero 
the usual dynamic stability conditions result 
For this case of dynamic stability, it would seem certain that 


Q(0) <0. This follows from the fact that Cy _ is always negative 


and that Cy, + Cy, is usually taken to represent a damping 
moment and therefore be always negative. Surprisingly, positive 
values of Cu, + Cu, for simple cone cylinders’ with forward 
centers of mass at Mach Number of 0.8 have been measured to 
be as large as 5.6 and the resulting Q(0) has been positive 

Condition (2) for statically unstable missiles places a lower 
bound on spin when 0 < 3(0) < 2 and rejects the possibility of 
stabilization by spin when $ is outside this interval. For statically 
stable missiles it places no restriction on spin when § is inside this 
interval, but imposes an upper bound on spin when § is outside 
Thus we see that dynamic instability is insured for bodies of 
revolution and is possible for finned missiles when § is either nega- 
tive or greater than two. 

An examination of the definition of 5 reveals that the most likely 
for CM», 
This Magnus instability has been observed for both 
finned missiles* ° Table 1 the 
aerodynamic coefficients are listed for representative models. 
The finned missile will become dynamically unstable if (pb/2V) > 
0.11. 

§ can become greater than two when the contribution of the 


way for it to become negative would be , to become 


negativ e 


and bodies of revolution.’ In 


(b, the fin span, is 3d for this model.) 


damping moment Cy, + Car, in Q(0) is small with respect to that 


of the Magnus moment in P(0). Since Cu, + Car, is multiplied 


by k,~*, this would probably happen for long missiles where the 
transverse radius of gyration is large. Since the body of the 
finned missile of reference 8 is a cone cylinder with fineness ratio 
of ten, estimates of its aerodynamic coefficients were made and 


Note that a value of § = 2.3 resulted 


also listed in Table 1 

The effect that variation of center of mass location has on § is 
of interest. The following transformation relation, where x is 
the change in c.m. location in calibers, can be easily derived.* § 
(x is positive for forward motion of the c.m.) and starred quan- 
tities are for the new c.m. location.) 


. * 
C Mo - 


= 5 + xO 
Cm, xCyn 
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TABLE lt 


( Mo T 
Ml Va Cb Ma ( Mina Va 3 
5-cal. cone cyl 
inder 0.8 2.3 0.25 1.63 0.7 3.0 0.84 
Finned missile*® & | 14.0 0.63 20.4 44.1 250.0 7.16 
Body alone® BP 3.1 0.51 10.9 +3.3 150.0 2.2 


+ S, the reference area, is rd?*/4, and d, the characteristic length, is the 


body diameter Vis Mach Number 


ie : a 1c¢ rere. ; fer ee 
Cm T ( MG ( Mg 1 ( “Ma + Cre + C Vo 2¢ n) 


2x?Cy 
a 


CM pa CMa + XCN», 


For the body alone of reference 8, Cue + Cn _ = 15 and 
2.0. If x = 


Thus a rearward movement of the c.m. makes it possible 


Cyy. —1.2, § takes on the optimum value of 
unity 
to stabilize this model by spin. 
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Unpublished data on dynamic properties of the finned 


On the Transonic Drag of Accelerated Bodies 


John W. Miles 
Department of Engineering, University of California, Los Angeles 
April 21, 1954 


A LINEARIZED APPROXIMATION to the transonic drag on an 


accelerated body of revolution has been calculated by 
Cole.» * 


vious results of slender body theory* 4 lead to a rather more 


The purpose of the present note is to show that pre- 


general form of the result, paralleling the supersonic drag inte- 
grals of Ward® and von Karman.® 


Let ./ and \/p be the instantaneous and initial Mach Numbers 
and jf a dimensionless acceleration parameter defined by 
(acceleration): (body length 


M = = rs (1) 
(sonic velocity )? 


In practical applications .\f will be quite small (approximately 


10~* per g for a 3-ft. body at sea level). A consideration of the 
space-time domains of dependence of a uniformly accelerated 
body (see, e.g., references 1 and 2) leads naturally to the intro- 


duction of the acceleration parameter 
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for bodies of arbitrary cross sectional distribution. 
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p= (M—1!1 (2; M , , 9 


In particular, it is found that the transonic pressure distribution 


over an accelerated body is independent of the initial conditions jf 


Ma a Bw) < (yo? — 1)7? 3 


Using the results of references 3 and 4 (see especially Eqs 


(6)-(8) of reference 3) and assuming .V/ to be small, the drag on a 
slender body of revolution at zero incidence has been calculated 


in two sub-domains of Eq. (3), viz., 


(a) M>O, M<1 
19)9 ae me = 
pU?/? l 
D, = In x 
8x (Jo Jo (u2 + |x —€ + |p 
+] , 
S”’(x)S"(e)dtdx — 2S"(1 In ( ) S"(x)dx + 
J 0 l—-£z 
Sr ( 
S"2(1) In } 
(MV + 1) (247)'*S(01) JI 
b) M<0, M l 
» *1 ag | ° 
pur? § n° + t= et T pw 
Dy = — In : x 
Sia \. 0 0 t= $F 
"] 
- or a f= 2 + pb 
S”’(x)S"(e)dtdx — 2S'(1 in| 
J0 l—x 
Sr 
S"(x) dx + S’2(1) In 5 
We+ 1) (—247)'75(1) Jf 
where 
p = density of undisturbed air 
U = instantaneous flight velocity 
l = body length 
P = body cross section 
R = body radius /] 
x, £ = dimensionless body coordinates (= 0 at nose; 1 at tail 
The error term in these results is O( | 1/ if « = O(1) or 0(.1/) if 


uw>1. Inthe limits ~—~ *F © Eqs. (4) and (5) approach zero 
and the known, supersonic result,® respectively 

If, in addition to the usual slender body assumptions? (par- 
ticularly the restriction to a pointed nose), it is assumed (follow- 
ing Ward*) that either (i) S(1) = 0, so that the body is pointed 
at both ends, or (ii) the body is of general cylindrical form near 
x = 1, sothat S’(1) = 0, and the generators of the cylinder are 
parallel to the flight path (so that there is no drag due to lift) the 


results of Eqs. (4) and (5) simplify to 


979 od i 
pU? [ 
yy = — In ; x 
Sr Jo Jo Liu? Yr ie | Tr |e 
S”"(x)S"(é)dtdx 6 
ane 71 7) e ; 
pl? p? + —¢ + bu 
D, = se In xX 
Sr Jo 0 pian. ois 
S"(x)S”"(e)dtdx (7) 


valid 


In particular, 


Moreover, subject to the stated restrictions, Eqs. (6) and (7) are 


the sonic drag of such a body, obtained by setting « = O in Eqs. 
(6) and (7), is given by 


l 
D, = 1 Dek» VW 0 M = ] 8 
D, = Ky WV 0, V= 1 9 


where Dx is the drag of the same body at definitely supersonic 


speeds, as given by von Karman's result® 


*] vy 
U2 l a ‘ 
Dix = p In | S"(x)S"(é)didx 10) 
47 0 0 ee" 
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No extensive numerical applications of Eqs. (6) and (7) to 
specific bodies have been carried out, although comparative 
checks of Eqs. (5) and (8 


for the sonic drags of a decelerated cone and accelerated, parabolic 


have been made against Cole’s results 


rc body, respectively. Some idea of the variation of Eqs. (6) 


7) with J may be obtained by noting that Eq. (6) ap 


ind 
proaches zero as uw — — ©, while Eq. (7) approaches Eq. (10) as 


“> + 


Moreover, the derivatives of Eqs. (6) and (7) are 


both given by the positive definite expression 


oD, oD pul f[' [' S%x)S"(e)dtdx - 
Ou Ou Sr Jo > wer jz —s 


Accordingly, it appears that the drag of a slender, accelerated 
body of type (1) or (11) follows a symmetrical, hysteresis loop of the 
type sketched in Fig. 1, the solid portions of the curve being given 
by Eqs. (6) and (7). 

The calculation of the dotted portions of Fig. 1 is now under 
way, and it is hoped to present numerical results for representa 
tive bodies in the near future. 
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Compressible Fluid 


Thermal Loads on Joints 


B. E. Gatewood 

Mechanics Department, USAF Institute of Technology, 
Patterson AFB, Ohio 

April 23, 1954 


ws TWO BARS of different materials are bolted or riveted 
together and subjected to a temperature change, or two bars 


Wri ght- 


FORUM 645 


of same material are subjected to different temperature changes, 
What makes 
the problem serious is that, if the bolts do not deflect, the end 
bolt must take all the load while the others take 
lb, bolt or rivet 1 takes all the load 
wr bolts when their 


the bolts must take loads to hold the bars together 


none (in Fig 


This paper investigates 


the loads on the rivets deflection is con 
sidered. 

Consider two bars, bar 1 with area A;, modulus of elasticity 
k, temperature change 7), coefficient of thermal expansion a, 


and bar 2 with A», Eo, 7s, a Assume that the bars do not 
buckle, bow, or yi ld 
First, take the case of one loose fitting bolt, Fig. la, with gap 


2e. From the elongation and load equations 


aLT, + ((oL)/Ei| = azlhT., + [(o2L)/E2| + « l 


there results 


1 + [(4,2)) /(A2E 
From Eqs. (2) and (3), it is apparent that the larger e is and the 


smaller LZ is, the smaller the load ? In fact, for no load 


Eq. (4), with 7, = 0, 7; = 400°F, a = 12(10~) for aluminum 


alloy, L = 2 in., gives e = 0.0096 in., which is equivalent to the 
usual tolerances in ordinary bolts, 

Next, consider the distribution of load to several rivets, Fig 
lb Assume that each rivet deflects in proportion to the load it 
carries (the deflection may be due to bearing on bar, shear on 


rivet, bending on rivet, ete 


is deflection of rivet number », é. is allowable deflection 


is ratio of load on rivet to 


where ¢ 
of the rivets at allowable load P., 7 


total load P 1 equations can be written 


From Eq. (1), V 4+ 


af P rrecP 
a T; - = aol T 
1,F, Aolk aP 
P< Pic 
a7 > r, = al T > r T 
4,E) ; 1 Ask so ( 
e-P 
\ sP 


< a1 » a x 


and solve for r,/ry to get 
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Notation for joints. 
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ic. 2. Load P, on end rivet. 


aie + HY — — n = 1,2, 


IN IN rx 
sP ( ) ® 
1 + ; 
eA Ay A>»E» 


P can be written as 


H = 


From Eq. (6), 


| 
P =PiJ, J = 
, SV 
all 
asl» { (8S) 
5 = al 
Py =A Fiay,T 4 
7 4 A\EF, 
Ask, 


where 7, is the load without any rivet deflection. Further, the 
load on each rivet becomes 
r.P = rirJPo = qrPo, Qn = YrJ (9) 


» 


more load than rivet 2, ete., the pri 


mary value needed is g; for different values of N, s/a, and H, 
where // may be regarded as the joint stiffness parameter. Fig. 2 
shows a plot of g; against H with s/a and N as parameters 

The results shown by Fig. 2 indicate that for the usual values 
of H, the end rivet does not take all the load P)—the deflection 
of the rivets reduces the total load and divides it 
For example, take aluminum alloy bars with E; = Ey = 


Since the end rivet 1 carries; 


among all the 


rivets. 
107 psi, Ay = Ae = 0.50 in.?2, ands = 1.00 in., P. = 1500 Ib., 
e- = 0.01 in., whence 7 = 0.06, Hence, for s/a = 0.01 and N = 


2,q: = 0.46P, ete. For this case, Eq. (7) gives gz = 0.43Py. If 
é- approaches zero, then 7 becomes large and the curves show 
all the load 
result from the curves is that for four or more rivets, the load on 
and N, 


differences occurring in the region 0 < H < 0.10 


that the end rivet carries nearly An unexpected 


the end rivet is nearly independent of s/a some small 


Design of Tubes for Maximum Torsional 
Frequency 


P. Cicala 
Professor, Politecnico di Torino, Italy 
April 16, 1954 


Epes RAPID PROCEDURE for minimum weight design of a 
torsion tube to attain a prescribed torsional primary fre- 
quency, was indicated by E. P. MacDonough.' In this con- 
nection, it seems in order to present some results from a previous 
paper? which may provide more general information 


With reference to the equation for natural torsional vibrations 


w?) 1,0 l 


(GJ, torsional stiffness, 0’ gradient of twist angle 6, w frequency; 
all partial moments of 


GJ’ = 


the summation or integral extends to 
inertia /, from current section to wing tip) the minimum weight 
condition yields 

») 


= kV (dw/ds) /(OGJ,/ds) 2 


where & is a constant factor, w is the structural weight per unit 


span, and s the varied cover thickness 
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rhese equations, which agree with reference | in the particular 
case considered there, may lead to important corrections in the 
case of stiffened monocoque wing structures, where an increase 
reduction of stringer material 


f covering weight carries a 


Therefore, Ow/Os may show marked variations along the span 


In this case, when the right-hand side of Eq. (2) depends upon s, 
the solution can be obtained by iteration, starting from a pre 


distribution and successively calculating 06’, 


Similar 


sumed thickness 


4, GJ, |from Eq. (1)| and hence a new distribution for s 
considerations apply when determining optimum thickness dis 
tributions to meet aileron effectiveness requirements.” 

rhe bending stiffness of the spars was also taken into account 
in reference 2 under simplifying assumptions, viz., ribs and spar 
webs rigid, parallel spars. The structure of minimum weight 
presents an uncovered part near the root, whereas in the remain- 


The 


boundary section and the corresponding values of 6/6’ depend 


ing covered part Eq. (2) still holds, and gives 6’ = const 


1, where w, is the fundamental torsional 
If it is desired to find the 


upon the ratio w/w 
frequency for the uncovered structure 
least-weight sheet distribution giving a prescribed frequency for 
an overtone, more separate covered zones would appear along 
the span, with Eq. (2) holding with alternate signs for the square 
For instance, for the uniform section structure with cover 


root 
distribution giving the highest second frequency with a fixed 
sheet weight, only one outer zone is covered until 1 (w/w) 


1.8 (w: second natural frequency for bare structure). For higher 
w/w. values, another covered zone, starting at 17 per cent of the 
beam length from the root, extends to the first third part of the 
beam for wo/w — 0; the other zone tends to cover the remaining 
two-thirds, with maximum sheet thickness at its midpoint 
From these results it may be inferred that the reduction of 


cover thickness towards wing tip indicated by the first simplified 


The Second-Order Compressibility Rule for 
Airfoils 


Milton D. Van Dyke 

Aeronautical Research Scientist, National Advisory Committee for 
Aeronautics, Ames Aeronautical Laboratory, Moffett Field, 
Calif. 

May 3, 1954 

3 IS THE PURPOSE of this note to point out a remarkable rule 
that expresses the second-order subsonic pressure distribution 

over any thin two-dimensional airfoil directly in terms of the 

incompressible pressures. With this rule it is easy to reproduce 

the results of all the various special second-order solutions ob 

tained previously by laborious calculations. The rule can be 

deduced from an extension of transonic similitude to include all 

second-power terms that was initiated by Imai! and completed 

by Hayes.2. Hayes has shown that in either subsonic or super- 

sonic flow the ratio of the second-order pressure term to the first- 


order term at any point on the surface of an airfoil is proportional 
to 
T y¥+ 1 , 
Vi + 201 \/?) 
(1 — AJ?) 2 
where 7 is the thickness parameter. Now for subsonic flow the 


first-order pressure term is related to its value in incompressible 
flow by the Prandtl-Glauert rule. Combining these two results 
yields the second-order compressibility correction 

Thus, consider first an airfoil at zero angle of attack 
pose that the incompressible pressure coefficient on one surface is 


Sup- 


expanded in powers of the thickness ratio 7 to give 


rh(x) + r7g(x) +... (1 


Cp, = 


647 


FORUM 


analysis may be excessive, specially for engine-carrying wings, 
for which flutter modes more nearly approach an overtone of still 
air vibrations.‘ It should be noted that a uniform section tube 
with covering designed according to the proce lure of reference 
1, would be more than 20 per cent heavier than a constant thick 
ness cover giving the same second frequency Higher weight 


penalties might be encountered when the beam carries concen 


trated midspan masses. This makes evident the opportunity 
of limiting the cover thickness reduction to a small zone near 
the wing tip, when important masses are carried by the wing, 


unless a monocoque wing structure is used which allows greater 
cover thickness in the first part of the structure without important 


increase in total weight 
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* For the uniform section tube, if the cover thickness is halved from 0.6 


beam length to tip, a 12 per cent weight saving is obtained with respect to 


the constant thickness cover giving the same primary frequency, while a 17 


per cent excess weight is met with refpect to the constant thickness cover 
giving the same secondary frequency Halving the thickne from O85 
beam length to tip gives a 7 per cent and 3 per cent weight saving, respec 


tively Effects of spar stiffness are here disregarded 


Then at any subsonic Mach Number .J/ the pressure distribution 


is given by 


Cp, = Kirf(x + Kor*g(x) 4 2 
where 

. l l 
AK, = 

V1 VW 3 
, (y+ 1)M!+ 48 
K. = 

13 


If the airfoil is at angle of attack a, there will be a term linear 
in @ that follows the rule for the term in 7, and two second-order 
terms proportional to 7a and a? that follow the rule for the term 
It should be noted that these series break down near stag 
for 


in 7? 
nation points, and that they are believed to 
Mach the 


converge only 


subcritical Numbers, for which flow is everywhere 


subsonic 
Integration leads to corresponding rules for the lift and mo 


ment. Thus the incompressible lift coefficient has the form 
Cl = 2ra(l + Ar) 4 . (3 
where the constant \ can be determined from the profile. Then 


at subsonic speeds, to second order in thickness and angle of 
attack 

Cig = 2ra(K, + AeAr) 4 ; ! 
the Prandtl-Glauert 


the 


In these rules the first term represents 
correction. The factor A 
(aside from a factor of 2) just that arising in Busemann’s second 


Thus essentially the same rule 


associated with second term is 


order theory for supersonic flow 
holds 
flow it happens that the functions f(x) and g(x 


at both subsonic and supersonic speeds. In supersoni 


depend only upon 
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the local airfoil slope, whereas in subsonic flow they are given 
by integrals over the entire surface (see Lighthill’s paper*) which 
can, if necessary, be evaluated numerically 

As an example, consider the pressures on an ellipse of thickness 
ratio 7 at zero angle of attack. If the ellipse extends over x? < 1, 


the incompressible pressure distribution is given by 


27(1 2x?) 


It follows that in subsonic flow 
(y + 1) M4 + 4p? 


134 


(6) 


and this agrees with the results of analyses by Hantzsche and 
Wendt# and by Kaplan.® 

This example illustrates that the second-order thin-airfoil 
solution is not valid near a stagnation point, where even for in 
compressible flow it predicts infinite values of pressure coefficient 
In a subsequent paper the writer will show, as Lighthill did for 
incompressible flow,* how the subsonic solution can be rendered 
uniformly valid near stagnation points 

As a second example, consider the lift on a symmetrical Jou 


kowski airfoil. In incompressible flow the lift coefficient is® 


(y + 1) M4 + 46? 
7 
3°V 3,3 
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and this agrees with the result that Hantzsche and Wendt found 
by more elaborate computation 

The question naturally arises whether the rule for subsonic 
flow can be extended to still higher approximations using the 
coefficients from supersonic theory, which are known to fourth 
order.’ This question has been answered in the negative by con- 
sidering the known third-order solution for subsonic flow past a 


sinusoidal wall.*? 
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The Secondary Flow About Struts and Airfoils (Continued from page 608) 


and Eq. (7.9) becomes from Eqs. (A.10) and (A.18) 


6r n° , aT 2a 29 
ay = cosec- t= cot 
c(dU /dz) 2 2 n n 


we) G- 5)0-G)) +a 


r® 9 6 
oar 
14,175 \n 


4 nT “*) T (- 
cot sedi 
2 2 2 n 


For calculations when n — 2 this may be expanded in powers of [1 — (7/2) ]. 


yG- 


n 24 l ny 
— cot (A.19) 
rs n so k&2 


The first three terms are: 
_sJn\* n\? 
-7(*)'-0(*))4 
2 2 
a \* " ny n \" 
} + 50 — 84 (A.20) 
9 >) » 


The Effect of Tip Removal on the Natural Vibrations of Uniform Cantilevered Triangular Plates 


(Continued from page 633) 


quite well with both the numerical solutions of Barton* 
and the experimental work of Dalley and Ripper- 
ger.' 

The measured beam frequencies also agree closely 
with the theoretical values based on elementary beam 
theory.® 
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